洛谷 P1262 间谍网络 强连通分量 缩点 tarjan bfs

题目链接：

https://www.luogu.org/problem/P1262

说明：

特别奇怪，自己的代码确实自己感觉思路正确，提交也ac，但是就是在cb上运行样例2时，不通过，可能是此题样例给错了

思路：

1：tarjan缩点

2：经过步骤1缩点之后，每一个节点就代表一个强连通分量，此时节点之间都是链式结构，没有环

3：首先初步处理可以被策反的间谍中所属的强连通分量，如果一个强连通分量中有多个可以被策反的间谍，则取花费最少的那一个花费

int en[cnt+1][2];
    for(int i=1;i<=cnt;i++)
    {
        en[i][1]=1e9;
    }
    for(int i=1;i<=p;i++)
    {
        en[color[money[i][0]]][1]=min(en[color[money[i][0]]][1],money[i][1]);
    }

en[i][0]存强连通分量编号即重构后的图的节点编号，en[i][1]存这个强连通分量可以被策反的最少费用

4：关于重构后的图，要明白，如果要使每一个节点都可以被遍历到，那么就必须保证入度为0的所有节点都可以被策反

     1：如果入度为0的所有节点都可以被策反，那么策反这些节点的费用和即为最少花费

     2：只要发现有一个入度为0的节点不可以被策反，那么就意味着，不可以保证所有节点都可以被策反了，那么此时就去对所有可以被策反的节点运行bfs，然后顺序遍历所有节点，发现第一个所属强连通分量没有被遍历过的节点就是不能控制的间谍中，编号最小的间谍编号

int sum=0;
    for(int i=1;i<=cnt;i++)
    {
        if(id[i]==0)
        {
            if(en[i][1]!=1e9)
            {
                bfs(i);
                sum+=en[i][1];
            }
            else
            {
                for(int j=1;j<=cnt;j++)
                {
                    if(en[j][1]!=1e9)
                    {
                        bfs(j);
                    }
                }
                for(int j=1;j<=n;j++)
                {
                         if(colorcnt[color[j]]==0)
                         {
                             cout<<"NO"<<endl<<j<<endl;
                             break;
                         }
                }
                return 0;
            }
        }
    }
    cout<<"YES"<<endl<<sum<<endl;
    return 0;

#include <bits/stdc++.h>

using namespace std;
const int maxn=3e3+1;
vector<int>e[maxn];
set<int>ee[maxn];
int dfn[maxn],low[maxn],ins[maxn],color[maxn],timing,cnt,n,p,m,id[maxn],colorcnt[maxn];
stack<int>s;
vector<pair<int,int> >mapp[maxn];
int ing[maxn];
int money[maxn][2];

void tarjan(int x)
{
    low[x]=dfn[x]=++timing;
    s.push(x);
    ins[x]=1;
    for(int i=0;i<e[x].size();i++)
    {
        int v=e[x][i];
        if(!dfn[v])
        {
            tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if(ins[v])
        {
            low[x]=min(low[x],dfn[v]);
        }
    }
    if(dfn[x]==low[x])
    {
        cnt++;
        while(s.top()!=x)
        {
            color[s.top()]=cnt;
            ins[s.top()]=0;
            s.pop();
        }
        color[s.top()]=cnt;
        ins[s.top()]=0;
        s.pop();
    }
}

int bfs(int x)
{
    queue<int>q;
    q.push(x);
    ing[x]=1;
    while(!q.empty())
    {
        int now=q.front();
        colorcnt[now]=1;
        q.pop();
        ing[now]=0;
        for(int i=0;i<mapp[now].size();i++)
        {
                int v=mapp[now][i].first;
                if(ing[v])
                    continue;
                q.push(v);
                ing[v]=1;
        }
    }
}

int main()
{
    ios::sync_with_stdio(0);
    cin>>n>>p;
    for(int i=1;i<=p;i++)
    {
        cin>>money[i][0]>>money[i][1];
    }
    cin>>m;
    for(int i=1;i<=m;i++)
    {
        int a,b;
        cin>>a>>b;
        e[a].push_back(b);
    }
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i])
            tarjan(i);
    }
    int en[cnt+1][2];
    for(int i=1;i<=cnt;i++)
    {
        en[i][1]=1e9;
    }
    for(int i=1;i<=p;i++)
    {
        en[color[money[i][0]]][1]=min(en[color[money[i][0]]][1],money[i][1]);
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<e[i].size();j++)
        {
            if(color[i]!=color[e[i][j]]&&ee[color[i]].find(color[e[i][j]])==ee[color[i]].end())
            {
                ee[color[i]].insert(color[e[i][j]]);
                id[color[e[i][j]]]++;
            }
        }
    }
    int sum=0;
    for(int i=1;i<=cnt;i++)
    {
        if(id[i]==0)
        {
            if(en[i][1]!=1e9)
            {
                bfs(i);
                sum+=en[i][1];
            }
            else
            {
                for(int j=1;j<=cnt;j++)
                {
                    if(en[j][1]!=1e9)
                    {
                        bfs(j);
                    }
                }
                for(int j=1;j<=n;j++)
                {
                         if(colorcnt[color[j]]==0)
                         {
                             cout<<"NO"<<endl<<j<<endl;
                             break;
                         }
                }
                return 0;
            }
        }
    }
    cout<<"YES"<<endl<<sum<<endl;
    return 0;
}

