动态规划(牛客上的几个题)

1.被3整除的

01背包问题

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
string s;
int dp[57];
int main()
{
    cin>>s;
    int len=s.length(),ans=0,sum=0;
    for(int i=0;i<len;i++)
    {
        sum+=s[i]-'0';
    }
    dp[0]=1;
    for(int i=0;i<len;i++)
    for(int j=sum;j>=s[i]-'0';j--)
    {
        dp[j]=(dp[j-s[i]+'0']%mod+dp[j]%mod)%mod;
    }
    dp[0]-=1;
    for(int i=0;i<=sum;i++)
    {
        if(i%3==0&&dp[i])
        ans=(ans%mod+dp[i]%mod)%mod;
     }
    printf("%d\n",ans);
    return 0;
}

2.删括号

3维dp，dp{i][j][k]表示当前匹配到s1第i个位置与s2第j个位置删除的‘(’与‘)’的差值为k
只有当k等于0时我们才可能能继续往下更新j，因为删括号是成对删的

#include<bits/stdc++.h>
using namespace std;
int dp[107][107][57];
string s1,s2;
int len1,len2;
int main()
{
    cin>>s1>>s2;
    len1=s1.length();
    len2=s2.length();
    dp[0][0][0]=1;
    for(int i=0;i<len1;i++)
    for(int j=0;j<len2;j++)
    for(int k=0;k<len1/2;k++)
    {
       if(dp[i][j][k])
       {
           if(k==0&&s1[i+1]==s2[j+1]) dp[i+1][j+1][k]=1;
           if(s1[i+1]=='(') dp[i+1][j][k+1]=1; //之所以不是else if是因为如果k为0且下一个能匹配上，不是一定要保留的，因为下下一个可能就匹配不上了，所以我们也可能将这个位置的字符删除
           else if(s1[i+1]==')') dp[i+1][j][k-1]=1;
       }
    }
    if(dp[len1-1][len2-1][0]) printf("Possible\n");
    else printf("Impossible\n");
    return 0;
}

3.美丽序列

这个题就是纯暴力，这是我做过的最暴力的一个动态规划题，思维dp[i][j][k][z]表示当前处理到第n个数，该位置之前以及该位置的连续的递减数目为j，k是当前的处理状态的最后一个数的值，z是这i个数的前缀和，然后暴力枚举，一重循环枚举i，遇到-1，三重循环分别枚举-1要替换的数、-1前面的数以及符合条件的前缀和

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
long long n,dp[47][3][47][1607];
long long a[47];
int main()
{
    scanf("%lld",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&a[i]);
    }
    if(a[1]==-1) for(int i=0;i<=40;i++) dp[1][1][i][i]=1;
    else dp[1][1][a[1]][a[1]]=1;
    for(int i=2;i<=n;i++)
    {
        if(a[i]==-1)
        {
             for(int j=0;j<=40;j++)
             for(int k=0;k<=40;k++)
             for(int z=(i-1)*j;z<=(i-1)*40;z++){
             if(j>=k){
             dp[i][1][j][z+j]=(dp[i][1][j][z+j]+dp[i-1][1][k][z])%mod;
             dp[i][1][j][z+j]=(dp[i][1][j][z+j]+dp[i-1][2][k][z])%mod;
             }
             else
             dp[i][2][j][z+j]=(dp[i][2][j][z+j]+dp[i-1][1][k][z])%mod;
         }
        }
         else
         {
             for(int k=0;k<=40;k++)
             for(int z=(i-1)*a[i];z<=(i-1)*40;z++)
             {
             if(k<=a[i])
             {
             dp[i][1][a[i]][z+a[i]]=(dp[i][1][a[i]][z+a[i]]+dp[i-1][1][k][z])%mod;
             dp[i][1][a[i]][z+a[i]]=(dp[i][1][a[i]][z+a[i]]+dp[i-1][2][k][z])%mod;
             }
             else
             dp[i][2][a[i]][z+a[i]]=(dp[i][2][a[i]][z+a[i]]+dp[i-1][1][k][z])%mod;
             }
         }
    }
    long long ans=0;
    for(int i=0;i<=40;i++)
    for(int j=i*n;j<=40*n;j++)
    {
        ans=(ans+dp[n][1][i][j])%mod;
        ans=(ans+dp[n][2][i][j])%mod;
    }
    cout<<ans<<endl;
    return 0;
}

4. codeforce
   
   这个题需要首先判断这n个题目要做的先后顺序
   对于两道题t1,t2来说，有两种做题顺序，C12,C21,我们定义P1为t1的每分钟减小的分数，T1为做题需要的时间对于题1与2有两种顺序
   1、2得分为C12=sum-T1P1-(T1+T2)P2;
   2、1得分为C12=sum-T2P2-(T1+T2)P1;
   则 C12-C21=T2P1-T1P2
   要判断两道题的大小，这上式同时除T1*T2的P1/T1-P2/T2 这样我们就发现了优先级函数为P/T即每分钟减小分数除需要的做题时间
   排好序之后，进行01背包的dp算法

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=100007;
int t,n;
ll dp[maxn];
struct node
{
    int maxpoints;
    int pointsperminute;
    int requiredtime;
    bool operator < (const node &w)const
    {
        return requiredtime*1.0/pointsperminute<w.requiredtime*1.0/w.pointsperminute;
    }
}a[57];
int main()
{
    scanf("%d%d",&n,&t);
    for(int i=1;i<=n;i++)
    scanf("%d",&a[i].maxpoints);
    for(int i=1;i<=n;i++)
    scanf("%d",&a[i].pointsperminute);
    for(int i=1;i<=n;i++)
    scanf("%d",&a[i].requiredtime);
    sort(a+1,a+1+n);
    long long ans=0;
    for(int i=1;i<=n;i++)
    for(int j=t;j>=a[i].requiredtime;j--)
    {
        dp[j]=max(dp[j-a[i].requiredtime]+a[i].maxpoints-j*a[i].pointsperminute,dp[j]);
        ans=max(ans,dp[j]);
    }
    cout<<ans<<endl;
    return 0;
}