第六届福建省赛 - G - Simple String Problem （状压DP）

最新推荐文章于 2018-10-15 17:38:54 发布

原创最新推荐文章于 2018-10-15 17:38:54 发布 · 392 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

其他DP

状态压缩

本文介绍了一种使用状态压缩动态规划（状压DP）的方法来解决特定字符串问题的技术。该问题要求从由前k个字母组成的长度为n的字符串中找出两个不包含相同字母的子串，使这两个子串长度的乘积最大化。文中提供了两种实现思路：一种是通过广度优先搜索（BFS）进行记忆化搜索；另一种是采用普通的动态规划方法。

题目

题意：

给你n长的字符串，是由前k个字母组成的。

问你挑出两个子串，他们不能拥有相同的字母。

问你len(串a)*len(串b)最大是多少。

POINT：

状压。

考虑串长度只有2000。可以用2000*2000的效率算出每个状态最大的长度是多少。

然后如果枚举每个状态再和其他状态相比较得出答案，效率是（1<<k）*（1<<k）。显然是不行的

所以预处理出dp[i]，代表i这个状态和他的所有子集的最大长度是多少。

比如1111这个状态，包含了1111,0111,1011,1101,1110,0011,0101…… 等等到0001所有状态的最大长度。

然后只要枚举每个状态，找出反状态的dp就OK了。

可以普通dp，

记忆化bfs搜索, 我每次都把dp写成记忆化搜索

bfs：

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <queue>
using namespace std;

char s[2222];

int dp[70000],Max[70000];

int n,k;

void bfs()
{
	queue<int>q;
	int vis[70000];
	q.push(0);
	for(int i=0;i<(1<<k);i++)
		dp[i]=Max[i],vis[i]=0;
	while(!q.empty()){
		int now = q.front();
		q.pop();
		for(int i=0;i<k;i++){
			if((now&(1<<i))==0){
				int a=now|(1<<i);
				dp[a]=max(dp[now],dp[a]);
				if(!vis[a])
					q.push(a),vis[a]=1;
			}
		}
	}


}


int main()
{
	int T;
	scanf("%d",&T);
	while(T--){
		scanf("%d%d",&n,&k);
		scanf("%s",s+1);
		memset(Max,0,sizeof Max);
		for(int i=1;i<=n;i++){
			int now = 0;
			for(int j=i;j<=n;j++){
				now=now|(1<<(s[j]-'a'));
				Max[now]=max(Max[now],j-i+1);
			}
		}
		bfs();
		int ans = 0;
		int b=(1<<k)-1;
		for(int i=1;i<(1<<k);i++){
			ans=max(ans,dp[i]*dp[i^b]);
		}
		printf("%d\n",ans);
	}
}

普通dp：

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <queue>
using namespace std;

char s[2222];

int dp[70000],Max[70000];

int n,k;

int main()
{
	int T;
	scanf("%d",&T);
	while(T--){
		scanf("%d%d",&n,&k);
		scanf("%s",s+1);
		memset(Max,0,sizeof Max);
		for(int i=1;i<=n;i++){
			int now = 0;
			for(int j=i;j<=n;j++){
				now=now|(1<<(s[j]-'a'));
				Max[now]=max(Max[now],j-i+1);
			}
		}
		for(int i=0;i<(1<<k);i++)
			dp[i]=Max[i];
		int b=(1<<k)-1;
		for(int 	i=1;i<=b;i++){
			for(int j=0;j<k;j++){
				if(i&(1<<j)){
					dp[i]=max(dp[i],dp[i-(1<<j)]);
				}
			}
		}
		int ans = 0;

		for(int i=1;i<(1<<k);i++){
			ans=max(ans,dp[i]*dp[i^b]);
		}
		printf("%d\n",ans);
	}
}