bzoj1194

O(S²)枚举2个诅咒机, 然后O(n²)BFS去判断. 构成一个有向图, tarjan缩点, 然后就是求DAG的最长路——题解来自lsj

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#define clr(a,x) memset(a,x,sizeof(a))
#define rep(i,l,r) for(int i=(l);i<(r);i++)
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
#define mkp(a,b) make_pair(a,b)
int read(){
    int ans=0,f=1;
    char c=getchar();
    while(!isdigit(c)){
        if(c=='-') f=-1;
        c=getchar();
    }
    while(isdigit(c)){
        ans=ans*10+c-'0';
        c=getchar();
    }
    return ans*f;
}
const int maxs=59,maxn=59,maxm=59;
struct edge{
    int v;
    edge*next;
}e[maxn*maxn],*fir[maxn],*pt=e,E[maxn*maxn],*Fir[maxn],*Pt=E;
void addedge(int u,int v){
    pt->v=v;pt->next=fir[u];
    fir[u]=pt++;
}
void Addedge(int u,int v){
    Pt->v=v;Pt->next=Fir[u];
    Fir[u]=Pt++;
}
int s,n[maxs],m[maxs],v[maxs][maxn][2];
bool vis[maxn][maxn],f[maxs][maxn];
bool bfs(int x,int y){
    queue<pii>Q;Q.push(mkp(0,0));
    clr(vis,0);
    while(!Q.empty()){
        int a=Q.front().first,b=Q.front().second;Q.pop();
        rep(i,0,2){
            int U=v[x][a][i],V=v[y][b][i];
            if(vis[U][V]) continue;
            if(f[x][U]&&!f[y][V]) return 0;
            Q.push(mkp(U,V));vis[U][V]=1;
        }
    }
    return 1;
}
int ans,dfstime,scccnt,size[maxn],pre[maxn],low[maxn],scc[maxn];
stack<int>S;
void dfs(int x){
    pre[x]=low[x]=++dfstime;S.push(x);
    for(edge*e=fir[x];e;e=e->next){
        if(!pre[e->v]){
            dfs(e->v);
            low[x]=min(low[x],low[e->v]);
        }else if(!scc[e->v]) low[x]=min(low[x],pre[e->v]);
    }
    if(pre[x]==low[x]){
        int t=-1;++scccnt;
        while(t!=x){
            t=S.top();S.pop();
            scc[t]=scccnt;size[scccnt]++;
        }
    }
}
void tarjan(){
    rep(i,0,s) if(!scc[i]) dfs(i);
    clr(vis,0);
    rep(i,0,s){
        for(edge*e=fir[i];e;e=e->next){
            if(scc[i]!=scc[e->v]&&!vis[scc[i]][scc[e->v]]){
                vis[scc[i]][scc[e->v]]=1;Addedge(scc[i],scc[e->v]);
            }
        }
    }
}
bool in[maxn];
int w[maxn];
void dp(int x){
    if(w[x]) return;
    for(edge*e=Fir[x];e;e=e->next){
        dp(e->v);w[x]=max(w[x],w[e->v]);
    }
    w[x]+=size[x];
}
void solve(){
    tarjan();
    rep(i,1,scccnt+1){
        for(edge*e=Fir[i];e;e=e->next) in[e->v]=1;
    }
    rep(i,1,scccnt+1) if(!in[i]) dp(i);
    rep(i,1,scccnt+1) ans=max(ans,w[i]);
    printf("%d\n",ans);
}
int main(){
    s=read();
    rep(i,0,s){
        n[i]=read();m[i]=read();
        rep(j,1,m[i]+1){
            int t=read();
            f[i][t]=1;
        }
        rep(j,0,n[i]){
            v[i][j][0]=read();v[i][j][1]=read();
        }
    }
    rep(i,0,s){
        rep(j,0,s) if(i!=j&&bfs(i,j)) addedge(i,j);
    }
    solve();
    return 0;
}

1194: [HNOI2006]潘多拉的盒子

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 255  Solved: 133
[Submit][Status][Discuss]

Description

 

Input

第一行是一个正整数S，表示宝盒上咒语机的个数，（1≤S≤50）。文件以下分为S块，每一块描述一个咒语机，按照咒语机0,咒语机1„„咒语机S－1的顺序描述。每一块的格式如下。 一块的第一行有两个正整数n,m。分别表示该咒语机中元件的个数、咒语源输出元的个数（1≤m≤n≤50）。 接下来一行有m个数，表示m个咒语源输出元的标号（都在0到n-1之间）。接下来有n行，每一行两个数。第i行（0≤i≤n-1）的两个数表示pi,0和pi,1（当然，都在0到n-1之间）。

Output

第一行有一个正整数t,表示最长升级序列的长度。

Sample Input

4
1 1
0
0 0
2 1
0
1 1
0 0
3 1
0
1 1
2 2
0 0
4 1
0
1 1
2 2
3 3
0 0

Sample Output

3

HINT

 

Source

 

[Submit][Status][Discuss]