hdu5739

以前从来没写过求点双连通分量，现在写一下……

这题还用到了一个叫做block forest data structure，第一次见过……

——对于每一个点双联通分量S, 新建一个节点s, 向S中每个节点v连边. 这样一来, 新增的点和原来图中的点会构成一个森林

主要是这个性质很厉害：

【对于这个森林F, 删掉一个关键点或者一个叶子ii之后, 会得到一个新森林Fi, 这个Fi​​对应的连通块集合和Gi对应的连通块集合其实是一样的(不考虑那些新增的点)】

更具体的题解见http://www.cnblogs.com/WABoss/p/5696926.html

顺便学了一个扩大栈空间的方法：在编译器中加入命令--stack=16777216 （16777216相当于16MB空间）

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<vector>
 #include<stack>
 typedef long long ll;
 using namespace std;
 const int mo=;
 struct edge{int x,y;};
 vector<int> g[];
 stack<edge> st;
 int dfn[],low[],f[],w[],s[],a[],root[],be[],fa[];
 bool cut[];
 int cas,n,m,t,h,x,y,sum;

 ll quick(ll x)
 {
   int y=mo-; ll s=;
   while (y)
   {
     if (y&) s=s*x%mo;
     x=x*x%mo;
     y>>=;
   }
   return s;
 }

 void dfs(int x)
 {
   dfn[x]=low[x]=++h;
   int chi=;
   for (int i=;i<g[x].size();i++)
   {
     int y=g[x][i];
     if (!dfn[y])
     {
       st.push((edge){x,i});
       fa[y]=x;dfs(y);
       chi++;
       low[x]=min(low[x],low[y]);
       if (low[y]>=dfn[x])
       {
         cut[x]=;
         g[++t].clear();
         while ()
         {
           int u=st.top().x,j=st.top().y;
           int v=g[u][j]; st.pop();
           if (be[u]!=t) {g[t].push_back(u); be[u]=t;}
           if (be[v]!=t) {g[t].push_back(v); be[v]=t;}
           if (u==x&&j==i) break;
         }
       }
     }
     else if (dfn[y]<dfn[x]&&fa[x]!=y)
     {
       st.push((edge){x,i});
       low[x]=min(low[x],dfn[y]);
     }
   }
   if (fa<&&chi==) cut[x]=;
 }

 void dp(int x)
 {
   dfn[x]=root[h];
   if (x<=n) w[x]=a[x]; else w[x]=;
   s[x]=;
   for (int i=; i<g[x].size(); i++)
   {
     int y=g[x][i];
     if (!dfn[y])
     {
       dp(y);
       w[x]=1ll*w[x]*w[y]%mo;
       s[x]=(s[x]+w[y])%mo;
     }
   }
 }

 int main()
 {
   freopen("1006.in","r",stdin);
   freopen("1006.ans","w",stdout);
   scanf("%d",&cas);
   while (cas--)
   {
     scanf("%d%d",&n,&m);
     for (int i=; i<=n; i++)
     {
       dfn[i]=fa[i]=cut[i]=be[i]=;
       scanf("%d",&a[i]);
     }
     for (int i=; i<=m; i++)
     {
       scanf("%d%d",&x,&y);
       g[x].push_back(y);
       g[y].push_back(x);
     }
     t=n;
     for (int i=; i<=n; i++)
       if (!dfn[i])
       {
         h=;
         dfs(i);
       }
     for (int i=; i<=n; i++) g[i].clear();
     for (int i=n+; i<=t; i++)
       for (int j=; j<g[i].size(); j++) {int x=g[i][j]; g[x].push_back(i);}
     for (int i=; i<=t; i++) dfn[i]=f[i]=;
     h=sum=;
     for (int i=; i<=t; i++)
       if (!dfn[i])
       {
         root[++h]=i;
         dp(i);
         sum=(sum+w[i])%mo;
       }
     for (int i=; i<=n; i++)
     {
       if (dfn[i]==i) f[i]=; else f[i]=1ll*w[dfn[i]]*quick(w[i])%mo;
       f[i]=((ll)f[i]+sum-w[dfn[i]]+s[i]+mo)%mo;
     }
     int ans=;
     for (int i=; i<=n; i++) ans=(ans+1ll*f[i]*i%mo)%mo;
     printf("%d\n",ans);
     for (int i=; i<=t; i++) g[i].clear();
   }
   return ;
 }