【2022国赛模拟】书——状压DP、容斥

最新推荐文章于 2023-07-25 22:46:31 发布

偶耶XJX

最新推荐文章于 2023-07-25 22:46:31 发布

阅读量142

点赞数

分类专栏：信息竞赛解题文章标签：算法动态规划状态压缩容斥原理

本文链接：https://blog.csdn.net/weixin_43960287/article/details/125126680

版权

信息竞赛解题专栏收录该内容

195 篇文章 5 订阅

订阅专栏

这个人知道原题链接

题目描述

在这里插入图片描述

题解

首先我们可以把 $1\sim n$ 剖分成若干条 $x,2x,4x,...2^yx$ 的极长链，那么对于每一段连续的‘1’，一定会占用某条链中间的一段，然后把链断成两条短链。

所以我们有个朴素的想法，考虑记录每一种长度的链的数量，然后按照这个状态转移。由于链长最多为6以及一些其他的限制，合法的状态数非常少，直接 map 记录合法状态然后暴力转移就能跑出来。

但是只是这样还不行，因为就算保证了每一个‘1’的限制，但是我们这种DP方式没有考虑‘0’，会导致某些不该出现2倍关系的位置出现了2倍关系。这个很好办，容斥一下即可。

代码

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define lll __int128
#define uns unsigned
#define fi first
#define se second
#define IF (it->fi)
#define IS (it->se)
#define END putchar('\n')
#define lowbit(x) ((x)&-(x))
#define inline jzmyyds
using namespace std;
const int MAXN=233;
const ll INF=1e18;
ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[50],lpt;
void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt>0)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}
const ll MOD=1e9+7;
int n,k,lm[6];
ll o[45][45],ans,g,mi[41];
char a[233];
unordered_map<ll,ll>dp[41];
void AD(ll&a,ll b){a+=b;if(a>=MOD)a-=MOD;}
int siz(ll s,int x){return s/mi[x-1]%g;}
ll cad(ll s,int x,int d){return s+d*mi[x-1];}
int main()
{
	freopen("book.in","r",stdin);
	freopen("book.out","w",stdout);
	n=read(),*new(int)=scanf("%s",a+1);
	mi[0]=1,g=n+1;
	for(int i=1;i<=6;i++)mi[i]=mi[i-1]*g%MOD;
	ll stg=0;
	for(int i=1;i<=n;i++)if(i&1){
		int y=1;
		while((i<<y)<=n)y++;
		stg=cad(stg,y,1);
	}
	dp[0][stg]=1;
	for(int i=1;i<=n;i++){
		if(i-7>=0)swap(dp[i],dp[i-7]),dp[i].clear();//为了省空间，毕竟unordered_map
		if(a[i]=='1')continue;
		for(int j=i,b=0;j>0&&j>i-6;j--,b+=(a[j]=='0')){
			int x=i-j+1;
			if(a[j-1]=='1')continue;
			for(auto it=dp[j-1].begin(),ED=dp[j-1].end();it!=ED;it++){
				ll s=IF;
				for(int k=x;k<=6;k++){
					int ct=siz(s,k);
					if(!ct)continue;
					const ll d=IS*ct%MOD*(x>1?2:1)%MOD*((b&1)?MOD-1:1)%MOD;
					for(int l=0;l<=k-x;l++){
						ll to=cad(s,k,-1);
						if(l>0)to=cad(to,l,1);
						if(k-x-l>0)to=cad(to,k-x-l,1);
						AD(dp[i][to],d);
					}
				}
			}
		}
	}
	print(dp[n][0]);
	return 0;
}