[HNOI2014]江南乐

题目地址:https://www.luogu.com.cn/problem/P3235

题目意思是这样的:
给定(n)堆石子,每次操作选择一堆石子数目(x)((x>=f)),把它分成(m)堆((2<=m<=x)),分出来的每一堆中石子数目最大值和最小值相差不超过(1)(也就是均分的意思),问你先手是否必胜。

其实这个问题可以转化成(n)个子问题,每一个子问题都只对一堆石子进行操作。

然后就可以转化成(SG)函数的问题了。

SG定理:(SG=SG(a_1) igoplus SG(a_2) igoplus SG(a_3)... igoplus SG(a_n)),及总问题的SG值为分问题的SG值的异或和,本文就不进行证明了。

想法1

考虑到(x)(SG)值只跟(f)有关,题目在最开始已经把(f)告诉我们了,我们可以考虑预处理。
我们用枚举每个(x)分成的石子堆数(m)((2<=m<=x)),将这(m)堆石子的(SG)异或起来,然后对于所有的(m)的异或和取(mex)
这显然是(O(x^2))的做法,会直接TLE。

所以我们不能预处理,我们就想到了记忆化搜索。

想法2

受到想法1的启发,我们可以枚举(m),然后异或和取(mex),这不过这一次我们不预处理出(SG)值,而是边暴力递归边存(SG)值,遇到已经算过的(SG)值直接返回即可。
那么如何计算异或和呢?
我们用到异或的一个性质(x igoplus x=0)
看到下面这张图

我们发现把(x)分成(m)份的异或和只有可能是(0,SG(x/m),SG(x/m+1),SG(x/m) igoplus SG(x/m+1))
那么接下来该如何计算(mex)值呢
我们每一次算完异或和以后,用一个数字打上一个特殊的标记,然后最后统计(mex)值的时候找到第一个没有这种标记的值即可

然后我们就可以拿到了70分(TLE30)

#include<bits/stdc++.h>
using namespace std;
int T,F,sg[200000];
int n,a,po[200000],q=0;
int SG(int x,int pos)
{
	if(sg[x]>-1)return q--,sg[x];//已经算出来过
	if(x<F)return q--,0;//不可再分
	sg[x]=0;
	for(int i=2;i<=x;i++)//计算异或和
	{
		int x1=x/i,x2=x/i+1,n2=x-x1*i,n1=i-n2;//计算分出来的两个值和两个值所占的数量
		if(n1%2==0&&n2%2==1)po[SG(x2,++q)]=pos;
		if(n1%2==0&&n2%2==0)po[0]=pos;
		if(n1%2==1&&n2%2==1)po[SG(x2,++q)^SG(x1,++q)]=pos;
		if(n1%2==1&&n2%2==0)po[SG(x1,++q)]=pos;
	}
	int pp=0;
	while(po[pp]==pos)pp++;//计算mex值
	return sg[x]=pp;
}
int main()
{
	scanf("%d%d",&T,&F);
	if(F==1)sg[1]=0;
	memset(sg,-1,sizeof(sg));
	while(T--)
	{
		memset(po,0,sizeof(po));
		scanf("%d",&n);
		int res=0;//res是总的SG值
		for(int i=1;i<=n;i++)
		{
			scanf("%d",&a);
			res^=SG(a,++q);
		}
		if(res==0)printf("0 ");
		else printf("1 ");
	}
	
	
	return 0;
}

想法3

我们可以观察到其实同一个(x)分成(m)份时会出现(t)(t+1)两个值,我们可以把分出来的(t)(t+1)都相同的(m)一起算,此时答案只有可能是(SG(t),SG(t) igoplus SG(t+1))两种,而且存在一定的奇偶关系

我们可以发现,同一个t和t+1,m0与m0+2的答案是一样的,因此我们可以枚举t,然后算出对应的前两个m即可

可以看下面这个图,尤其是1和2的蓝色部分,规律非常明显

AC Code

#include<bits/stdc++.h>
using namespace std;
int T,F,sg[200000];
int n,a,po[200000],q=0;
int SG(int x,int pos)
{
	if(sg[x]>-1)return q--,sg[x];
	if(x<F)return q--,0;
	sg[x]=0;
	int i=2;
	while(i<=x)
	{
		if((x/i)==(x/(i-1))&&(x/i)==(x/(i-2)))i=(x/(x/i));//算完前两个后跳过
		int x1=x/i,x2=x/i+1,n2=x-x1*i,n1=i-n2;
		if(n1%2==0&&n2%2==1)po[SG(x2,++q)]=pos;
		if(n1%2==0&&n2%2==0)po[0]=pos;
		if(n1%2==1&&n2%2==1)po[SG(x2,++q)^SG(x1,++q)]=pos;
		if(n1%2==1&&n2%2==0)po[SG(x1,++q)]=pos;
		i++;
	}
	int pp=0;
	while(po[pp]==pos)pp++;
	return sg[x]=pp;
}
int main()
{
	scanf("%d%d",&T,&F);
	if(F==1)sg[1]=0;
	memset(sg,-1,sizeof(sg));
	while(T--)
	{
		memset(po,0,sizeof(po));
		scanf("%d",&n);
		int res=0;
		for(int i=1;i<=n;i++)
		{
			scanf("%d",&a);
			res^=SG(a,++q);
		}
		if(res==0)printf("0 ");
		else printf("1 ");
	}
	
	
	return 0;
}
原文地址:https://www.cnblogs.com/Laoli-2020/p/14285390.html