[关闭]
@ysner 2018-07-31T09:22:27.000000Z 字数 1989 阅读 1679

[HAOI2015]数组游戏

博弈论 分块


题面

有一个长度为的数组,每个格子颜色为黑或白。两人轮流操作。每次操作选择一个白格,假设它的下标为。接着,选择一个大小在之间的整数,然后将下标为、...、的格子都进行颜色翻转。不能操作的人输。问先手是否有必胜策略。

解析

算法

状态压缩+记忆化搜索

算法

该问题很像翻硬币问题。
既然不能翻黑色格子,说明黑色格子的值对答案没有直接影响。
由性质
可知整个游戏的值为所有可行操作的异或和,即所有白格的值异或和。

值,就是对 所有后继状态(操作后状态)的值 的异或和取
再应用一下上面那个性质,可得一个递推式

模拟一发即可。(看不懂上面可以去蹭蹭我的博弈论总结
然而我因没注意到存在和一定有得怀疑人生。。。
复杂度为

  1. const int mod=1e9+7,N=1e7+100;
  2. int n,q,w,SG[N],ans,viss[N];
  3. il void getSG(re int x)
  4. {
  5. re int t=1,tot=0;
  6. viss[0]=x;
  7. while(233)
  8. {
  9. ++t;
  10. if(x*t>n) break;
  11. viss[SG[x*t]^tot]=x;
  12. tot^=SG[x*t];
  13. }
  14. re int tmp=0;while(viss[tmp]==x) ++tmp;
  15. SG[x]=tmp;
  16. }
  17. int main()
  18. {
  19. n=gi();q=gi();
  20. fq(i,n,1) getSG(i);
  21. while(q--)
  22. {
  23. ans=0;
  24. w=gi();
  25. fp(i,1,w) ans^=SG[gi()];
  26. puts(ans?"Yes":"No");
  27. }
  28. return 0;
  29. }

算法

经过愉快地打函数表后,发现只要,则
于是把值相同的值合为一块来计算即可。
注意存值的技巧:分为两半

  1. #include<iostream>
  2. #include<cstdio>
  3. #include<cstdlib>
  4. #include<cstring>
  5. #include<cmath>
  6. #include<algorithm>
  7. #include<vector>
  8. #define re register
  9. #define il inline
  10. #define ll long long
  11. #define max(a,b) ((a)>(b)?(a):(b))
  12. #define min(a,b) ((a)<(b)?(a):(b))
  13. #define fp(i,a,b) for(re int i=a;i<=b;i++)
  14. #define fq(i,a,b) for(re int i=a;i>=b;i--)
  15. using namespace std;
  16. const int mod=1e9+7,N=1e6+100;
  17. int n,q,w,SG[2][N],ans,blk,b[N],id,vis[N];
  18. il ll gi()
  19. {
  20. re ll x=0,t=1;
  21. re char ch=getchar();
  22. while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
  23. if(ch=='-') t=-1,ch=getchar();
  24. while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
  25. return x*t;
  26. }
  27. il int getSG(re int x)
  28. {
  29. x=n/(n/x);
  30. if(x<=blk) return SG[0][x];else return SG[1][n/x];
  31. }
  32. il void Pre()
  33. {
  34. for(re int i=1;i<=n;i=n/(n/i)+1) b[++b[0]]=n/(n/i);
  35. fq(i,b[0],1)
  36. {
  37. re int x=b[i],now=0;++id;vis[0]=id;
  38. for(re int j=x+x;j<=n;)
  39. {
  40. re int t=(n/(n/j))/x*x,tmp=(t-j)/x+1;
  41. vis[now^getSG(j)]=id;
  42. if(tmp&1) now^=getSG(j);
  43. j=t+x;
  44. }
  45. re int tmp=0;while(vis[tmp]==id) ++tmp;
  46. if(x<=blk) SG[0][x]=tmp;else SG[1][n/x]=tmp;
  47. }
  48. }
  49. int main()
  50. {
  51. n=gi();q=gi();blk=sqrt(n)+1;
  52. Pre();
  53. while(q--)
  54. {
  55. ans=0;w=gi();
  56. fp(i,1,w) ans^=getSG(gi());
  57. puts(ans?"Yes":"No");
  58. }
  59. return 0;
  60. }
添加新批注
在作者公开此批注前,只有你和作者可见。
回复批注