标签：ch frac 选数 复杂度 状压 times HNOI2012 while isdigit

题意：

题面

分析：

我们发现每个数 \(n\) 是否被选，只与 \(\frac{n}{3},\frac{n}{2},2n,3n\) 有关，那么我们考虑建一张表，表上每一行按照 \(\times 3\) 的方式递增，每一列按照 \(\times 2\) 的方式递增，那么对于同一张表，任意上下左右相邻的数都是不能选的，那么这样的表一共有 \(n-\frac{n}{2}-\frac{n}{3}+\frac{n}{2\times 3}\) 种，每一张表的转移是 \(\log_2 n\times 2^{2\times \log_3 n}\) 虽然理论复杂度达到了 \(5e10\) 但是实际跑下来十分优秀，可能因为后来的表变的极小，转移的复杂度大幅降低

代码:

#include<bits/stdc++.h>

using namespace std;

namespace zzc
{
	int read()
	{
		int x=0,f=1;
		char ch=getchar();
		while(!isdigit(ch))
		{
			if(ch=='-') f=-1;
			ch=getchar();
		}
		while(isdigit(ch))
		{
			x=(x<<1)+(x<<3)+ch-48;
			ch=getchar();
		}
		return x*f;
	}
	
	const int maxn = 1e5+5;
	const int maxm = 5005;
	const int mod = 1e9+1;
	int n,r,sta[maxm],num[20],f[20][maxm];
	long long ans=1,sum;
	bool vis[maxn];
	
	void update(int x)
	{
		r=1;sum=0;
		memset(num,0,sizeof(num));
		memset(f,0,sizeof(f));
		for(int i=x;i<=n;i<<=1,r++) for(int j=i;j<=n;j*=3) num[r]++,vis[j]=true;
		for(int i=0;i<(1<<num[1]);i++) f[1][i]=sta[i];
		for(int i=2;i<r;i++)
		{
			for(int j=0;j<(1<<num[i-1]);j++)
			{
				for(int k=0;k<(1<<num[i]);k++)
				{
					if(sta[j]&&sta[k]&&!(j&k)) f[i][k]=(f[i][k]+f[i-1][j])%mod;
				}
			}
		}
		for(int i=0;i<(1<<num[r-1]);i++) sum=(sum+f[r-1][i])%mod;
		ans=ans*sum%mod;
	}
	
	void work()
	{
		n=read();
		for(int i=0;i<=(2<<11);i++) sta[i]=((i<<1)&i)?0:1;
		for(int i=1;i<=n;i++) if(!vis[i]) update(i);
		printf("%lld\n",ans);
	}

}

int main()
{
	zzc::work();
	return 0;
}

标签：ch,frac,选数,复杂度,状压,times,HNOI2012,while,isdigit
来源： https://www.cnblogs.com/youth518/p/13900334.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。