题目链接

https://atcoder.jp/contests/agc024/tasks/agc024_f

题解

考虑对每个串都求出它是几个给定的串的子序列。
在 \(O(4^n)\) 暴力的基础上考虑优化，状态可以进行如下的压缩：用有序 01 字符串对 \((S,T)\) 代表当前已有字符串 \(S\), 后面还需要加字符串 \(T\). 转移就枚举下一个选的是什么即可。
可以发现转移是一个DAG，并且每条路径对应的字符串不同。因此直接在DAG上进行路径计数即可。
时间复杂度 \(O(2^nn)\).

代码

#include<bits/stdc++.h>
#define llong long long
#define mkpr make_pair
#define riterator reverse_iterator
using namespace std;

inline int read()
{
	int x = 0,f = 1; char ch = getchar();
	for(;!isdigit(ch);ch=getchar()) {if(ch=='-') f = -1;}
	for(; isdigit(ch);ch=getchar()) {x = x*10+ch-48;}
	return x*f;
}

const int N = 20;
char a[(1<<N+1)+3];
int f[(1<<N+1)+3][N+3];
int bitcnt[(1<<N+1)+3],lg2[(1<<N+1)+3];
int nxt[(1<<N+1)+3][2],cnt[(1<<N+1)+3];
int n,m;

int main()
{
	for(int i=1; i<(1<<N+1); i++) bitcnt[i] = bitcnt[i>>1]+(i&1);
	lg2[1] = 0; for(int i=2; i<(1<<N+1); i++) lg2[i] = lg2[i>>1]+1;
	scanf("%d%d",&n,&m);
	for(int i=0; i<=n; i++)
	{
		scanf("%s",a); for(int j=0; j<(1<<i); j++) f[j|(1<<i)][0] = a[j]-48;
	}
	for(int i=1; i<(1<<n+1); i++)
	{
		int len = lg2[i]; int sta = i^(1<<len);
		for(int x=0; x<2; x++)
		{
			nxt[i][x] = -1;
			for(int j=len-1; j>=0; j--)
			{
				if(((i>>j)&1)==x) {nxt[i][x] = (i&((1<<j)-1))|(1<<j); break;}
			}
//			printf("nxt[%d][%d]=%d\n",i,x,nxt[i][x]);
		}
	}
	for(int j=0; j<=n; j++)
	{
		for(int i=(1<<j); i<(1<<n+1); i++)
		{
			int len = lg2[i]; int sta = i^(1<<len);
			if(f[i][j]==0) continue;
			cnt[i>>len-j] += f[i][j];
			for(int x=0; x<2; x++)
			{
				if(nxt[(sta&((1<<len-j)-1))|(1<<len-j)][x]==-1) continue;
				int sta2 = nxt[(sta&((1<<len-j)-1))|(1<<len-j)][x],len2 = lg2[sta2]; sta2 ^= (1<<len2);
				sta2 |= (((i>>len-j)<<1|x)<<len2);
//				printf("(%d,%d)+%d->(%d,%d)\n",i,j,x,sta2,j+1);
				f[sta2][j+1] += f[i][j];
			}
		}
	}
	int ans = 0;
	for(int i=1; i<(1<<n+1); i++)
	{
//		printf("cnt[%d]=%d\n",i,cnt[i]);
		if(cnt[i]>=m)
		{
			if(lg2[i]>lg2[ans]||(lg2[i]==lg2[ans]&&i<ans)) {ans = i;}
		}
	}
	int len = lg2[ans]; for(int i=len-1; i>=0; i--) printf("%d",(ans>>i)&1); puts("");
	return 0;
}