1194: [HNOI2006]潘多拉的盒子 - BZOJ

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

终于AC了

网上说的那个二元组的确很神奇，用bfs判断出两个咒语机的关系，初始状态为（0，0），状态（x，y）表示在A咒语机上x格停留，在B咒语机上y格停留

当有一个状态（x，y）x是输出元，而y不是的话，则B咒语机有一个咒语不能产生而A能产生

现在我们求出了一个有向图，要求最长链，首先就是要消去环，我们发现咒语机不可能出现A-->B-->C-->A但是A不是C的升级，相当于A>=B>=C>=A那么A，B，C相等

所以我们对点进行染色，对于一个没有染色的点，我们把和它等价的点涂成一样的颜色（也就是一个强连通分量），只要i到j有边，j到i有边，i和j就是一个强连通分量里的，否则就不是

染完色之后缩点，就是一个森林了，于是dp出最大权值链（这时点上有权，就是涂成这个颜色的点的个数）

主要是因为如果打tarjan太长了，根本不想打，而且只有50个点，明显就是给我们乱搞的嘛

 var
     p:array[..,..,..]of longint;
     flag,l:array[..,..]of boolean;
     n:longint;

 function max(x,y:longint):longint;
 begin
     if x>y then exit(x);
     exit(y);
 end;

 procedure init;
 var
     i,j,a,b,k:longint;
 begin
     read(n);
     for i:= to n do
       begin
         read(a,b);
         for j:= to b do
           begin
             read(k);
             flag[i,k]:=true;
           end;
         for j:= to a- do
           read(p[i,j,],p[i,j,]);
       end;
 end;

 var
     d:array[..,..]of longint;
     v:array[..,..]of boolean;

 procedure work;
 var
     i,j,k,head,tail:longint;
     flagi,flagj:boolean;
 begin
     for i:= to n- do
       for j:=i+ to n do
         begin
           head:=;
           tail:=;
           d[,]:=;
           d[,]:=;
           fillchar(v,sizeof(v),true);
           v[,]:=false;
           flagi:=true;
           flagj:=true;
           while head<=tail do
             begin
               if flag[i,d[head,]]<>flag[j,d[head,]] then
               begin
                 if flag[i,d[head,]] then flagj:=false;
                 if flag[j,d[head,]] then flagi:=false;
               end;
               if (flagi=false)and(flagj=false) then break;
               for k:= to  do
                 if v[p[i,d[head,],k],p[j,d[head,],k]] then
                 begin
                   v[p[i,d[head,],k],p[j,d[head,],k]]:=false;
                   inc(tail);
                   d[tail,]:=p[i,d[head,],k];
                   d[tail,]:=p[j,d[head,],k];
                 end;
               inc(head);
             end;
           if flagi then l[i,j]:=true;
           if flagj then l[j,i]:=true;
         end;
 end;

 var
     map:array[..,..]of boolean;
     c,f,du,s:array[..]of longint;
     vis:array[..]of boolean;
     col,num,ans:longint;

 procedure get;
 var
     i,j:longint;
 begin
     for i:= to n do
       if c[i]= then
       begin
         inc(col);
         c[i]:=col;
         inc(s[col]);
         for j:= to n do
           if l[i,j] and l[j,i] then
           begin
             c[j]:=col;
             inc(s[col]);
           end;
       end;
     for i:= to n do
       for j:= to n do
         if c[i]<>c[j] then
         if l[i,j] then map[c[i],c[j]]:=true;
     for i:= to col do
       for j:= to col do
         if map[i,j] then inc(du[j]);
     fillchar(vis,sizeof(vis),true);
     num:=col;
     while num> do
       begin
         for i:= to col do
           if (vis[i])and(du[i]=) then
           begin
             vis[i]:=false;
             dec(num);
             for j:= to col do
               if map[i,j] then
               begin
                 dec(du[j]);
                 f[j]:=max(f[j],f[i]+s[i]);
               end;
           end;
       end;
     for i:= to col do
       ans:=max(ans,f[i]+s[i]);
     write(ans);
 end;

 begin
     init;
     work;
     get;
 end.