BZOJ1491: [NOI2007]社交网络(Floyd 最短路计数)
Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 2343 Solved: 1266
[Submit][Status][Discuss]
Description
.png)
Input
Output
输出包括n行,每行一个实数,精确到小数点后3位。第i行的实数表示结点i在社交网络中的重要程度。
Sample Input
1 2 1
2 3 1
3 4 1
4 1 1
Sample Output
1.000
1.000
1.000
HINT
社交网络如下图所示。
Source
最短路计数。。
mdzz想到一个做法,应该是$N^3$的,不过与边权有关,然后被卡成$90$分卡成一下午。。
就是直接dfs求最短路计数的时候统计答案,但是不能写记忆化。会wa
// luogu-judger-enable-o2
#include<cstdio>
#include<vector>
#include<algorithm>
#include<cstring>
#include<map>
#include<cmath>
using namespace std;
const int MAXN = 1e5 + , INF = 1e8 + ;
inline int read() {
char c = getchar(); int x = , f = ;
while(c < '' || c > '') {if(c == '-') f = -; c = getchar();}
while(c >= '' && c <= '') x = x * + c - '', c = getchar();
return x * f;
}
int N, M;
int dis[][], f[][], num[][][], w[][];
void Floyed() {
for(int k = ; k <= N; k++) {
dis[k][k] = ;
for(int i = ; i <= N; i++)
for(int j = ; j <= N; j++)
dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
}
}
int GetAns(int bg, int now, int pre) {
int ans = ;
if(bg == now) return ;
for(int i = ; i <= N; i++) {
if(now != i && dis[bg][now] == dis[bg][i] + w[i][now]) {
int x = GetAns(bg, i, pre);
num[bg][pre][i] += x;
ans += x;
}
}
return ans;
}
double ans[MAXN];
main() {
#ifdef WIN32
freopen("a.in", "r", stdin);
//freopen("a.out", "w", stdout);
#endif
N = read(); M = read();
memset(dis, 0x3f, sizeof(dis));
memset(w, 0x3f, sizeof(w));
for(int i = ; i <= M; i++) {
int x = read(), y = read(), z = read();
w[x][y] = w[y][x] = dis[x][y] = dis[y][x] = z;
}
Floyed(); for(int i = ; i <= N; i++) {
f[i][i] = ;
for(int j = ; j <= N; j++)
f[i][j] = GetAns(i, j, j);
} for(int i = ; i <= N; i++) {
double ans = ;
for(int s = ; s <= N; s++) {
for(int t = ; t <= N; t++)
if(s != i && t != i && s != t)
ans += (double)num[s][t][i] / f[s][t];
}
printf("%.3lf\n", ans);
}
return ;
}
90
标算是Floyd最短路计数,
// luogu-judger-enable-o2
// luogu-judger-enable-o2
#include<cstdio>
#include<vector>
#include<algorithm>
#include<cstring>
#include<map>
#include<cmath>
#define LL long long
using namespace std;
const int INF = 1e9 + ;
inline int read() {
char c = getchar(); int x = , f = ;
while(c < '' || c > '') {if(c == '-') f = -; c = getchar();}
while(c >= '' && c <= '') x = x * + c - '', c = getchar();
return x * f;
}
int N, M;
LL dis[][], num[][];
double ans[];
void Floyed() {
for(int k = ; k <= N; k++) {
for(int i = ; i <= N; i++)
for(int j = ; j <= N; j++) {
int to = dis[i][k] + dis[k][j];
if(to == INF) continue;
if(to < dis[i][j])
dis[i][j] = to, num[i][j] = num[i][k] * num[k][j];
else if(to == dis[i][j])
num[i][j] += num[i][k] * num[k][j];
}
}
for(int i = ; i <= N; i++) {
for(int s = ; s <= N; s++) {
for(int t = ; t <= N; t++) {
if(s == i || t == i || (dis[s][i] + dis[i][t] != dis[s][t]) || s == t) continue;
ans[i] += (double)(1.0 * num[s][i] * num[i][t]) / num[s][t];
}
}
}
}
int main() {
N = read(); M = read();
for(int i = ; i <= N; i++)
for(int j = ; j <= N; j++)
dis[i][j] = INF;
for(int i = ; i <= M; i++) {
int x = read(), y = read(), z = read();
dis[x][y] = dis[y][x] = z;
num[x][y] = num[y][x] = ;
}
Floyed();
for(int i = ; i <= N; i++)
printf("%.3lf\n", ans[i]);
return ;
}
BZOJ1491: [NOI2007]社交网络(Floyd 最短路计数)的更多相关文章
- 1491. [NOI2007]社交网络【最短路计数】
Description 在社交网络(socialnetwork)的研究中,我们常常使用图论概念去解释一些社会现象.不妨看这样的一个问题. 在一个社交圈子里有n个人,人与人之间有不同程度的关系.我们将这 ...
- [BZOJ1491][NOI2007]社交网络 floyd
1491: [NOI2007]社交网络 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 2196 Solved: 1170[Submit][Status ...
- 【BZOJ1491】[NOI2007]社交网络 Floyd
[BZOJ1491][NOI2007]社交网络 Description 在社交网络(socialnetwork)的研究中,我们常常使用图论概念去解释一些社会现象.不妨看这样的一个问题. 在一个社交圈子 ...
- BZOJ 1491: [NOI2007]社交网络( floyd )
floyd...求最短路时顺便求出路径数. 时间复杂度O(N^3) ------------------------------------------------------------------ ...
- [NOI2007]社交网络(最短路)
[NOI2007]社交网络 Description 在社交网络(socialnetwork)的研究中,我们常常使用图论概念去解释一些社会现象.不妨看这样的一个问题. 在一个社交圈子里有n个人,人与人之 ...
- 洛谷P2047||bzoj1491 [NOI2007]社交网络
https://www.luogu.org/problemnew/show/P2047 https://www.lydsy.com/JudgeOnline/problem.php?id=1491 也可 ...
- BZOJ1491 [NOI2007]社交网络[最短路计数]
$n$非常的小,结合题目计算式可以想到$O(n^3)$暴枚$s,t,v$,看$v$在不在$s\to t$最短路上($dis_{s,v}+dis_{v,t}=dis_{s,v}$是$v$在两点最短路上的 ...
- BZOJ1491 [NOI2007]社交网络 【floyd】
题目 在社交网络(socialnetwork)的研究中,我们常常使用图论概念去解释一些社会现象.不妨看这样的一个问题. 在一个社交圈子里有n个人,人与人之间有不同程度的关系.我们将这个关系网络对应到一 ...
- 【Floyd】BZOJ1491: [NOI2007]社交网络
Description Solution n<=100自然联想Floyd 设两个数组d[n][n]存最短距离,t[n][n]存最短路径条数 更新d的时候顺便更新t,乘法原理 if(d[i][ ...
随机推荐
- Hadoop学习笔记(3) Hadoop文件系统二
1 查询文件系统 (1) 文件元数据:FileStatus,该类封装了文件系统中文件和目录的元数据,包括文件长度.块大小.备份.修改时间.所有者以及版权信息.FileSystem的getFileSta ...
- C# 将外部exe程序 嵌入到自己的窗体界面
将别人开发的exe程序,放到自己的窗体里面来运行. 1.基本功能实现 首先,在自己的窗体后面加上代码: [DllImport("User32.dll", EntryPoint = ...
- 自定义input文件上传 file的提示文字及样式
简单记录一下 效果图: 代码: <input class="aload" type='button' value='上传附件' onClick='javascript:$(& ...
- bundle绑定资源表
1.注册绑定资源表 在application_Start函数中: (注意不要加拓展名,否则压缩时出问题) BundleTable.Bundles.Add(new ScriptBundle(" ...
- 服务器断电后 redis重启后启动不起来
服务器断电后 redis 重启后启动不起来 原因:db持久化失败 1. 先查询redis的进程 ps -ef|grep redis 2. 查询redis的缓存文件在哪 whereis dump.rdb ...
- java webServer(一)
java webServer实现 浏览器和服务器使用的是http协议,http协议使用的是tcp 这里主要在服务器端监听端口号 实现功能 通过浏览器向服务器发送http请求:http://localh ...
- intellijidea课程 intellijidea神器使用技巧1-4 idea安装
idea安装 run intellij idea==>finish ==>do not import settings==>ok==>evaluate for free==&g ...
- Brackets - 前端神器
做了几年的 .Net 项目开发,后来公司转 Java 语言开发,Java 做了还没一年,公司准备前后端分离开发,而我被分到前端! Brackets是一款基于web(html+css+js)开发的web ...
- 视频SharePoint 2010 大局观 By 陈希章[zt]
SharePoint 2010 大局观 By 陈希章 http://hi.baidu.com/jinzesudawei/item/9d0adace8cbcaf2c47d5c0fc
- Anaconda上安装Tensorflow并在jupyter上运行
博客原文地址:https://blog.csdn.net/index20001/article/details/73555182 https://www.cnblogs.com/HongjianChe ...