补充：

1：以上对思路的分析都正确，但是确实自己的代码是不正确的，但是不知道洛谷为什么就ac了，整个代码少了一行，加上这一行才是真正的正确

if(color[i]!=color[e[i][j]]&&ee[color[i]].find(color[e[i][j]])==ee[color[i]].end())
            {
                ee[color[i]].insert(color[e[i][j]]);
                id[color[e[i][j]]]++;
                mapp[color[i]].push_back(make_pair(color[e[i][j]],0));
            }

#include <bits/stdc++.h>

using namespace std;
const int maxn=3e3+1;
vector<int>e[maxn];
set<int>ee[maxn];
int dfn[maxn],low[maxn],ins[maxn],color[maxn],timing,cnt,n,p,m,id[maxn],colorcnt[maxn];
stack<int>s;
vector<pair<int,int> >mapp[maxn];
int ing[maxn];
int money[maxn][2];

void tarjan(int x)
{
    low[x]=dfn[x]=++timing;
    s.push(x);
    ins[x]=1;
    for(int i=0;i<e[x].size();i++)
    {
        int v=e[x][i];
        if(!dfn[v])
        {
            tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if(ins[v])
        {
            low[x]=min(low[x],dfn[v]);
        }
    }
    if(dfn[x]==low[x])
    {
        cnt++;
        while(s.top()!=x)
        {
            color[s.top()]=cnt;
            ins[s.top()]=0;
            s.pop();
        }
        color[s.top()]=cnt;
        ins[s.top()]=0;
        s.pop();
    }
}

void bfs(int x)
{
    queue<int>q;
    q.push(x);
    ing[x]=1;
    while(!q.empty())
    {
        int now=q.front();
        colorcnt[now]=1;
        q.pop();
        ing[now]=0;
        for(int i=0;i<mapp[now].size();i++)
        {
                int v=mapp[now][i].first;
                if(ing[v])
                    continue;
                q.push(v);
                ing[v]=1;
        }
    }
}

int main()
{
    ios::sync_with_stdio(0);
    cin>>n>>p;
    for(int i=1;i<=p;i++)
    {
        cin>>money[i][0]>>money[i][1];
    }
    cin>>m;
    for(int i=1;i<=m;i++)
    {
        int a,b;
        cin>>a>>b;
        e[a].push_back(b);
    }
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i])
            tarjan(i);
    }
    int en[cnt+1][2];
    for(int i=1;i<=cnt;i++)
    {
        en[i][1]=1e9;
    }
    for(int i=1;i<=p;i++)
    {
        en[color[money[i][0]]][1]=min(en[color[money[i][0]]][1],money[i][1]);
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<e[i].size();j++)
        {
            if(color[i]!=color[e[i][j]]&&ee[color[i]].find(color[e[i][j]])==ee[color[i]].end())
            {
                ee[color[i]].insert(color[e[i][j]]);
                id[color[e[i][j]]]++;
                mapp[color[i]].push_back(make_pair(color[e[i][j]],0));
            }
        }
    }
    int sum=0;
    for(int i=1;i<=cnt;i++)
    {
        if(id[i]==0)
        {
            if(en[i][1]!=1e9)
            {
                bfs(i);
                sum+=en[i][1];
            }
            else
            {
                for(int j=1;j<=cnt;j++)
                {
                    if(en[j][1]!=1e9)
                    {
                        bfs(j);
                    }
                }
                for(int j=1;j<=n;j++)
                {
                         if(colorcnt[color[j]]==0)
                         {
                             cout<<"NO"<<endl<<j<<endl;
                             break;
                         }
                }
                return 0;
            }
        }
    }
    cout<<"YES"<<endl<<sum<<endl;
    return 0;
}