思路1
用warshall算法求传递闭包,之后遍历邻接矩阵计算有几个部门能够访问到N个部门
时间:O(n^3),显然超时,但是代码量短,水一点分数还是很值的。
思路2
在G中分别对所有结点进行DFS,再将所有边反向得到G^T,再次进行DFS
DFS过程中能够连接到的点记录在标记数组中,得到传递闭包+对称闭包。
时间:n*(n+m) = O(m*n),500ms,正好卡过
Warshall
有两个版本,其中第一个版本60分,第二个版本65分
#include <iostream>
using namespace std;
int n,m,a,b;
bool adj[1001][1001];
int main()
{
cin>>n>>m;
for(int i=1; i<=n; ++i)
adj[i][i] = 1;
while(m--)
{
cin>>a>>b;
adj[a][b] = 1;
}
// 版本1:
// for(int k=1; k<=n; ++k)
// for(int i=1; i<=n; ++i)
// for(int j=1; j<=n; ++j)
// adj[i][j] = adj[i][j] | (adj[i][k] & adj[k][j]);
// 版本2:
for(int j=1; j<=n; j++)
for(int i=1; i<=n; i++)
{
if(adj[i][j])
for(int k=1; k<=n; k++)
adj[i][k] |= adj[j][k];
}
int cnt = 0;
bool flag = 1;
for(int i=1; i<=n; ++i)
{
flag = 1;
for(int j=1; j<=n; ++j)
{
if(!((adj[i][j] || adj[j][i])))
{
flag = 0;
break;
}
}
if(flag) ++cnt;
}
cout<<cnt;
return 0;
}
两次DFS(100分)
#include <iostream>
#include <iomanip>
#include <cmath>
#include <cstring>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <algorithm>
#include <sstream>
#include <list>
#define endl "\n"
using namespace std;
int n,m,a,b,cur;
vector<list<int> > g1;
vector<list<int> > g2;
bool g[1005][1005];
bool vis[1005];
void dfs(int x, vector<list<int> >& G)
{
for(list<int>::iterator p = G[x].begin(); p != G[x].end(); ++p)
{
if(vis[*p] == 0)
{
vis[*p] = 1;
g[cur][*p] = 1;
dfs(*p, G);
}
}
}
inline bool isOk(int a)
{
for(int i=1;i<=n;++i)
if(g[a][i] == 0) return 0;
return 1;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
cin>>n>>m;
g1.assign(n+1, list<int>());
g2.assign(n+1, list<int>());
for(int i=0; i<m; ++i)
{
cin>>a>>b;
g1[a].push_back(b);
g2[b].push_back(a);
}
for(cur=1; cur<=n; ++cur)
{
g[cur][cur] = 1;
memset(vis,0,sizeof(vis));
vis[cur] = 1;
dfs(cur, g1);
memset(vis,0,sizeof(vis));
vis[cur] = 1;
dfs(cur, g2);
}
int cnt = 0;
for(int i=1;i<=n;++i)
if(isOk(i)) ++cnt;
cout<<cnt;
return 0;
}