bzoj3576[HNOI2014]江南乐
首先整个局面的SG值等于各个独立子局面的SG值异或和,于是只要求SG(1~100000)了。考虑将一堆i个石头分成j堆时,它的后继状态是 {\((j-i \ mod\ j)\) 个 \((i/j)\),\((i\ mod \ j)\) 个 \((i/j+1)\)},这个后继状态的SG值同样等于这j个局面的SG异或和,于是要计算SG(i)就只需枚举将它分成j=1~i堆(j>i是没有意义了的,在j=i时已经考虑到了),然后把后继状态的SG取mex即可。
于是就有\(O(n^2)\)的算法:
for(int i=1;i
然后容易想到这个可以分块优化,对于一段\(i/j\)相同的j,\(j-i\ mod \ j\) 和 \(i\ mod \ j\)有一个是公差为偶数的等差数列,一个是公差为奇数的等差数列,所以只需要考虑一段\(i/j\)相同的j的其中两个即可。
#include
#include
#include
#include
#include
#include
#include
#define P puts("lala")
#define cp cerr<<"lala"<'9') {if(ch=='-')g=-1;ch=getchar();}
while(ch<='9'&&ch>='0') re=(re<<1)+(re<<3)+(ch^48),ch=getchar();
return re*g;
}
typedef long long ll;
typedef pair pii;
const int N=125;
const int M=100050;
int a[N],F,n;
int SG[M],stk[M],top=0;
bool vis[M];
int main()
{
#ifndef ONLINE_JUDGE
freopen("1.in","r",stdin);freopen("1.out","w",stdout);
#endif
int T=read(); F=read();
for(int i=1;i