bzoj 1194

http://www.lydsy.com/JudgeOnline/problem.php?id=1194

我们看看怎么判断B是否是A的升级。

我们新建一个图G，每个点为一个二元组(i,j)，表示A的i号点和B的j号点。

设A中i号点的0出边为p0[i]，B中j号点的0出边为p0[j]，那么在图G中(i,j)连一条有向边到(p0[i],p0[j])

同理，设A中i号点的1出边为p1[i]，B中j号点的1出边为p1[j]，那么在图G中(i,j)连一条有向边到(p1[i],p1[j])

我们在图G中从(0,0)点开始BFS，如果到了某个点(i,j)，在A中i号是输出点，但是在B中j号点不是输出点，那么B不是A的升级。

如果没出现上面那种情况，那么B是A的升级。

现在已知道了各个咒语机的升级关系了。

如果B是A的升级，那么连有向边A->B。

我们要在这个有向图中找一条最长路径。

先缩点，变成有向无环图，然后拓扑排序。

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<fstream>
#include<algorithm>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<utility>
#include<set>
#include<bitset>
#include<vector>
#include<functional>
#include<deque>
#include<cctype>
#include<climits>
#include<complex>
//#include<bits/stdc++.h>适用于CF,UOJ，但不适用于poj

using namespace std;

typedef long long LL;
typedef double DB;
typedef pair<int,int> PII;
typedef complex<DB> CP;

#define mmst(a,v) memset(a,v,sizeof(a))
#define mmcy(a,b) memcpy(a,b,sizeof(a))
#define fill(a,l,r,v) fill(a+l,a+r+1,v)
#define re(i,a,b)  for(i=(a);i<=(b);i++)
#define red(i,a,b) for(i=(a);i>=(b);i--)
#define ire(i,x) for(typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define fi first
#define se second
#define m_p(a,b) make_pair(a,b)
#define SF scanf
#define PF printf
#define two(k) (1<<(k))

template<class T>inline T sqr(T x){return x*x;}
template<class T>inline void upmin(T &t,T tmp){if(t>tmp)t=tmp;}
template<class T>inline void upmax(T &t,T tmp){if(t<tmp)t=tmp;}

const DB EPS=1e-;
inline int sgn(DB x){if(abs(x)<EPS)return ;return(x>)?:-;}
const DB Pi=acos(-1.0);

inline int gint()
  {
        int res=;bool neg=;char z;
        for(z=getchar();z!=EOF && z!='-' && !isdigit(z);z=getchar());
        if(z==EOF)return ;
        if(z=='-'){neg=;z=getchar();}
        for(;z!=EOF && isdigit(z);res=res*+z-'',z=getchar());
        return (neg)?-res:res;
    }
inline LL gll()
  {
      LL res=;bool neg=;char z;
        for(z=getchar();z!=EOF && z!='-' && !isdigit(z);z=getchar());
        if(z==EOF)return ;
        if(z=='-'){neg=;z=getchar();}
        for(;z!=EOF && isdigit(z);res=res*+z-'',z=getchar());
        return (neg)?-res:res;
    }

const int maxN=;

int S;
struct Tdata
  {
      int N,M,flag[maxN+],p[maxN+][];
      inline void input()
        {
            int i;
            N=gint();M=gint();
            mmst(flag,);
            re(i,,M)flag[gint()+]=;
            re(i,,N)p[i][]=gint()+,p[i][]=gint()+;
        }
  }data[maxN+];

int vis[maxN+][maxN+];
int head,tail;
PII que[maxN*maxN+];
inline int check(int a,int b)
  {
      mmst(vis,);
      que[head=tail=]=PII(,);
      vis[][]=;
      if(data[a].flag[] && !data[b].flag[])return ;
      while(head<=tail)
        {
            int i,x=que[head].fi,y=que[head].se;head++;
            re(i,,)
              {
                  int tx=data[a].p[x][i],ty=data[b].p[y][i];
                  if(vis[tx][ty])continue;
                  vis[tx][ty]=;
                  que[++tail]=PII(tx,ty);
                  if(data[a].flag[tx] && !data[b].flag[ty])return ;
              }
        }
      return ;
  }

int now,first[maxN+];
struct Tedge{int v,next;}edge[maxN*maxN+];
inline void addedge(int u,int v){now++;edge[now].v=v;edge[now].next=first[u];first[u]=now;}
int tot,idx[maxN+],val[maxN+],f[maxN+],form[maxN+],ans;
vector<int> G[maxN+];

int cnt,low[maxN+],dfn[maxN+],vis2[maxN+];
int top,sta[maxN+];
inline void Tarjan(int u)
  {
      dfn[u]=low[u]=++cnt;
      vis2[sta[++top]=u]=;
      int i,v;
      for(i=first[u],v=edge[i].v;i!=-;i=edge[i].next,v=edge[i].v)
        if(!vis2[v])
          {
              Tarjan(v);
              upmin(low[u],low[v]);
          }
        else
          if(vis2[v]==) upmin(low[u],dfn[v]);
      if(low[u]==dfn[u])
        {
            ++tot;
            while()
                  {
                      idx[sta[top]]=tot;
                        val[tot]++;
                        G[tot].push_back(sta[top]);
                        vis2[sta[top]]=;
                        top--;
                      if(sta[top+]==u)break;
                    }
        }
  }

int du[maxN+];
int ge,h[maxN+];
Tedge t[maxN*maxN+];
inline void addedge2(int u,int v){ge++;t[ge].v=v;t[ge].next=h[u];h[u]=ge;}

int Q[maxN+];

int ds,out[maxN+];

int main()
  {
      /*freopen("bzoj1194.in","r",stdin);
      freopen("bzoj1194.out","w",stdout);*/
      int i,j;
      S=gint();
      re(i,,S)data[i].input();
      mmst(first,-);now=-;
      re(i,,S)re(j,,S)if(i!=j && check(i,j))addedge(i,j);/*,PF("%d->%d\n",i,j);*/
      re(i,,S)if(!vis2[i])Tarjan(i);
      ge=-;mmst(h,-);
      int u,v;
      re(u,,S)for(i=first[u],v=edge[i].v;i!=-;i=edge[i].next,v=edge[i].v)if(idx[u]!=idx[v])addedge2(idx[u],idx[v]),du[idx[v]]++;
      head=;tail=-;
      re(i,,tot)if(du[i]==)Q[++tail]=i,f[i]=val[i],form[i]=;
        while(head<=tail)
          {
              int u=Q[head++],v;
                for(i=h[u],v=t[i].v;i!=-;i=t[i].next,v=t[i].v)
                  {
                      if(val[v]+f[u]>f[v])f[v]=val[v]+f[u],form[v]=u;
                        du[v]--;
                        if(du[v]==)Q[++tail]=v;
                    }
          }
        int p;
        ans=;
        re(i,,tot)if(f[i]>ans)ans=f[i],p=i;
        cout<<ans<<endl;
        /*while(p)out[++ds]=p,p=form[p];
        red(i,ds,1)re(j,0,int(G[out[i]].size())-1)PF("%d ",G[out[i]][j]-1);PF("\n");*/
        return ;
  }