【BZOJ5109】【CodePlus2017】大吉大利，晚上吃鸡！

【题目链接】

点击打开链接

【思路要点】

补档博客，无题解。

【代码】

#include<bits/stdc++.h>
using namespace std;
#define MAXN    50005
#define INF 1e18
#define P   1000000007
template <typename T> void read(T &x) {x = 0; int f = 1;char c = getchar();for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';x *= f;
}
struct info {long long value; int pos; };
bool operator < (info a, info b) {return a.value > b.value; }
int n, m, s, t;
int x[MAXN], y[MAXN], z[MAXN];
vector <int> a[MAXN];
vector <int> b[MAXN];
priority_queue <info> Heap;
long long dist[MAXN];
bool finished[MAXN];
bitset <MAXN> seta[MAXN], setb[MAXN];
map <int, bitset <MAXN> > f;
long long cnta[MAXN], cntb[MAXN];
void worka(int pos) {cnta[pos] = (pos == t);seta[pos].set(pos);if (pos == t) return;for (unsigned i = 0; i < a[pos].size(); i++) {worka(a[pos][i]);cnta[pos] += cnta[a[pos][i]];cnta[pos] %= P;if (cnta[a[pos][i]]) seta[pos] |= seta[a[pos][i]];}
}
void workb(int pos) {cntb[pos] = (pos == s);setb[pos].set(pos);if (pos == s) return;for (unsigned i = 0; i < b[pos].size(); i++) {workb(b[pos][i]);cntb[pos] += cntb[b[pos][i]];cntb[pos] %= P;if (cntb[b[pos][i]]) setb[pos] |= setb[b[pos][i]];}
}
int main() {read(n), read(m), read(s), read(t);for (int i = 1; i <= m; i++) {read(x[i]), read(y[i]), read(z[i]);if (z[i] == 0) continue;a[x[i]].push_back(y[i]);a[y[i]].push_back(x[i]);b[x[i]].push_back(z[i]);b[y[i]].push_back(z[i]);}for (int i = 1; i <= n; i++)dist[i] = INF;dist[s] = 0;Heap.push((info) {0, s});while (!Heap.empty()) {while (!Heap.empty() && finished[Heap.top().pos])Heap.pop();if (Heap.empty()) break;info tmp = Heap.top();Heap.pop();for (unsigned i = 0; i < a[tmp.pos].size(); i++)if (tmp.value + b[tmp.pos][i] < dist[a[tmp.pos][i]]) {dist[a[tmp.pos][i]] = tmp.value + b[tmp.pos][i];Heap.push((info) {dist[a[tmp.pos][i]], a[tmp.pos][i]});}}for (int i = 1; i <= n; i++)a[i].clear();for (int i = 1; i <= n; i++)b[i].clear();for (int i = 1; i <= m; i++) {if (z[i] == 0) continue;if (dist[x[i]] + z[i] == dist[y[i]]) a[x[i]].push_back(y[i]), b[y[i]].push_back(x[i]);if (dist[y[i]] + z[i] == dist[x[i]]) a[y[i]].push_back(x[i]), b[x[i]].push_back(y[i]);}worka(s);workb(t);for (int i = 1; i <= n; i++) {cnta[i] = cnta[i] * cntb[i] % P;seta[i] |= setb[i];seta[i].set(i);f[cnta[i]].set(i);}long long ans = 0;for (int i = 1; i <= n; i++) {int tmp = (cnta[t] - cnta[i] + P) % P;ans += (f[tmp] & (~seta[i])).count();}cout << ans / 2 << endl;return 0;
}

【BZOJ5109】【CodePlus2017】大吉大利，晚上吃鸡！相关推荐

HTML5+css+JS实现页面打枪声音特效减压玩法 PUBG 大吉大利今晚吃鸡
HTML5+css+JS实现页面打枪声音特效减压玩法 PUBG 大吉大利今晚吃鸡 css部分 *{margin: 0;padding:0;}body{min-height:100vh;overflo ...
老司机带你在MySQL领域“大吉大利，晚上吃鸡”
关注↑↑↑我们获得更多精彩内容! 作者 | 张甦, 数据库领域的专家和知名人士.图书<MySQL王者晋级之路>作者,51CTO 专家博主.近10年互联网线上处理及培训经验,专注于 MySQ ...
大吉大利-今晚吃鸡-枪械
题目描述在绝地求生(吃鸡)游戏里,不同的枪支有不同的威力,更是可以搭配不同的配件,以提升枪支的性能. 每一把枪都有其威力及其可装备配件种类.每一个配件有其所属种类,可以为枪支提供威力的百分比加成.每 ...
大吉大利-今晚吃鸡-跑毒
题目描述现在有一款很火的游戏playerunknown's battlegrounds,人称"吃鸡",在里面经常面临跑毒(从安全区外跑进安全区内)的问题,在安全区外,人们会处于中 ...
大吉大利今晚吃鸡之跑毒篇
现在有一款很火的游戏playerunknown's battlegrounds,人称"吃鸡",在里面经常面临跑毒(从安全区外跑进安全区内)的问题,在安全区外,人们会处于中毒状态,每 ...
大吉大利今晚吃鸡——枪械篇
题目链接:点击打开链接这题难点带就是输入?我是因为用float精度不够才不过的,改成double就ok了 #include <cstdio> using namespace std;in ...
大吉大利今晚吃鸡——跑毒篇
题目链接:点击打开链接问题 1. 边界部分恰好擦着的时候到底是死是活.......写到程序里就是有没有=的问题. 首先题目里明确指出打包的时候擦着就是死. 然后后面有个带着你理解第二个样例的部分,这 ...
计算机内存不足吃鸡怎么办,Win10玩吃鸡游戏提示虚拟内存不足怎么办？
Win10玩吃鸡游戏提示"虚拟内存不足"怎么办? 最近吃鸡游戏非常火,大家都想大吉大利晚上吃鸡!因为吃鸡游戏对电脑配置的要求非常高,所以很多用户在玩吃鸡手游时会遇到各种问题,最常见 ...
吃鸡一直显示服务器,吃鸡一直显示服务器未连接 | 手游网游页游攻略大全
发布时间:2017-05-12 怎么吃鸡?相信很多新手玩家在经过几十甚至上百小时的游戏洗礼之后,心中对于吃鸡的疑惑也是越来越大,为什么有时候明明开局不顺却最后吃鸡,而有时候明明很肥却送了快递.所以今天 ...
为什么吃鸡显示连接不到服务器,为什么吃鸡进游戏显示连接不上 | 手游网游页游攻略大全...
发布时间:2015-09-28 原因:显卡太老,不支持DX9.0c+SM3.0 解决方法:换新显卡标签: 发布时间:2017-05-12 怎么吃鸡?相信很多新手玩家在经过几十甚至上百小时的游戏洗礼之 ...

【BZOJ5109】【CodePlus2017】大吉大利，晚上吃鸡！

【BZOJ5109】【CodePlus2017】大吉大利，晚上吃鸡！相关推荐

最新文章

热门文章