时间限制: 1 s 空间限制: 128000 KB 题目等级 : 大师 Master
2007年NOI全国竞赛
题目描述 Description
在社交网络(social network)的研究中,我们常常使用图论概念去解释一些 社会现象。 不妨看这样的一个问题。在一个社交圈子里有 n 个人,人与人之间有不同程 度的关系。我们将这个关系网络对应到一个 n 个结点的无向图上,两个不同的人 若互相认识,则在他们对应的结点之间连接一条无向边,并附上一个正数权值 c, c 越小,表示两个人之间的关系越密切。 我们可以用对应结点之间的最短路长度来衡量两个人 s 和 t 之间的关系密切 程度,注意到最短路径上的其他结点为 s 和 t 的联系提供了某种便利,即这些结 点对于 s 和 t 之间的联系有一定的重要程度。我们可以通过统计经过一个结点 v 的最短路径的数目来衡量该结点在社交网络中的重要程度。 考虑到两个结点 A 和 B 之间可能会有多条最短路径。我们修改重要程度的定 义如下: 令 Cs,t表示从 s 到 t 的不同的最短路的数目, Cs,t(v)表示经过 v 从 s 到 t 的最短 路的数目;则定义I(v)为结点 v 在社交网络中的重要程度。 为了使 I(v)和 Cs,t(v)有意义,我们规定需要处理的社交网络都是连通的无向 图,即任意两个结点之间都有一条有限长度的最短路径。
现在给出这样一幅描述社交网络的加权无向图,请你求出每一个结点的重要程度。
输入描述 Input Description
输入文件中第一行有两个整数,n 和 m,表示社交网络中结点和无向边的数 目。在无向图中,我们将所有结点从 1 到 n 进行编号。 接下来 m 行,每行用三个整数 a, b, c 描述一条连接结点 a 和 b,权值为 c 的 无向边。注意任意两个结点之间最多有一条无向边相连,无向图中也不会出现自 环(即不存在一条无向边的两个端点是相同的结点)。
输出描述 Output Description
输出文件包括 n 行,每行一个实数,精确到小数点后 3 位。第 i 行的实数表 示结点 i 在社交网络中的重要程度
样例输入 Sample Input
4 4
1 2 1
2 3 1
3 4 1
4 1 1
样例输出 Sample Output
1.000
1.000
1.000
1.000
题解
——《NOI 2007 总结与解题 余林韵》
- 最短路数目 long long!double!
代码
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int MN=102;
typedef long long ll;
ll n,m,t[MN][MN],a[MN][MN];
int main(){
//fre//open("1796.in","r",stdin);
cin>>n>>m;
memset(t,0x3f,sizeof(t));
ll u,v,w;
for(int i=1;i<=m;i++){
cin>>u>>v>>w;
t[u][v]=t[v][u]=w;
a[u][v]=a[v][u]=1;
}
for(int i=1;i<=n;i++)t[i][i]=0,a[i][i]=0;
ll s,num;
for(int k=1;k<=n;k++){
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
s=t[i][k]+t[k][j];
num=a[i][k]*a[k][j];
if(s<t[i][j]){
t[i][j]=s;
a[i][j]=num;
}
else if(s==t[i][j]){
a[i][j]+=num;
}
}
}
}
//for
double I;
/*
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){cout<<t[i][j]<<' ';}cout<<endl;}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){cout<<a[i][j]<<' ';}cout<<endl;}
*/
for(int k=1;k<=n;k++){
I=0;
for(int i=1;i<=n;i++){
if(i==k)continue;
for(int j=1;j<=n;j++){
if(j==k)continue;
//cout<<i<<' '<<k<<' '<<j<<' '<<t[i][k]<<' '<<t[j][k]<<' '<<t[i][j]<<endl;
if((t[i][k]+t[k][j] == t[i][j]) && (t[i][j]!=0x3f3f3f3f)){
I+=(double)((double)(a[i][k]*a[k][j])/(double)a[i][j]);//!!!
}
}
}
printf("%.3f\n",I);
}
return 0;
}