poj 1187 陨石的秘密

动态规划

一开始定义为dp[l1][l2][l3][dd],表示用l1个{},l2个[],l3个(),深度为dd的方案数,后来不行。参考别人后,dp[l1][l2][l3][dd]表示深度小于等于dd的方案数,那么答案是

dp[l1][l2][l3][dd]-dp[l1][l2][l3][dd-1]

可以用分块和整块的思想来看这个问题

好像[[()][]] 这类的我们可以说它是整块,[][](),这类的我们可以说它是分块,我们在构建的dp数组的时候就是把其当做是分块来处理,然后使用乘法原理

例如dp[l1][l2][l3][dd],我们可以先构建一个小分块,用小分块的方案数*其实括号能组成的方案数,而这个小分块,其实本身是一个整块

我们来看这个小分块,最外层可以是(),[],{}

如果是(),那么里面必定全部是(),由题意的优先级可知,所以假设用i个()组成小分块,则dp[l1][l2][l3][dd]=dp[0][0][i-1][dd-1]*dp[l1][l2][l3-i][dd]

也就是用i个()组成小分块,已经用掉一个做外层,所以只剩下i-1个,并且深度也应该下降为dd-1,而剩下l1个{},l2个[],l3-i个(),他们再组合,两者相乘

同样的,如果小分块使用[]作为外层,那么里面可以使用[]和(),整个小分块假设使用了j个[],i个(),则dp[l1][l2][l3][dd]=dp[0][j-1][i][dd-1]*dp[l1][l2-j][l3-i][dd]

意思就是,用j个[],i个()去构建小分块,已经用掉一个[]作为外层,因此深度也应该下降为dd-1,剩下l1个{},l2-j个[],l3-i个(),他们再组合,两者相乘

同样的,如果小分块使用{}作为外层,那么里面可以使用{},[],(),假设分别用了k,j,i个,那么dp[l1][l2][l3][dd]=dp[k-1][j][i][dd-1]*dp[l1-k][l2-j][l3-i][dd]

小分块用掉一个{}作为外层,深度也应该下降为dd-1,再用剩下的括号去组合,两者相乘

整个算法中需要按照优先级去枚举小分块最外层的类型,应该()在最外层,接着是[],{},这个也是 值得思考的地方

接着就是记忆化实现,感觉这东西递推不好写

搜索的话就要知道边界

1.l1=0 && l2=0 && l3=0  , dp值为1

1.dd=0时,dp值为0

#include <cstdio>
#include <cstring>
#define MOD 11380
//#define LOCAL


int dp[15][15][15][35];
bool vis[15][15][15][35];

int dfs(int l1,int l2,int l3,int dd)
{
    if(!l1 && !l2 && !l3)
    {
        vis[l1][l2][l3][dd] = true;
        return dp[l1][l2][l3][dd] = 1;
    }
    if(!dd)
    {
        vis[l1][l2][l3][dd] = true;
        return dp[l1][l2][l3][dd] = 0;
    }
    if(vis[l1][l2][l3][dd])
        return dp[l1][l2][l3][dd];

    int ans = 0;
    for(int i=0;i<=l3;i++)
    {
        if(i)
        {
            ans = (ans+dfs(0,0,i-1,dd-1) * dfs(l1,l2,l3-i,dd)) % MOD;
            ans %= MOD;
        }
        for(int j=0;j<=l2;j++)
        {
            if(j)
            {
                ans = (ans+dfs(0,j-1,i,dd-1) * dfs(l1,l2-j,l3-i,dd)) % MOD;
                ans %= MOD;
            }
            for(int k=1;k<=l1;k++)
            {
                ans = (ans+dfs(k-1,j,i,dd-1) * dfs(l1-k,l2-j,l3-i,dd)) % MOD;
                ans %= MOD;
            }
        }
    }
    vis[l1][l2][l3][dd] = true;
    return dp[l1][l2][l3][dd] = ans;
}



int main()
{

    #ifdef LOCAL
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);
    #endif

    int l1,l2,l3,dd;

    memset(vis,0,sizeof(vis));
    scanf("%d%d%d%d",&l1,&l2,&l3,&dd);

    dp[l1][l2][l3][dd]=dfs(l1,l2,l3,dd);
    if(dd) dp[l1][l2][l3][dd-1]=dfs(l1,l2,l3,dd-1);
    //dp[l1][l2][l3][dd-1]=dfs(l1,l2,l3,dd-1);

    if(!dd) { printf("%d\n",dp[l1][l2][l3][dd]); return 0; }
    int a1=dp[l1][l2][l3][dd] , a2=dp[l1][l2][l3][dd-1];

    //printf("%d %d\n",a1,a2);
    printf("%d\n",((a1-a2)%MOD+MOD)%MOD);
    return 0;
}
原文地址:https://www.cnblogs.com/scau20110726/p/2941322.html