【BZOJ】【2208】【JSOI2010】连通数

题解：

　　1.Tarjan缩点以后对每个连通分量进行深搜，看能到哪些连通分量，能到达的所有连通分量的size之和记为sum。则第i个连通分量对答案的贡献为size[i]*sum（到其他连通分量）+size[i]*size[i]（本身互相可达）

　　2.在网上搜了一下……这题可以直接dfs过……汗。“正解”是Tarjan缩点+拓扑排序+状态压缩

 /**************************************************************
     Problem: 2208
     User: Tunix
     Language: C++
     Result: Accepted
     Time:9476 ms
     Memory:64772 kb
 ****************************************************************/

 //BZOJ 2208
 #include<vector>
 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 using namespace std;
 const int N=;
 typedef long long LL;
 void read(int &v){
     v=; int sign=; char ch=getchar();
     while(ch<''||ch>'') {if (ch=='-') sign=-; ch=getchar();}
     while(ch>=''&&ch<='') {v=v*+ch-''; ch=getchar();}
     v*=sign;
 }
 /*******************tamplate********************/
 int to1[N*N],head1[N],next1[N*N],cnt=;
 int to2[N*N],head2[N],next2[N*N];
 void ins1(int x,int y){
     to1[++cnt]=y; next1[cnt]=head1[x]; head1[x]=cnt;
 }
 void ins2(int x,int y){
     to2[++cnt]=y; next2[cnt]=head2[x]; head2[x]=cnt;
 }
 /***********************************************/
 int n,m;
 int dfn[N],low[N],dfs_clock=,belong[N],num,size[N];
 int st[N],top=;
 bool in[N];
 void tarjan(int x){
     int y;
     dfn[x]=low[x]=++dfs_clock;
     st[top++]=x;
     in[x]=;
     for(int i=head1[x];i;i=next1[i]){
         y=to1[i];
         if (!dfn[y]){
             tarjan(y);
             low[x]=min(low[x],low[y]);
         }
         else if (in[y]) low[x]=min(low[x],dfn[y]);
     }
     if (low[x]==dfn[x]){
         num++; size[num]=;
         for(y=;y!=x;){
             y=st[--top];
             in[y]=;
             belong[y]=num;
             size[num]++;
         }
     }
 }
 void rebuild(){
     cnt=;
     F(x,,n)
         for(int i=head1[x];i;i=next1[i])
             if (belong[x]!=belong[to1[i]])
                 ins2(belong[x],belong[to1[i]]);
 }
 /***********************************************/
 bool vis[N];
 int sum=;
 void dfs(int x){
     int y;
     for(int i=head2[x];i;i=next2[i]){
         y=to2[i];
         if (!vis[y]){
             vis[y]=;
             sum+=size[y];
             dfs(y);
         }
     }
 }
 void solve(){
     LL ans=;
     F(i,,num){
         memset(vis,,sizeof vis);
         sum=;
         dfs(i);
         ans+=size[i]*sum+size[i]*size[i];
     }
     printf("%lld\n",ans);
 }
 int main(){
     #ifndef ONLINE_JUDGE
     freopen("input.txt","r",stdin);
     #endif
     read(n);
     char s[N];
     F(i,,n){
         scanf("%s",s);
         rep(j,strlen(s))
             if (s[j]=='') ins1(i,j+);
     }
     F(i,,n) if (!dfn[i]) tarjan(i);
     rebuild();
     solve();
     return ;
 }