题目链接: CodeForces - 402E
题目大意
给你一个n*n的矩阵, 每个元素>=0, 可以对这个矩阵做任意次方的乘方操作, 求能不能的到一个所有元素>0的矩阵
思路
将这个矩阵看作图的邻接矩阵, 非零为有边, 那么对这个矩阵做乘方就相当于离散数学里面求传递闭包的方法嘛…最后的到的矩阵 aij 为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