HDU 4899

在看了一位大佬的代码过后总算慢慢写出来了。
发现这道题是我对LCS并不熟悉，对他的状态转移方程不够了解，看来要回去看看这个LCS这类题目了。
题目还是状压的思想，问题是status是什么？
我们这里用ai来表示 输入字符串s[1…i] 和用来比对的字符串t[1…j]时候的lcs值。
我们可以知道的是an-a(n-1)一定==0||==1的。通过LCS转移方程可以看出
if s[i]==t[j]
dp[i][j]=dp[i-1][j-1]+1;
else
dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
所以我们以差分的值作为状态，可以看出来，从0-len（s）,一共有len个01串，即为状态。
所以a[i]=a[i-1]+s>>（i-1)&1。枚举status的值，然后还原status。
然后这里还要思考，新状态，即t[1…j+1]添加了一个新的字母后，状态成了什么样子。
这里就成了新的状态转移了，即bi。
所以bi就用lcs来表示出来，表示出来后再重新算回去新的status并保存就ok。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<climits>
#define up(i,a,b)  for(int i=a;i<b;i++)
#define dw(i,a,b)  for(int i=a;i>b;i--)
#define upd(i,a,b) for(int i=a;i<=b;i++)
#define dwd(i,a,b) for(int i=a;i>=b;i--)
#define modulo 1e9+7
//#define local
typedef long long ll;
const double esp = 1e-6;
const double pi = acos(-1.0);
const long long INF = 0x3f3f3f3f;
const int maxn = 1e9 + 7;
using namespace std;
int t;
char str[100];
int m;
int dp[2][1<<18];
int len;
int trans[5][1 << 18];
int ans[20];
char exmp[] = "AGCT";
void add(int &x, int y)//这里开始模1e9+7
{
	x = (x + y) % maxn;
}
int countone(int x)//这里算的是有多少个1
{
	int sum=0;
	while (x)
	{
		if (x & 1)
			sum++;
		x >>= 1;
	}
	return sum;
}
void init()//初始化预处理转移的状态
{
	up(s, 0, (1 << len))
	{
		int a[20] = { 0 };
		upd(i, 1, len)
		{
			a[i] = a[i - 1] + ((s >> (i - 1)) & 1);
		}
		up(j, 0, 4)
		{
			int b[20] = { 0 };
			upd(i,1, len)
			{
				if (str[i] == exmp[j])
				{
					b[i] = a[i - 1] + 1;
				}
				else
					b[i] = max(b[i - 1], a[i]);
			}
			int new1=0;
			upd(i, 1, len)
				new1 |= b[i] - b[i - 1] << (i - 1);
			trans[j][s] = new1;
		}
	}
}
int main()
{
	cin >> t;
	while (t--)
	{
		memset(ans, 0, sizeof(ans));	
		cin >> str + 1;
		cin >> m;
		len = strlen(str + 1);
		init();
		int l = 0, r = 1;//这里学习大佬的做法用滚动数组
		memset(dp, 0, sizeof(dp));
		dp[0][0] = 1;
		up(i, 0, m)
		{
			up(j, 0, 4)
			{
				up(s, 0, (1 << len))
				{
					add(dp[r][trans[j][s]] , dp[l][s]);//新状态等于就状态加上
				}
			}
			swap(l, r);
			memset(dp[r], 0, sizeof(dp[r]));
		}
		up(s, 0, (1 << len))
		{
			add(ans[countone(s)], dp[l][s]);//这里算ans
		}
		upd(i, 0, len)
		{
			cout << ans[i];
                         cout<<endl;
		}
	}
	return 0;
}