AT_abc343_G [ABC343G] Compress Strings 题解

分析

水题,评分能有 2100 可能是因为很多人卡 E 了。

我说真的,E 好难啊。

n 只有 20,考虑从状压的角度入手。定义状态函数 fs,i 表示当某个字符串 T 包含了所有 s 的二进制中为 1 的下标 SjT 末尾包含的子串为 SiT 的最小长度。那很显然的就有转移方程:fs2j,j=min(fs2j,fs,i+|Sj|Δ)。其中 Δ 表示 SjSi 拼起来时,Si 的后缀与 Sj 的前缀的最大匹配长度。

有个问题,当 S1=aab,S2=bc,S3=abc 时,上面的 DP 得到的答案是 7,即 T=aabcabc。但是很显然在 S1S2 拼起来之后,S3 与其的 Δ3。此时的 S2 已经被 S3 完全包含了,就相当于是 S1S3 拼起来的情况。所以当 SiSj(ji) 包含时,Si 对答案的代价为 0(可以不管 Si)。

任意 Si,SjΔ 可以暴力哈希,判断 Si 是否被 Sj 包含也可以暴力哈希。这两个预处理出来之后跑状压就行了。

复杂度 O(n22n),实测删掉被包含的之后飞快。

代码

#include<bits/stdc++.h>
using namespace std;
#define int long long
#define re register
#define il inline
#define pii pair<int,string>
#define x first
#define y second
#define gc getchar()
#define rd read()
#define debug() puts("------------")

namespace yzqwq{
	il int read(){
		int x=0,f=1;char ch=gc;
		while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc;}
		while(ch>='0'&&ch<='9') x=(x<<1)+(x<<3)+(ch^48),ch=gc;
		return x*f;
	}
	il int qmi(int a,int b,int p){
		int ans=1;
		while(b){
			if(b&1) ans=ans*a%p;
			a=a*a%p,b>>=1;
		}
		return ans;
	}
	il auto max(auto a,auto b){return (a>b?a:b);}
	il auto min(auto a,auto b){return (a<b?a:b);}
	il int gcd(int a,int b){
		if(!b) return a;
		return gcd(b,a%b);
	}
	il int lcm(int a,int b){
		return a/gcd(a,b)*b;
	}
	il void exgcd(int a,int b,int &x,int &y){
		if(!b) return x=1,y=0,void(0);
		exgcd(b,a%b,x,y);
		int t=x;
		x=y,y=t-a/b*x;
		return ;
	}
	mt19937 rnd(time(0));
}
using namespace yzqwq;

const int N=21,M=2e5+10;
int n,m;
string s_[N],s[N];
int f[1<<N][N];
int b[M],c[N][N];
int hs_,hs[N][M],bk=12347;
bool vis[N];

il int get(int now,int l,int r){
	return hs[now][r]-hs[now][l-1]*b[r-l+1];
}
il bool check(int now,string a,int hs_){
	if(s_[now].size()<=a.size()) return 0;
	for(re int l=0;l+a.size()-1<s_[now].size();++l){
		int r=l+a.size()-1;
		if(get(now,l+1,r+1)==hs_) return 1;
	}
	return 0;
}

il void solve(){
	n=rd,b[0]=1;
	for(re int i=1;i<M;++i) b[i]=b[i-1]*bk;
	for(re int i=1;i<=n;++i){
		cin>>s_[i];
		for(re int j=0;j<s_[i].size();++j)
			hs[i][j+1]=hs[i][j]*bk+s_[i][j]-'A'+1;
	}
	for(re int i=1;i<=n;++i){
		hs_=0;
		for(re int j=0;j<s_[i].size();++j)
			hs_=hs_*bk+s_[i][j]-'A'+1;
		for(re int j=1;j<=n;++j)
			if(check(j,s_[i],hs_)){
				vis[i]=1;break;
			}
	}
	for(re int i=1;i<=n;++i)
		if(!vis[i]){
			s[++m]=s_[i];
			for(re int j=0;j<s[m].size();++j)
				hs[m][j+1]=hs[m][j]*bk+s[m][j]-'A'+1;	
		}
	for(re int i=1;i<=m;++i)
	for(re int j=1;j<=m;++j){
		int Mx=0;
		for(re int len=1;len<=min(s[i].size(),s[j].size());++len)
			if(hs[j][len]==get(i,s[i].size()-len+1,s[i].size())) Mx=len;
		c[i][j]=Mx;		
	}
	for(re int i=0;i<(1<<m);++i)
	for(re int j=1;j<=m;++j)
		f[i][j]=1e9;
	for(re int i=1;i<=m;++i) f[1<<(i-1)][i]=s[i].size();
	for(re int i=0;i<(1<<m);++i)
	for(re int lst=1;lst<=m;++lst)
		if((i>>(lst-1))&1)
			for(re int now=1;now<=m;++now){
				if((i>>(now-1))&1) continue;
				f[i+(1<<(now-1))][now]=min(f[i+(1<<now-1)][now],f[i][lst]+s[now].size()-c[lst][now]);
			}
	int ans=1e9;
	for(re int i=1;i<=m;++i) ans=min(ans,f[(1<<m)-1][i]);
	printf("%lld\n",ans);
	return ;
}

signed main(){
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
	int t=1;while(t--)
	solve();
	return 0;
}

注:P2322 是这题的输出方案版。

posted @   harmis_yz  阅读(11)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 开源Multi-agent AI智能体框架aevatar.ai,欢迎大家贡献代码
· Manus重磅发布:全球首款通用AI代理技术深度解析与实战指南
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· AI技术革命,工作效率10个最佳AI工具
点击右上角即可分享
微信分享提示