暑假没过,今天打了下表发现了规律 ,然后就A
然后for 一边把Q^=Pi
然后我把前面的打表出来发现竟然有规律
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
int main(){
int n;
int a[100][100];
scanf("%d",&n);
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
a[i][j]=i%j;
}
}
for(int i=1;i<n;i++){
for(int j=0;j<=n;j++){
printf("%d ",a[i][j]);
}
puts("");
}
}
/*
10
0 0 1 1 1 1 1 1 1 1 1
0 0 0 2 2 2 2 2 2 2 2
0 0 1 0 3 3 3 3 3 3 3
0 0 0 1 0 4 4 4 4 4 4
0 0 1 2 1 0 5 5 5 5 5
0 0 0 0 2 1 0 6 6 6 6
0 0 1 1 3 2 1 0 7 7 7
0 0 0 2 0 3 2 1 0 8 8
0 0 1 0 1 4 3 2 1 0 9
*/
发现从第三列每次都有循环节,第零列是0(0)第一列是1(0),第二列是2(1 0),第三列是3(1 2 0)。。。。。
就可按照循环节算出前面的(for i=1->n)(i mod 1)^...(i mod n);,因为他就是这个矩阵
。。balabala。。。
代码如下
:
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
#define N 1100000
ll ans;
int n;
ll a[N];
ll mod[N];
int num[N];
void init(){
mod[0]=0;
for(int i=1;i<N;i++){
mod[i]=mod[i-1]^i;
}
}
void work(int x){
memset(num,0,sizeof(num));
if(x==0||x==1) return ;
for(int i=2;i<=x;i++){
num[i-1]+=x/i;
num[x%i]++;
}
}
int main(){
init();
ans=0;
scanf("%d",&n);
work(n);
for(int i=1;i<n;i++){
if(num[i]&1){
ans^=mod[i];
}
}
ll b;
while(n--){
scanf("%I64d",&b);
ans^=b;
}
printf("%I64d\n",ans);
}