7-12（图） 社交网络图中结点的“重要性”计算（30 分）

在社交网络中，个人或单位（结点）之间通过某些关系（边）联系起来。他们受到这些关系的影响，这种影响可以理解为网络中相互连接的结点之间蔓延的一种相互作用，可以增强也可以减弱。而结点根据其所处的位置不同，其在网络中体现的重要性也不尽相同。

“紧密度中心性”是用来衡量一个结点到达其它结点的“快慢”的指标，即一个有较高中心性的结点比有较低中心性的结点能够更快地（平均意义下）到达网络中的其它结点，因而在该网络的传播过程中有更重要的价值。在有N个结点的网络中，结点v​i​​的“紧密度中心性”Cc(v​i​​)数学上定义为v​i​​到其余所有结点v​j​​ (j≠i) 的最短距离d(v​i​​,v​j​​)的平均值的倒数：

对于非连通图，所有结点的紧密度中心性都是0。

给定一个无权的无向图以及其中的一组结点，计算这组结点中每个结点的紧密度中心性。

输入格式:

输入第一行给出两个正整数N和M，其中N（≤10​4​​）是图中结点个数，顺便假设结点从1到N编号；M（≤10​5​​）是边的条数。随后的M行中，每行给出一条边的信息，即该边连接的两个结点编号，中间用空格分隔。最后一行给出需要计算紧密度中心性的这组结点的个数K（≤100）以及K个结点编号，用空格分隔。

输出格式:

按照Cc(i)=x.xx的格式输出K个给定结点的紧密度中心性，每个输出占一行，结果保留到小数点后2位。

输入样例:

9 14
1 2
1 3
1 4
2 3
3 4
4 5
4 6
5 6
5 7
5 8
6 7
6 8
7 8
7 9
3 3 4 9

输出样例:

Cc(3)=0.47
Cc(4)=0.62
Cc(9)=0.35
解析：这个题如果掌握迪杰斯特拉算法并不难，希望同学们好好看课本，
我还是变参考课本边做的，有个难点就是，你如果想偷懒，静态分配一个
a【n】【n】的数组，那么肯定会数组太大直接崩掉，所以动态分配就好了，
（希望仅供参考，不要直接复制到作业）
代码：

#include<iostream>
#include<cstring>
#include<string>
#include<bits/stdc++.h>
#include<cstdio>
#include<cmath>
#define maxn 10001
using namespace std;
int n,m;
int d[maxn],visit[maxn];
int Dijkstra(int *a[],int v){
    for(int i=1;i<=n;i++){

        visit[i]=0;
        d[i]=a[v][i];

}

visit[v]=1;
 d[v]=0;
     int sum=0;
     int x;
    for(int k=1;k<=n;k++){
            int min=maxn;
     for(int w=1;w<=n;w++)
       if(!visit[w]&&d[w]<min){
        x=w;min=d[w];
       }
    visit[x]=1;
    for(int w=1;w<=n;w++){
        if(!visit[w]&&d[x]+a[x][w]<d[w]){
            d[w]=d[x]+a[x][w];
        }
    }

}for(int i=1;i<=n;i++)
  {
    if(d[i]==maxn)return -1;
    sum+=d[i];
    //cout<<"@@@@@@@@@"<<d[i]<<endl;
    }

return sum;
}

int main(){
int x,y;
scanf("%d%d",&n,&m);
int **a=(int**)malloc(sizeof(int*)*(n+1));
 for(int i=1; i<=n; i++) {
       a[i]=(int*)malloc(sizeof(int)*(n+1));

 }


for(int i=1;i<=n;i++){
    d[i]=maxn;
    for(int j=1;j<=n;j++){
        a[i][j]=maxn;
    }
}

//输入
for(int i=1;i<=m;i++){
        cin>>x>>y;
        a[x][y]=a[y][x]=1;
}
int k;
cin>>k;int v;
for(int i=0;i<k;i++){
    cin>>v;
  int su=  Dijkstra(a,v);
  if(su==-1) printf("Cc(%d)=%.2lf\n",v,0.00);
  else
  printf("Cc(%d)=%.2lf\n",v,(n-1)*1.0/su*1.0);
}

return 0;


}