1. 程式人生 > >【hdu6275】【2017ccpc杭州L】Mod, Xor and Everything 題解

【hdu6275】【2017ccpc杭州L】Mod, Xor and Everything 題解

題目大意

       ~~~~~~ ( n  

m o d   1 )   x o r
  ( n   m o d   2 )
  x o r   . . .   x o r   ( n   m o d   n ) (n~mod~1)~xor~(n~mod~2)~xor~...~xor~(n~mod~n)
       n < = 1 e 12 ~~~~~~n<=1e12

題解

       ~~~~~~ 此題跟超級綿羊異或異曲同工。

       ~~~~~~ 求異或和有一個思路是:求出每一位分別是什麼。
       ~~~~~~ 從低到高第 k k 位(從0開始)的答案其實就是
i = 1 n n   m o d   i 2 k ( m o d   2 ) \sum_{i=1}^n\lfloor\frac{n~mod~i}{2^k}\rfloor(mod~2)
= i = 1 n n n i i 2 k ( m o d   2 ) =\sum_{i=1}^n\lfloor\frac{n-\lfloor\frac{n}{i}\rfloor i}{2^k}\rfloor(mod~2)

       ~~~~~~ 注意到 n i \lfloor\frac{n}{i}\rfloor 是可分塊的,於是把它分塊,這就成類歐啦~

       ~~~~~~ 當然可能還需要卡卡常:比如調整迴圈順序以剪枝,比如前 1e7 個直接暴力。。。

程式碼

#include<bits/stdc++.h>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;

typedef long long LL;

LL n;

bool f(LL a,LL b,LL c,LL n)
{
	if (!a) return (((n+1)&(b/c))&1)>0;
	if (a>=c || b>=c)
	{
		LL sqr=(n&1) ?(n+1)/2*n :n/2*(n+1) ;
		return ((f(a%c,b%c,c,n)+(a/c)*sqr+(n+1)*(b/c))&1)>0;
	} else
	{
		LL m=(a*n+b)/c;
		return (((m*n)^f(c,c-b-1,a,m-1))&1)>0;
	}
}

int T;
int main()
{
	scanf("%d",&T);
	while (T--)
	{
		scanf("%lld",&n);
		
		LL ans=0,sqrtn=min(30000000ll,n);
		fo(i,1,sqrtn) ans^=n%i;
		for(LL i=sqrtn+1, j; i<=n; i=j+1)
		{
			j=n/(n/i);
			LL lim=n/i*(j-i)+n%j, ans1=0;
			for(LL k=1; k<=lim; k<<=1) ans1+=f(n/i,n%j,k,j-i)*k;
			ans^=ans1;
		}
		
		printf("%lld\n",ans);
	}
}