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

转载于:https://www.cnblogs.com/thkkk/p/8677671.html