NKOI 3747 宠物

Description
xxxxxyt学姐经常一个人在家，难免会感到寂寞，于是学姐养了n只可爱的宠物，比
如皮皮虾、大蟒蛇、藏狐、安康鱼…但即便如此学姐还是感到无聊。突然有一天，学姐想到
了让宠物们互相对战的消遣方法（请不要给动物保护协会打电话！）。学姐让宠物们两两进
行对战，n*(n-1)/2场对战后，学姐得到了一张相生相克图，然后又根据自己的喜好，把
宠物们分成了一队与二队。就在队伍分好后，学姐的强迫症又犯了，她希望自己的两支队伍
都满足这样一个性质：存在某种排列，使得排在后面的宠物能够击败排在前面的所有宠物。
但学姐的懒惰大家都是知道的，所以她找到了你，希望你能告诉她这两支队伍是否均满足要
求，如果是，她还希望你告诉她最多可以从二队中抽出多少只宠物放在一队，使得两支队伍
仍然满足要求。努力解决问题吧，而xxxxxyt学姐，瘫躺。
Input
第一行输入两个数字n和m，分别表示学姐有n只宠物，其中被分到一队的宠物有m
只。
接下来n行每行n个数字，ai,j
表示第i只宠物是否能战胜第j只宠物，保证ai,i=0
且ai,j=!aj,i
。
接下来一行m个数字，表示有哪些宠物被分到了一队。
Output
如果两支队伍均不能让xxxxxyt满意，则输出“NO”；否则输出“YES”，并输出一个
最大的k，使得从二队中非任意地抽出k只宠物放入一队后，两支队伍仍然满足条件。详细
格式见样例输出。
Sample Input 1
3 2
0 1 1
0 0 1
0 0 0
3 1
Sample Output 1
YES 1
Sample Input 2
4 3
0 1 0 1
0 0 1 1
1 0 0 1
0 0 0 0
1 2 3
Sample Output 2
NO
Sample Input 3
4 2
0 1 0 1
0 0 1 1
1 0 0 1
0 0 0 0
1 2
Sample Output 3
YES 1
Hint
注意：
宠物们的实力是相对的，也就是可能会出现A战胜B，B战胜C，C又战胜A的
情况。
数据范围：
20%的数据1<=m<n<=10
60%的数据1<=m<n<=100

100%的数据1<=m<n<=100

显然一个序列的规则是满足拓扑排序的，故可以用拓扑排序check，且合法的排列方法唯一，即拓扑序。现在问题就变成了有一队二队两条有向链，问二队链中至多有多少个可以合法地插入一队链，用f[i][j]表示一队链讨论到第i个二队
链讨论到第j个，二队链中至多插入f[i][j]个，转移为:

f[i][j]=max{f[i-1][j],f[i][j-1],f[i][j-1]+1(二队中第j个宠物弱于一队中1...i强于i+1...m)}
直接转移暴力判断复杂度为O(n^3)
因此对二队中第j个宠物是否弱于一队中1...i强于i+1...m进行预处理，便可在复杂度O(n2)下得到答案

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
const int maxn=1005,inf=1e9;
inline void _read(int &x){
    char t=getchar();bool sign=true;
    while(t<'0'||t>'9')
    {if(t=='-')sign=false;t=getchar();}
    for(x=0;t>='0'&&t<='9';t=getchar())x=x*10+t-'0';
    if(!sign)x=-x;
}
int n,m,type[maxn],du[maxn],f[maxn][maxn];
int s1[maxn],s2[maxn],l[maxn][maxn],r[maxn][maxn];
int map[maxn][maxn];
bool topsort(int s[],int kind,int sum){
	queue<int>q;
	int i,j,t=0;
	bool mark[maxn];
	memset(mark,0,sizeof(mark));
	for(i=1;i<=n+m;i++)
	    if(type[i]==kind&&(!du[i])){
	    	q.push(i);
	    	s[++t]=i;
	    	mark[i]=1;
		}
	while(q.size()){
		int x=q.front();q.pop();
		for(i=1;i<=n+m;i++)
		    if(!mark[i]&&type[i]==kind){
		    	if(map[x][i]&&du[i])du[i]--;
		    	if(!du[i]){
		    		q.push(i);
		    		mark[i]=1;
		    		s[++t]=i;
				}
			}
	}
	return t==sum;
}
int main(){
	int i,j,x;
	_read(n);_read(m);
	for(i=1;i<=n;i++){
		for(j=1;j<=n;j++)_read(map[i][j]);
		type[i]=2;
	}
	for(i=1;i<=m;i++){
		_read(x);
		type[x]=1;
	}
	for(i=1;i<=n;i++)
	    for(j=1;j<=n;j++)
	        if(type[i]==type[j]&&map[i][j])du[j]++;
	n-=m;swap(n,m);   
	if(!topsort(s1,1,n)||(!topsort(s2,2,m))){
		puts("NO");
		return 0;
	}
	for(j=1;j<=m;j++){
		l[0][j]=r[n+1][j]=1;
		for(i=1;i<=n;i++)l[i][j]= map[s1[i]][s2[j]]&&l[i-1][j];
		for(i=n;i;i--)r[i][j]= map[s2[j]][s1[i]]&&r[i+1][j];
	}
	for(i=0;i<=n;i++)
	    for(j=0;j<=m;j++){
	    	if(j)f[i][j]=max(f[i][j],f[i][j-1]);
	    	if(i)f[i][j]=max(f[i][j],f[i-1][j]);
	    	if(j&&l[i][j]&&r[i+1][j])f[i][j]=max(f[i][j],f[i][j-1]+1);
		}
	cout<<"YES"<<" "<<f[n][m];
}