【模拟】XMU 1599 斐波那契汤
题目链接：
http://acm.xmu.edu.cn/JudgeOnline/problem.php?id=1599
题目大意：
给k,m,q以及f[1]...f[k],当n<m时,f[n]=f[1]/2+f[2]/2...f[n-1]/2,
n>=m时 F(n)=F(n-1) XOR F(n-2) XOR……XOR F(n-m); n>m.
求F(L) xor F(L+1) xor …… xor F(R)。
(1 =< k <=m <=1e5,且m-k<64；q<=10000，0< L <= R <=1e18,f[i]在int64范围内）
题目思路：
【模拟】
m-k<64，前面暴力做到m，之后发现f是一个以m+1为循环节的循环函数。所以就可以做了。
1 // 2 //by coolxxx 3 // 4 #include<iostream> 5 #include<algorithm> 6 #include<string> 7 #include<iomanip> 8 #include<memory.h> 9 #include<time.h> 10 #include<stdio.h> 11 #include<stdlib.h> 12 #include<string.h> 13 #include<stdbool.h> 14 #include<math.h> 15 #define min(a,b) ((a)<(b)?(a):(b)) 16 #define max(a,b) ((a)>(b)?(a):(b)) 17 #define abs(a) ((a)>0?(a):(-(a))) 18 #define lowbit(a) (a&(-a)) 19 #define sqr(a) ((a)*(a)) 20 #define swap(a,b) ((a)^=(b),(b)^=(a),(a)^=(b)) 21 #define eps 1e-8 22 #define J 10 23 #define MAX 0x7f7f7f7f 24 #define PI 3.1415926535897 25 #define inf 10000000 26 #define N 100004 27 using namespace std; 28 int n,m,lll,ans,cas; 29 long long l,r,f; 30 long long a[N],b[N]; 31 int main() 32 { 33 #ifndef ONLINE_JUDGE 34 // freopen("1.txt","r",stdin); 35 // freopen("2.txt","w",stdout); 36 #endif 37 int i,j,k; 38 // while(~scanf("%s",s1)) 39 while(~scanf("%d",&n)) 40 // for(scanf("%d",&cas),l=1;l<=cas;l++) 41 { 42 memset(a,0,sizeof(a)); 43 scanf("%d%d",&m,&cas); 44 for(i=1;i<=n;i++) 45 scanf("%lld",&a[i]); 46 for(i=n+1;i<=m;i++) 47 { 48 for(j=1;j<i;j++) 49 { 50 a[i]+=a[j]/2; 51 a[j]=(a[j]+1)/2; 52 } 53 } 54 for(i=1;i<=m;i++) 55 a[m+1]^=a[i]; 56 m++; 57 for(i=1;i<=m;i++) 58 b[i]=b[i-1]^a[i]; 59 for(i=1;i<=cas;i++) 60 { 61 scanf("%lld%lld",&l,&r); 62 f=(r-l)%(m+m);l=(l-1)%m+1; 63 if(l+f<=m)printf("%lld\n",b[l+f]^b[l-1]); 64 else if(l+f<=m+m)printf("%lld\n",b[m]^b[l-1]^b[l+f-m]); 65 else printf("%lld\n",b[l+f-m-m]^b[l-1]); 66 } 67 } 68 return 0; 69 } 70 71 /* 72 // 73 74 // 75 */
