Morass

Largest Clique

Mar 26th, 2016
145
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.78 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define PB push_back
  4. #define ZERO (1e-10)
  5. #define INF (1<<29)
  6. #define CL(A,I) (memset(A,I,sizeof(A)))
  7. #define DEB printf("DEB!\n");
  8. #define D(X) cout<<"  "<<#X": "<<X<<endl;
  9. #define EQ(A,B) (A+ZERO>B&&A-ZERO<B)
  10. typedef long long ll;
  11. typedef long double ld;
  12. typedef pair<ll,ll> pll;
  13. typedef vector<int> vi;
  14. typedef pair<int,int> ii;
  15. typedef vector<ii> vii;
  16. #define IN(n) int n;scanf("%d",&n);
  17. #define FOR(i, m, n) for (int i(m); i < n; i++)
  18. #define REP(i, n) FOR(i, 0, n)
  19. #define F(n) REP(i, n)
  20. #define FF(n) REP(j, n)
  21. #define FT(m, n) FOR(k, m, n)
  22. #define aa first
  23. #define bb second
  24. void ga(int N,int *A){F(N)scanf("%d",A+i);}
  25. #define MX (1024)
  26. vi g[MX],G[MX],e[MX],R[MX];
  27. int cn[MX],C[MX],I[MX],H[MX],L,S[MX];
  28. void CLR(int N=MX){CL(S,L=0);F(N)g[i].clear(),e[i].clear(),G[i].clear(),R[i].clear(),cn[i]=I[i]=0;}
  29. void dfs1(int x){
  30.     cn[x]=1;
  31.     for(auto h:g[x])if(!cn[h])dfs1(h);
  32.     H[L++]=x;
  33. }
  34. void dfs2(int x,int c){
  35.     cn[x]=0,C[x]=c;
  36.     for(auto h:e[x])if(cn[h])dfs2(h,c);
  37. }
  38. #define ADD(A,B) (g[A].PB(B),e[B].PB(A))
  39. int ini(int N){
  40.     F(N)if(!cn[i])dfs1(i);
  41.     int I(0);
  42.     for(int i(L-1);~i;--i)if(cn[H[i]])dfs2(H[i],++I);
  43.     F(N)++S[C[i]-1];
  44.     F(N)for(auto&h:g[i])if(C[i]!=C[h])G[C[i]-1].PB(C[h]-1),R[C[h]-1].PB(C[i]-1);
  45.     F(I)sort(G[i].begin(),G[i].end()),unique(G[i].begin(),G[i].end());//+/-
  46.     return I;
  47. }
  48. int N,M,a,b,X,dp[MX];
  49. int dfs(int u){
  50.     int &v(dp[u]);
  51.     if(~v)return v;
  52.     v=S[u];
  53.     for(auto&h:G[u])v=max(v,dfs(h)+S[u]);
  54.     return v;
  55. }
  56. int main(void){
  57.     IN(tt)F(tt){
  58.         scanf("%d%d",&N,&M);
  59.         CLR(N),X=0,CL(dp,-1);
  60.         F(M)scanf("%d%d",&a,&b),--a,--b,ADD(a,b);
  61.         N=ini(N);
  62.         F(N)X=max(X,dfs(i));
  63.         printf("%d\n",X);
  64.     }
  65.     return 0;
  66. }
Advertisement
Add Comment
Please, Sign In to add comment