Strictly Positive Matrix CodeForces - 402E 强联通分量 模版

题目链接: CodeForces - 402E

题目大意

给你一个n*n的矩阵, 每个元素>=0, 可以对这个矩阵做任意次方的乘方操作, 求能不能的到一个所有元素>0的矩阵

思路

将这个矩阵看作图的邻接矩阵, 非零为有边, 那么对这个矩阵做乘方就相当于离散数学里面求传递闭包的方法嘛…最后的到的矩阵$a_{ij}$为1表示原来图中ij有路径相连, 要整个矩阵都为1, 说明原图是强联通图, 也就是强联通分量为1

代码

强联通分量(矩阵)

效率$O(|V|+|E|)$
在CF上耗时800ms, 原来用的cin直接超时了, 改成scanf才勉强过得

#include <bits/stdc++.h>
using namespace std;
const int maxv = 2e3+10;
int V;
bool G[maxv][maxv], rG[maxv][maxv], used[maxv];
int vs[maxv], tot;
int cmp[maxv];

void add_edge(int from, int to)
{
    G[from][to] = 1;
    rG[to][from] = 1;
}

void dfs(int v)
{
    used[v] = true;
    for(int i=0; i<V; ++i)
    {
        if(G[v][i] && !used[i]) dfs(i);
    }
    vs[tot++] = v;
}

void rdfs(int v, int k)
{
    used[v] = true;
    cmp[v] = k;
    for(int i=0; i<V; ++i)
    {
        if(rG[v][i] && !used[i]) rdfs(i, k);
    }
}

int scc()
{
    for(int v=0; v<V; ++v) if(!used[v]) dfs(v);
    memset(used, 0, sizeof(bool)*V);
    int k=0;
    for(int i=tot-1; i>=0; --i)
        if(!used[vs[i]]) rdfs(vs[i], k++);
    return k;
}

int main()  
{
    scanf("%d", &V);
    int t;
    for(int i=0; i<V; ++i)
    {
        for(int j=0; j<V; ++j)
        {
            scanf("%d", &t);
            if(t) add_edge(i, j);
        }
    }
    if(scc() == 1) puts("YES");
    else puts("NO");

    return 0;
}

邻接矩阵

#include <bits/stdc++.h>
using namespace std;

const int maxv = 1e5+100;
int V;
vector<int> G[maxv], rG[maxv];
vector<int> vs;
bool used[maxv];
int cmp[maxv];

void add_edge(int from, int to)
{
    G[from].push_back(to);
    rG[to].push_back(from);
}
void dfs(int v)
{
    used[v] = true;
    for(auto ite : G[v])
    {
        if(!used[ite]) dfs(ite);
    }
    vs.push_back(v);
}
void rdfs(int v, int k)
{
    used[v] = true;
    cmp[v] = k;
    for(auto ite : rG[v])
    {
        if(!used[ite]) rdfs(ite, k);
    }
}

int scc()
{
    memset(used, 0, sizeof(used));
    vs.clear();
    for(int v=0; v<V; ++v) if(!used[v]) dfs(v);
    memset(used, 0, sizeof(used));
    int k = 0;
    for(int i=vs.size()-1; i>=0; --i)
        if(!used[vs[i]]) rdfs(vs[i], k++);
    return k;
}


int main()
{
    int E;
    cin >> V >> E;
    for(int i=0; i<E; ++i)
    {
        int u, v;
        cin >> u >> v;
        --u, --v;
        add_edge(u, v);
    }
    cout << scc() << endl;

    return 0;
}
//Input:
// 12 16
// 12 11
// 11 10
// 11 8
// 8 10
// 10 9
// 9 8
// 9 7
// 5 7
// 7 6
// 6 5
// 6 3
// 6 4
// 4 3
// 3 2
// 2 3
// 4 1
//Output:
//7