P428 交错匹配（dp-双进程）

最新推荐文章于 2019-09-06 10:33:12 发布

原创最新推荐文章于 2019-09-06 10:33:12 发布 · 449 阅读

2 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#dp

双进程。专栏收录该内容

4 篇文章

订阅专栏

博客内容介绍了交错匹配问题的解决思路，来源于某博客。该问题与轮船问题类似，但要求匹配必须且只能相交一次。通过动态规划（DP）的方法，定义f[i][j]表示up前i个数与down前j个数的最大匹配数。状态转移方程为f[i][j]=max(f[i-1][j-1],f[i-1][j],f[i][j-1])，并考虑当up[i]！=down[j]时的特殊情况，利用预处理找到最近的匹配项，避免O(n^4)的时间复杂度，转而实现O(n^2)的预处理方案。" 133453538,19671514,Android模拟器中模拟GPS位置编程指南,"['Android开发', '模拟器', 'GPS', '编程']

这里写图片描述

思路来源：http://blog.sina.com.cn/s/blog_51cea4040100fgo4.html

想到轮船问题了，那边是不允许相交，这里是必须且只能与一个匹配相交。
再仔细一想，两者什么关系都没有。

f[i][j]表示up取前i个数，down取前j个数时的最大匹配数，注意，不是i和j匹配。
仍然先写状态转移方程：
f[i][j]=max(f[i-1][j-1],f[i-1][j],f[i][j-1]);
因为k！=l，即两个交叉的匹配的数字不能相同，故当up[i]==down[j]时不执行匹配。

假设存在kup，满足在doan的前j个数中，down[kup]==up[i]且kup距离j最近。
假设存在kdown，满足在up的前i个数中，up[kdown]==down[i]且kdown距离i最近。

则，当up[i]！=down[j]时： f[i][j]=max(f[i][j],f[kup-1][kdown-1]+2);
即up[i]与down[kup]匹配，down[j]与up[kdown]匹配。

如何求kup以及kdown？最暴力的做法是在状态转移的过程当中每次搜一遍，O(n^4)，很危险。所以我们要预处理。
在这里采用了O(n^2)的预处理，看了几个同学的代码，思路好像还挺多的。

void first_()
{
    memset(Up,0,sizeof(Up));
    memset(Down,0,sizeof(Down));
    for(int i=1;i<=n;++i)//down的前j个元素中，和up[i]相等的最靠近j下标的kup 。 
        for(int j=2;j<=m;++j)
            if(down[j-1]==up[i]) Up[i][j]=j-1;//为什么与j-1比较而不与j比较？因为Up[i][j]==j是不合法的，我们不允许i与j匹配。 
            else Up[i][j]=Up[i][j-1];
    for(int i=1;i<=m;++i)//up的前j个元素中，与 down[i]相等的最靠近j下标的kdown 。 
        for(int j=2;j<=n;++j)
            if(up[j-1]==down[i]) Down[j][i]=j-1;
            else Down[j][i]=Down[j-1][i];//这里为了dp时的方便（第一维仍然属于up数组），所以j和i的位置有点乱。
} 

void init()
{
    read(n);read(m);
    for(int i=1;i<=n;++i) read(up[i]);
    for(int i=1;i<=m;++i) read(down[i]);
}

void work()
{
    first_();
    memset(f,0,sizeof(f));
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j)
        {
            f[i][j]=max(max(f[i-1][j],f[i][j-1]),f[i-1][j-1]);
            if(Down[i][j]&&Up[i][j]&&up[i]!=down[j]) //要特判掉kdown==0||kup==0的情况，因为下标还要-1。 
                f[i][j]=max(f[i][j],f[Down[i][j]-1][Up[i][j]-1]+2);

        }
    printf("%d",f[n][m]);
}