【模拟】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 */