UVa1625

/*************************************************************************
	> File Name: UVa1625.cpp
	> Author:HaoWei 
	> Mail:liang199611@outlook.com 
	> Created Time: 2016年06月21日 星期二 12时17分26秒
 ************************************************************************/

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<string>
#include<utility>

using namespace std;
const int maxn=5555;
#define me(a,b) memset(a,b,sizeof(a))
#define rep(i,n) for(int i=0;i<n;i++)
#define repl(i,n) for(int i=1;i<=n;i++)
int AS[27],AE[27],BS[27],BE[27];
int dp[2][maxn],c[2][maxn];

int main()
{
    string str1,str2;
    int T;cin>>T;
    while(T--)
    {
        cin>>str1>>str2;
        me(AS,0x3f);me(AE,-1);me(BS,0x3f);me(BE,-1);
        int n=str1.length(),m=str2.length();
        str1="#"+str1;str2="#"+str2;
        repl(i,n){
            int id=str1[i]-'A';
            if(AS[id]>i) AS[id]=i;
            if(AE[id]<i) AE[id]=i;
        }
        repl(i,m){
            int id=str2[i]-'A';
            if(BS[id]>i) BS[id]=i;
            if(BE[id]<i) BE[id]=i;
        }
        me(dp,0x3f);
        dp[0][0]=0,c[0][0]=0;
       // rep(i,26) printf("AS:%11d AE:%11d BS:%11d BE:%11d\n",AS[i],AE[i],BS[i],BE[i]);
        bool flag=false;
        for(int tt=0;tt<=n;tt++)
            for(int j=0;j<=m;j++)
            {
                int i=tt&1;
                if(tt>0)
                {
                    c[i][j]=c[i^1][j];
                    int id=str1[tt]-'A';
                    if(AS[id]==tt&&j<BS[id]) {c[i][j]++;}
                    if(AE[id]==tt&&j>=BE[id]) {c[i][j]--;}
                    dp[i][j]=dp[i^1][j]+c[i][j];
                    
                }
                if(j>0)
                {
                    c[i][j]=c[i][j-1];
                    int id=str2[j]-'A';
                    if(BS[id]==j&&tt<AS[id]) {c[i][j]++;}
                    if(BE[id]==j&&tt>=AE[id]) {c[i][j]--;}
                    dp[i][j]=min(dp[i][j-1]+c[i][j],dp[i][j]);
                }
           //   if(flag) cout<<"c ij"<<c[i][j]<<endl;
            }

        if(flag)
        rep(i,1)
        {
            cout<<endl;
            rep(j,m+1) printf("(%2d,%2d) ",dp[i][j],c[i][j]);
        }
        cout<<dp[n&1][m]<<endl;

    }
    return 0;
}

http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=51186

c[i][j]记录str1取i个字母str2取前j个字母时已经开始但未结束的字母的种类

dp(i,j)=min(dp[i-1][j]+c[i-1][j],dp[i][j-1]+c[i][j-1])

c[i][j]也是通过dp求得,T目测很多,需要用滚动数组以防memset时T掉