COGS 449. 网络病毒

★★   输入文件：virus.in   输出文件：virus.out   简单对比
时间限制：1 s   内存限制：128 MB

【题目描述】 
公元2008年10月31日星期五，笃志者所在的整个机房由于猖獗的病毒一片恐慌。经查证，病毒是由A1机器散播开来的。。这要追溯到29日，笃志者由于病毒被迫从A1机器撤离。 
一想到病毒是从自己的机器传开的，笃志者就心神不宁。他决定搞清楚病毒是怎么散播开来的。事实上，机房内的机器并不是全部都能够互相感染的。笃志者（ceeji）好不容易经过测试得到了机房中各机器间是否连通的图表，就在他马上就要得出结果的时候，大脑突然乱了！问题的严重性在于：如果他不在1s内搞清楚这个问题，机房就会整体瘫痪。现在笃志者求助于你，他需要知道病毒从未感染机房开始，最少入侵几台机器之后，机房就会整体感染。 
【输入格式】 
文件的第一行为一个整数n，第二行至第n+1行为n*n的矩阵（若第i行第j列为1，则机器i能对机器j进行ARP攻击（即感染机器j），若第i行第j列为0，则机器i不能感染机器j）。 
文件名为“2.in”。 
【输出格式】 
输出文件只有一行，为笃志者想知道的最少感染机器数。 
文件名为“2.out”。 
【输入样例】 
8 
0 0 1 0 0 0 0 0 
1 0 0 1 0 0 0 0 
0 1 0 1 1 0 0 0 
0 0 0 0 0 1 0 0 
0 0 0 1 0 0 0 0 
0 0 0 1 0 0 0 0 
0 0 0 1 0 0 0 1 
0 0 0 0 0 0 1 0 
【输出样例】 
2 
【数据范围】 
对于 100% 的数据，n<=1000

求入度为0的强连通分量的个数

一开始手残了。。顺道打上了拓扑排序

屠龙宝刀点击就送

#include <cstdio>
#include <vector>
#define N 10000
using namespace std;
vector<int>G[N];
bool instack[N];
int n,tim,cnt,top,ans,sumcol,c[N],in[N],to[N<<5],low[N],dfn[N],col[N],stack[N],head[N],nextt[N<<5];
inline int min(int a,int b) {return a>b?b:a;}
inline int max(int a,int b) {return a>b?a:b;}
void tarjan(int x)
{
    stack[++top]=x;
    low[x]=dfn[x]=++tim;
    instack[x]=true;
    for(int i=head[x];i;i=nextt[i])
    {
        int v=to[i];
        if(!dfn[v])
        {
            tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if(instack[v]) low[x]=min(low[x],dfn[v]);
    }
    if(low[x]==dfn[x])
    {
        sumcol++;
        int k;
        do{k=stack[top--];instack[k]=false;col[k]=sumcol;} while(k!=x);
    }
}
int main(int argc,char *argv[])
{
    freopen("virus.in","r",stdin);
    freopen("virus.out","w",stdout);
    scanf("%d",&n);
    for(int i=1;i<=n;++i)
     for(int a,j=1;j<=n;++j)
     {
         scanf("%d",&a);
        if(a) {nextt[++cnt]=head[i];to[cnt]=j;head[i]=cnt;} 
     }
    for(int i=1;i<=n;++i) if(!dfn[i]) tarjan(i);
    for(int i=1;i<=n;++i)
    {
        for(int j=head[i];j;j=nextt[j])
        {
            int v=to[j];
            if(col[i]==col[v]) continue;
            for(int k=0;k<G[col[i]].size();++k)
            if(G[col[i]][k]==col[v]) goto flag;
            G[col[i]].push_back(col[v]);
            in[col[v]]++; 
            flag:;
        }
    }
    for(int i=1;i<=sumcol;++i) if(!in[i]) ans++;
    printf("%d
",ans);
    return 0;
}