Codeforces 463D Gargari and Permutations：隱式圖dp【多串LCS】

題目連結：http://codeforces.com/problemset/problem/463/D

題意：

　　給你k個1到n的排列，問你它們的LCS（最長公共子序列）是多長。

 

題解：

　　因為都是1到n的排列，即每個串中，1到n每個數字恰好出現一次。

　　將相同的數字之間相連，可以得到下面的樣子(n = 4, k = 3)：

　　

　　顯然，要求的LCS就等於圖中互不相交的最多連線個數。

 

　　將每一個數字看做一個節點。

　　若i到j有一條有向邊，則代表：

　　　　數字j的連線在i的連線的後面，且互不相交。

　　即：

　　　　若i->j，則要滿足所有的pos[k][i] <= pos[k][j]。

　　　　其中pos[k][i]表示第k個串中，數字i出現的位置。

 

　　O(N^2*K)建圖，最終得到的一定是一個有向無環圖。

　　LCS就等於這個圖上的最長路徑長度。

　　所以dfs跑一邊dp就行了。

 

AC Code:

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <vector>
#define MAX_N 1005
#define MAX_K 10

using namespace std;

int n,k;
int dp[MAX_N];
int a[MAX_K][MAX_N];
int pos[MAX_K][MAX_N];
vector<int> edge[MAX_N];

void read()
{
    cin>>n>>k;
    for(int i=1;i<=k;i++)
    {
        for(int j=1;j<=n;j++)
        {
            cin>>a[i][j];
            pos[i][a[i][j]]=j;
        }
    }
}

bool is_valid(int x,int y)
{
    for(int i=1;i<=k;i++)
    {
        if(pos[i][x]>=pos[i][y]) return false;
    }
    return true;
}

void build()
{
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(is_valid(i,j)) edge[i].push_back(j);
        }
    }
}

void dfs(int now)
{
    dp[now]=1;
    for(int i=0;i<edge[now].size();i++)
    {
        int temp=edge[now][i];
        if(dp[temp]==-1) dfs(temp);
        dp[now]=max(dp[now],dp[temp]+1);
    }
}

void work()
{
    build();
    memset(dp,-1,sizeof(dp));
    int ans=0;
    for(int i=1;i<=n;i++)
    {
        if(dp[i]==-1) dfs(i);
        ans=max(ans,dp[i]);
    }
    cout<<ans<<endl;
}

int main()
{
    read();
    work();
}