51nod B 异或约束和

B 异或约束和

一 题目

定义 f(i) 为 i 的所有约数的异或和，给定 n(1≤n≤10^14) ，求 f(1) xor f(2) xor f(3) xor…xor f(n) （其中xor表示按位异或）

样例解释：

f(1) = 1

f(2) = 1 xor 2 = 3

f(3) = 1 xor 3 = 2

f(4) = 1 xor 2 xor 4 = 7

1 xor 3 xor 2 xor 7 = 7

输入
一行，输入一个整数n
输出
一行，一个整数为答案

输入样例
4
输出样例
7

二 思路

首先，简单了解一下异或：
异或是一种基于二进制的位运算，用符号XOR或者 ^ 表示，其运算法则是对 xor 两侧数的每一个二进制位，同值取0，异值取1。

它与布尔运算的区别在于，当 xor 两侧均为1时，布尔运算的结果为1，异或运算的结果为0。
简单理解就是不进位加法，如1+1=0，,0+0=0,1+0=1。

其次针对于这一道题，看到了数据那么大，肯定不可能是O(n^2)或者是O（n）的算法，肯定直接就卡了。说明肯定是有规律的

然后考虑到枚举因数x，然后算出他在n中出现的次数,即⌊n/x⌋,然后根据算出次数的奇偶性，判断是否要异或x

但是，这样复杂度是O(n)的，看到[n/x]很容易想到数论分块，反正我没想到 。

最后问题就变成了考虑如何快速查询一个连续区间的异或和。

针对于一个区间 [l,r] 的异或和,可以变换成[1,l-1] xor[1,r] 。【重点】

那么问题再一次转换成了求【1，x】的某一个区间的异或和。
这个时候我们就需要冷静思考 暴力打表

打表结果如下：

然后我们就可以很清晰的看到有规律 【重点】
$$f(0,n)=f(1,n)\{\begin{array}{l}nn\%4=0\\ 1n\%4=1\\ n+1n\%4=2\\ 0n\%4=3\end{array}$$
然后就可以根据这个规律，求解。
注：如果想要详细的证明的话，请移至 https://www.cnblogs.com/My-snowing/p/10439231.html

三 代码

打表代码：

#include<bits/stdc++.h>
using namespace std;
long long n;
long long sum=0,s=0,tot; 

int main()
{
	cin>>n;
	for(int i=1;i<=n;++i)
	{
		sum=sum^i;
		cout<<i<<' '<<sum<<endl;
	}
	
	return 0;
}

程序代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n;
ll sum=0;

ll f(ll x)//找到1~x的异或和 
{
	if(x%4==1)	return 1;
	else 	if(x%4==2)	return x+1;
	else	if(x%4==3)	return 0;
	return x;
}

int main()
{
	cin>>n;
	for(ll l=1,r;l<=n;l=r+1)
	{
		r=n/(n/l);
		if((n/l)&1)//判断奇偶
		sum^=f(l-1)^f(r);
	}
	cout<<sum<<endl;
	return 0;
}

完……