HDU 5811 Colosseo

首先判断一下两个集合是否能够拓扑排序，顺便记录下每个节点的拓扑序。

然后看T2中每个点在T1中能够放在哪一个位置，记录下这个位置Pi。

然后T2中（按拓扑序排好），计算Pi的一个非严格递增的LIS。LIS长度就是答案。

这题scanf读入都900+ms了，有时直接卡TLE。改用gets整行读入，时间上会有很大改进。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<iostream>
using namespace std;
typedef long long LL;
const double pi=acos(-1.0),eps=1e-;
void File()
{
    freopen("D:\\in.txt","r",stdin);
    freopen("D:\\out.txt","w",stdout);
}
inline int read()
{
    char c = getchar(); while(!isdigit(c)) c = getchar(); int x = ;
    while(isdigit(c)) { x = x *  + c - ''; c = getchar();  }
    return x;
}

const int maxn=;
int n,m,g[maxn][maxn],r[maxn],d[maxn],a[maxn],LIS[maxn],dp[maxn];
bool f[maxn];
struct Edge {int u,v,nx;}e[maxn*maxn];
int sz,h[maxn];
char s[];

bool Top(int x)
{
    queue<int>Q; sz=; memset(h,-,sizeof h); memset(r,,sizeof r);
    for(int i=;i<=n;i++) for(int j=;j<=n;j++)
        if(f[i]==x&&f[j]==x&&g[i][j]==)
                r[j]++, e[sz].u=i, e[sz].v=j, e[sz].nx=h[i], h[i]=sz++;
    sz=;
    for(int i=;i<=n;i++) if(f[i]==x&&r[i]==) Q.push(i);
    while(!Q.empty())
    {
        int t=Q.front(); Q.pop(); sz++; d[t]=sz;
        for(int i=h[t];i!=-;i=e[i].nx)
            { r[e[i].v]--; if(r[e[i].v]==) Q.push(e[i].v); }
    }
    if(x==) { if(sz!=m) return ; return ; }
    else { if(sz!=n-m) return ; return ; }
}

int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        getchar(); if(n==&&m==) break;
        for(int i=;i<=n;i++)
        {
            gets(s); int p=;
            for(int j=;j<=n;j++) { g[i][j]=s[p]-''; p=p+; }
        }
        memset(f,,sizeof f); for(int i=;i<=m;i++) {int x; scanf("%d",&x),f[x]=;}
        bool x1=Top(), x2=Top();
        if(x1==||x2==) { printf("NO\n"); continue; }
        for(int i=;i<=n;i++)
        {
            if(f[i]) continue; int L=m, R=;
            for(int j=;j<=n;j++)
            {
                if(!f[j]) continue;
                if(g[i][j]) L=min(L,d[j]-); else R=max(R,d[j]);
            }
            if(L!=R) a[i]=-; else a[i]=L;
        }
        memset(LIS,-,sizeof LIS);
        for(int i=;i<=n;i++) { if(f[i]) continue; LIS[d[i]]=a[i]; }
        int ans=; memset(dp,,sizeof dp);
        for(int i=;i<=n-m;i++)
        {
            if(LIS[i]==-) continue;
            int pre=;  for(int j=;j<i;j++) if(LIS[j]<=LIS[i]) pre=max(pre,dp[j]);
            dp[i]=pre+; ans=max(ans,dp[i]);
        }
        printf("YES %d\n",ans);
    }
    return ;
}