Color Length UVA - 1625 (线性dp)

题目连接:点击打开链接

题目思路:
这道题是紫书上的一道题,其实很明显是一道dp题,为什么呢,因为对于每一个字符什么时候去加入那个队列,都会产生一颗新的解答树,不用dp的话,大概复杂度,是o(2^n),其实我们仔细去发现对于父解答树,他只关心下一个子解答树的最优解,并不关心怎么来的,那么我们就可以想到用dp了,这道题定义状态还是非常简单的,dp(i,j)表示第一个移走了i个,第二个移走了j个,的最小花费,但是问题来了,怎么去计算当前花费呢,根据题目意思去直接计算?这样的话我们并不能计算当前状态的最小花费,而且不知道颜色什么时候结束,什么时候开始,这样子的话就不能去状态转移了,那么我们该怎么去定义呢,其实我们最大的难点就是不知道某个字符串什么时候结束,什么时候开始,从而不知道怎么去计算此状态的花费,那么我们尝试不去考虑这个状态什么时候开始,什么时候结束,我们对于某个字符加入新的队列的时候我们让当前已经开始但没有结束的字符的花费都加一,这样的话,我们就不用考虑什么时候开始什么呢时候结束,比如对于abb 和bab的当前花费是一样的,这样的话我们就只用按照这个计算方式去从dp(i-1,j)和dp(i,j-1)转移状态了,还剩下一个问题,我们怎么去记录当前状态开始但未结束的字符数量呢,这里也是这道题的难点,我们可以这样,先记录每一种字符在每一个数列开始和结束的位置,然后之dp的状态转移的时候开一个cnt【】【】去记录当前状态已经开始但未结束的字符数量,考率当前状态,如果这个字符是第一次出现那么cnt就++,如果不是第一次出现但是也没有结束,就不用管,因为只用考虑种类,如果结束了那么cnt--;

ac代码:

#include<cstdio>
#include<iostream>
#include<cstdio>
#include<sstream>
#include<algorithm>
#include<cstring>
#include<queue>
#include<set>
#include<map>
#define LL long long
#define INF 0x3f3f3f3f
using namespace std;
const int maxn = 5000+100;
int T;
char s1[maxn],s2[maxn];
int dp[maxn][maxn];
int cnt[maxn][maxn];
int sc1[27];
int sc2[27];
int ec1[27];
int ec2[27];
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%s",s1+1);
        scanf("%s",s2+1);
        int len1 = strlen(s1+1);
        int len2 = strlen(s2+1);
        for(int i = 0 ;i<=len1;i++)
        {
            for(int j = 0 ;j<=len2;j++)
            {
                dp[i][j] = INF;
                cnt[i][j] = 0;
            }
        }
        memset(sc1,INF,sizeof(sc1));
        memset(sc2,INF,sizeof(sc2));
        memset(ec2,0,sizeof(ec2));
        memset(ec1,0,sizeof(ec2));
        for(int i = 1 ;i<=len1;i++)
        {
            sc1[s1[i]-'A'] = min(i,sc1[s1[i]-'A']);
            ec1[s1[i]-'A'] =i;
        }
        for(int j = 1 ;j<=len2;j++)
        {
             sc2[s2[j]-'A'] = min(j,sc2[s2[j]-'A']);
            ec2[s2[j]-'A'] = j;
        }
        dp[0][0] = 0;
        int num = 0;
        int sum = -1;
        for(int i = 0;i<=len1;i++)
        {
            for(int j  = 0;j<=len2;j++)
            {
                if(i==0||j==0){
                if(i==0&&j==0)
                    continue;
                if(i==0&&j!=0)
                {
                    dp[i][j] = dp[i][j-1]+cnt[i][j-1];
                    cnt[i][j] = cnt[i][j-1];
                   if(sc2[s2[j]-'A']==j&&sc1[s2[j]-'A']>i)
                    cnt[i][j]++;
                   if(ec2[s2[j]-'A']==j&&ec1[s2[j]-'A']<=i)
                    cnt[i][j]--;
                }
                if(i!=0&&j==0)
                {
                    dp[i][j] = dp[i-1][j]+cnt[i-1][j];
                       cnt[i][j] = cnt[i-1][j];
                   if(sc1[s1[i]-'A']==i&&sc2[s1[i]-'A']>j)
                    cnt[i][j]++;
                   if(ec1[s1[i]-'A']==i&&ec2[s1[i]-'A']<=j)
                    cnt[i][j]--;
                }
                }
                else{
                if(dp[i-1][j]+cnt[i-1][j]<dp[i][j-1]+cnt[i][j-1]){
                        dp[i][j] =dp[i-1][j]+cnt[i-1][j];
                      cnt[i][j] = cnt[i-1][j];
                   if(sc1[s1[i]-'A']==i&&sc2[s1[i]-'A']>j)
                    cnt[i][j]++;
                   if(ec1[s1[i]-'A']==i&&ec2[s1[i]-'A']<=j)
                    cnt[i][j]--;
                }
                else
                {
                     dp[i][j] = dp[i][j-1]+cnt[i][j-1];
                    cnt[i][j] = cnt[i][j-1];
                   if(sc2[s2[j]-'A']==j&&sc1[s2[j]-'A']>i)
                    cnt[i][j]++;
                   if(ec2[s2[j]-'A']==j&&ec1[s2[j]-'A']<=i)
                    cnt[i][j]--;
                }
                }
              // cout<<"i:"<<i<<' '<<"j:"<<j<<' '<<dp[i][j]<<endl;
            }
        }
        cout<<dp[len1][len2]<<endl;
    }
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值