bzoj 1093: [ZJOI2007]最大半连通子图

2016-06-21

读一下题目，发现环内肯定在一个半联通子图内，于是先tarjan缩点，然后就是找这个图的最长链及数目。拓扑一下即可。

注意注意  此题一定要将边进行判重。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstdlib>
#include<map>
#define ll long long
#define M 2000009
using namespace std;
ll read()
{
    char ch=getchar();
    ll x=0,f=1;
    for(;ch<'0'||ch>'9';ch=getchar())
      if(ch=='-')
        f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar())
      x=x*10+ch-'0';
    return x*f;
}
int du[M],dfn[M],low[M],K,n,m,tot,head[M],next[M],u[M],belong[M],cnt,T,q[M],sum,he[M];
bool f[M];
int ss;
struct data
{
    int d,s;
}a[M],qq[M];
bool cmp(data a1,data a2)
{
    if(a1.s==a2.s)
      return a1.d<a2.d;
    return a1.s<a2.s;
}
void jia(int a1,int a2)
{
    next[++cnt]=head[a1];
    head[a1]=cnt;
    u[cnt]=a2;
}
void tarjin(int x)
{
    dfn[x]=low[x]=++T;
    q[++tot]=x;
    f[x]=1;
    for(int i=head[x];i;i=next[i])
      if(!dfn[u[i]])
        {
            tarjin(u[i]);
            low[x]=min(low[x],low[u[i]]);
        }
      else  if(f[u[i]])
             low[x]=min(low[x],dfn[u[i]]);
    if(low[x]==dfn[x])
      {
          sum++;
          for(;q[tot]!=x;tot--)
            {
                he[sum]++;
                belong[q[tot]]=sum;
                f[q[tot]]=0;
            }
        he[sum]++;
        belong[x]=sum;
        f[x]=0;
        tot--;
      }
}
int main()
{
    n=read();
    m=read();
    K=read();
    for(int i=1;i<=m;i++)
      {
          int a=read(),b=read();
          jia(a,b);
      }
    sum=n;
    for(int i=1;i<=n;i++)
      if(!dfn[i])
        tarjin(i);
    for(int i=1;i<=n;i++)
      for(int j=head[i];j;j=next[j])
          {
            qq[++ss].s=belong[i];
            qq[ss].d=belong[u[j]];
        }
    sort(qq+1,qq+ss+1,cmp);
    for(int i=1;i<=ss;i++)
      if(qq[i].s!=qq[i].d&&(qq[i].s!=qq[i-1].s||qq[i].d!=qq[i-1].d))
        {
          jia(qq[i].s,qq[i].d);
          du[qq[i].d]++;
        }
    int h=0,t=0;
    for(int i=n+1;i<=sum;i++)
      if(!du[i])
        {
            q[++t]=i;
            a[i].d=he[i];
            a[i].s=1;
        }
    for(;h<t;)
      {
          int p=q[++h];
          for(int i=head[p];i;i=next[i])
            {
                du[u[i]]--;
                if(!du[u[i]])
                  q[++t]=u[i];
                if(a[u[i]].d==a[p].d+he[u[i]])
                  a[u[i]].s=(a[u[i]].s+a[p].s)%K;
            if(a[u[i]].d<a[p].d+he[u[i]])
                  {
                    a[u[i]].d=a[p].d+he[u[i]];
                a[u[i]].s=a[p].s;
              }
          }
      }
    int mx=0,mxx=0;
    for(int i=n+1;i<=sum;i++)
      {
        if(a[i].d==mx)
          mxx=(mxx+a[i].s)%K;
        if(a[i].d>mx)
          {
              mx=a[i].d;
              mxx=a[i].s;
          }
      }
    printf("%d
%d",mx,mxx);
    return 0;
}