Bzoj5109: [CodePlus 2017]大吉大利，晚上吃鸡！

最新推荐文章于 2019-09-18 16:15:00 发布
最新推荐文章于 2019-09-18 16:15:00 发布
阅读量103 收藏
点赞数
原文链接：http://www.cnblogs.com/cjoieryl/p/8691806.html
版权

题面

传送门

Sol

先正反两遍\(Dijsktra\)算出经过某个点的\(S\)到\(T\)的最短路条数\(F\)

满足条件一就是要满足\(F(A)+F(B)=F(T)\)

条件二

标算比较简单

直接\(bitset\)存储不能到达它的和它不能到的点

然后开\(map\)把所有相同的\(F(B)\)变成\(bitset\)

然后每次枚举\(A\)，就直接查表然后用\(bitset\)里的\(count\)就好了

注意如果\(ST\)不连通输出\(C_n^2\)

空间\(1GB\)就可以过

注意！！！

\(AC\)的同志不要以为真的\(AC\)了！！！
数据真的水
\(S\)到\(T\)的路径条数不会大于\(1\)
然后我自己乱打一个，样例都没过就\(AC\)了

标算代码

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
const int _(50005);
typedef long long ll;

IL int Input(){
    RG int x = 0, z = 1; RG char c = getchar();
    for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
    for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
    return x * z;
}

int n, m, cnt, S, T, first[_], d[_], vis[_];
struct Edge{
    int to, next, w;
} edge[_ << 1];
ll dis[2][_], f[2][_], F[_], ans;
struct Point{
    int u; ll d;

    IL int operator <(RG Point B) const{
        return d > B.d;
    }
};
priority_queue <Point> Q;
queue <int> Que;
bitset <50000> G[2][_];
map <ll, bitset <50000> > M;

IL void Add_Graph(RG int u, RG int v, RG int w){
    edge[cnt] = (Edge){v, first[u], w}, first[u] = cnt++;
}

IL void Dijkstra(RG int op){
    for(RG int i = 1; i <= n; ++i) dis[op][i] = (ll)1e18, vis[i] = 0;
    if(!op) f[op][S] = 1, dis[op][S] = 0, Q.push((Point){S, 0});
    else f[op][T] = 1, dis[op][T] = 0, Q.push((Point){T, 0});
    while(!Q.empty()){
        RG Point x = Q.top(); Q.pop();
        if(vis[x.u]) continue;
        vis[x.u] = 1;
        for(RG int e = first[x.u]; e != -1; e = edge[e].next){
            RG int v = edge[e].to, w = edge[e].w;
            if(x.d + w < dis[op][v]){
                dis[op][v] = x.d + w, f[op][v] = f[op][x.u];
                Q.push((Point){v, x.d + w});
            }
            else if(x.d + w == dis[op][v]) f[op][v] += f[op][x.u];
        }
    }
}

IL int Check(RG int x, RG ll w, RG int y, RG int op){
    return dis[op][x] + w + dis[!op][y] == dis[0][T];
}

IL void TopSort(RG int op){
    for(RG int i = 1; i <= n; ++i)
        for(RG int e = first[i]; e != -1; e = edge[e].next){
            RG int v = edge[e].to, w = edge[e].w;
            if(Check(i, w, v, op)) ++d[v];
        }
    for(RG int i = 1; i <= n; ++i) if(!d[i]) Que.push(i);
    while(!Que.empty()){
        RG int u = Que.front(); Que.pop();
        for(RG int e = first[u]; e != -1; e = edge[e].next){
            RG int v = edge[e].to, w = edge[e].w;
            if(!Check(u, w, v, op)) continue;
            G[op][v] &= G[op][u];
            if(!--d[v]) Que.push(v);
        }
    }
}

int main(RG int argc, RG char* argv[]){
    Fill(first, -1), n = Input(), m = Input(), S = Input(), T = Input();
    for(RG int i = 1; i <= m; ++i){
        RG int u = Input(), v = Input(), w = Input();
        Add_Graph(u, v, w), Add_Graph(v, u, w);
    }
    Dijkstra(0);
    if(!f[0][T]) return printf("%lld\n", 1LL * n * (n - 1) >> 1), 0;
    Dijkstra(1);
    for(RG int i = 1; i <= n; ++i){
        G[0][i].set(), G[1][i].set();
        G[0][i][0] = G[0][i][i] = G[1][i][0] = G[1][i][i] = 0;
    }
    TopSort(0), TopSort(1);
    for(RG int i = 1; i <= n; ++i)
        if(dis[0][i] + dis[1][i] == dis[0][T]) F[i] = f[0][i] * f[1][i];
    for(RG int i = 1; i <= n; ++i) M[F[i]].set(i);
    for(RG int i = 1; i <= n; ++i)
        ans += (M[F[T] - F[i]] & G[0][i] & G[1][i]).count();
    return printf("%lld\n", ans >> 1), 0;
}

然而Bzoj的空间变成了512MB

TAT

然后还可以这样做

首先可以枚举一条最短路上的点\(i\)

然后枚举不在这条最短路上的点再来统计答案

把最短路弄出来到一个数组中

然后不在这条路上的点的合法的\(i\)是这个数组的下标的区间\([l, r]\)

然后可以正反两边拓扑排序求出这个区间

为什么是一个区间

首先前提是这个点走到\(j\)，然后\(j\)到\(T\)是最短路

  • 如果\(j\)能走到这个点，那么\(j\)的所有前驱都能到这个点

  • 如果这个点能走到\(j\)，那么这个点也能走到\(j\)的所有后继

然后枚举最短路上的点，每次到\(i\)时，就把左端点是\(i\)的加进\(map\)

然后把右端点是\(i\)的从\(map\)中减掉

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
const int _(50005);
typedef long long ll;

IL int Input(){
    RG int x = 0, z = 1; RG char c = getchar();
    for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
    for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
    return x * z;
}

int n, m, cnt, S, T, first[_], d[_], vis[_], l[_], r[_], p[_], len, pre[_];
struct Edge{
    int to, next, w;
} edge[_ << 1];
ll dis[2][_], f[2][_], F[_], ans;
struct Point{
    int u; ll d;

    IL int operator <(RG Point B) const{
        return d > B.d;
    }
};
queue <int> Que;
map <ll, int> M;
priority_queue <Point> Q;
vector <int> al[_], ar[_];

IL void Add_Graph(RG int u, RG int v, RG int w){
    edge[cnt] = (Edge){v, first[u], w}, first[u] = cnt++;
}

IL void Dijkstra(RG int op){
    for(RG int i = 1; i <= n; ++i) dis[op][i] = (ll)1e18, vis[i] = 0;
    if(!op) f[op][S] = 1, dis[op][S] = 0, Q.push((Point){S, 0});
    else f[op][T] = 1, dis[op][T] = 0, Q.push((Point){T, 0});
    while(!Q.empty()){
        RG Point x = Q.top(); Q.pop();
        if(vis[x.u]) continue;
        vis[x.u] = 1;
        for(RG int e = first[x.u]; e != -1; e = edge[e].next){
            RG int v = edge[e].to, w = edge[e].w;
            if(x.d + w < dis[op][v]){
                dis[op][v] = x.d + w, f[op][v] = f[op][x.u];
                Q.push((Point){v, x.d + w});
                if(op) pre[v] = x.u;
            }
            else if(x.d + w == dis[op][v]) f[op][v] += f[op][x.u];
        }
    }
}

IL int Check(RG int x, RG ll w, RG int y, RG int op){
    return dis[op][x] + w + dis[!op][y] == dis[0][T];
}

IL void TopSort(RG int op){
    for(RG int i = 1; i <= n; ++i)
        for(RG int e = first[i]; e != -1; e = edge[e].next){
            RG int v = edge[e].to, w = edge[e].w;
            if(Check(i, w, v, op)) ++d[v];
        }
    for(RG int i = 1; i <= n; ++i) if(!d[i]) Que.push(i);
    while(!Que.empty()){
        RG int u = Que.front(); Que.pop();
        for(RG int e = first[u]; e != -1; e = edge[e].next){
            RG int v = edge[e].to, w = edge[e].w;
            if(!Check(u, w, v, op)) continue;
            op ? r[v] = min(r[v], r[u]) : l[v] = max(l[v], l[u]);
            if(!--d[v]) Que.push(v);
        }
    }
}

int main(RG int argc, RG char* argv[]){
    Fill(first, -1), n = Input(), m = Input(), S = Input(), T = Input();
    for(RG int i = 1; i <= m; ++i){
        RG int u = Input(), v = Input(), w = Input();
        Add_Graph(u, v, w), Add_Graph(v, u, w);
    }
    Dijkstra(0);
    if(!f[0][T]) return printf("%lld\n", 1LL * n * (n - 1) >> 1), 0;
    Dijkstra(1);
    for(RG int i = S; i; i = pre[i]) p[++len] = i, l[i] = len + 1, r[i] = len - 1;
    for(RG int i = 1; i <= n; ++i) if(!(l[i] + r[i])) l[i] = 1, r[i] = len;
    TopSort(0), TopSort(1);
    for(RG int i = 1; i <= n; ++i){
        if(dis[0][i] + dis[1][i] == dis[0][T]) F[i] = f[0][i] * f[1][i];
        if(l[i] > r[i]) continue;
        al[l[i]].push_back(i), ar[r[i]].push_back(i);
    }
    for(RG int i = 1; i <= len; ++i){
        for(RG int l = al[i].size(), j = 0; j < l; ++j) ++M[F[al[i][j]]];
        ans += M[F[T] - F[p[i]]];
        for(RG int l = ar[i].size(), j = 0; j < l; ++j) --M[F[ar[i][j]]];
    }
    return printf("%lld\n", ans), 0;
}

转载于:https://www.cnblogs.com/cjoieryl/p/8691806.html

weixin_34208185
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
BZOJ_problem_rank:BZOJ问题等级按交流比排序
05-25
#BZOJ Problem Rankrank.cpp 程序文件data.dat bzoj题库数据done.dat AC过的题，初始可以把所有A过的题粘进去，正常退出的话自动维护。black.dat 黑名单。选题时会跳过。错题、神题、没题面、不想做等等。//Thank ...
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
BZOJ5109：[CodePlus 2017]大吉大利，晚上吃鸡！ （最短路+Hash表+二进制压位）
KsCla
12-04 2337
题目传送门：http://www.lydsy.com/JudgeOnline/problem.php?id=5109题目分析：过了挺久终于把这个坑填了。一开始以为是一道很难的题，后来发现也不难想。由于懒得打题解，直接引用出题人的题解好了（主要是来贴代码）QAQ： 虽然题目中给定的是无向图，但是实际上我们可以先从 SS 出发求一遍最短路，然后问题变成了：“在有向无环图上，求有多少个满足条件的点对
BZOJ 5109 [CodePlus 2017]大吉大利，晚上吃鸡！
无尽
08-19 303
最短路+bitset+DP 个人觉得这题的思路非常高妙。 首先肯定是要建出最短路DAG，这个图上任意一条路径都对应一条原图的最短路。 如果一个点a在S到T的必经之路上，那就会有 S到a的方案数 * a到T的方案数 = S到T的方案数 这个东西显然是充要的，这是一个巧妙的转化。套用这个想法，这题要求选出两个点，那就只需 S到a的方案数 * a到T的方案数 + S到b的方案数 * b到T的...
[BZOJ5109/CodePlus2017]大吉大利，晚上吃鸡！
a5163273的博客
10-05 208
Description 最近《绝地求生：大逃杀》风靡全球，皮皮和毛毛也迷上了这款游戏，他们经常组队玩这款游戏。在游戏中，皮皮和毛毛最喜欢做的事情就是堵桥，每每有一个好时机都能收到不少的快递。当然，有些时候并不能堵桥，皮皮和毛毛会选择在其他的必经之路上蹲点。K博士作为一个老年人，外加有心脏病，自然是不能玩这款游戏的，但是这并不能妨碍他对这款游戏进行一些理论分析，比如最近他就对皮皮和毛毛的战...
「CodePlus 2017 11 月赛」大吉大利，晚上吃鸡！（dij+bitset）
weixin_30469895的博客
12-26 129
　　从S出发跑dij，从T出发跑dij，顺便最短路计数。 　　令$F(x)$为$S$到$T$最短路经过$x$的方案数，显然这个是可以用$S$到$x$的方案数乘$T$到$x$的方案数来得到。 　　然后第一个条件就变成了满足$F(A)+F(B)=F(T)$，这个只要用map存一下点的状态，每次查$F(T)-F(A)$就可以得到$B$的状态了。 　　第二个条件实际上就是$A$无法到达$B$，怎么判...
#4【BZOJ5109】[CodePlus 2017]大吉大利，晚上吃鸡！（未完成）
weixin_30267785的博客
04-20 117
。。 bitset优化floyd 转载于:https://www.cnblogs.com/yinwuxiao/p/8893879.html
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
BZOJ3230相似子串测试数据
08-19
标题"BZOJ3230相似子串测试数据"指的是一个与在线判题平台BZOJ（BestCoder ZOJ，中国大学算法竞赛在线评测系统）相关的编程问题，具体问题编号为3230，涉及的主要知识点是"相似子串"的搜索和处理。在编程竞赛中，...
BZOJ上面的题目以及数据
10-13
【BZOJ】全称是“北京邮电大学在线评测系统”，它是一个为编程爱好者和参赛者提供在线编程练习和比赛的平台。该平台涵盖了多种算法和编程问题，旨在提升用户的编程技能和解决复杂问题的能力。BZOJ支持C、C++、Java等...
BZOJ 5109 大吉大利，晚上吃鸡 Code+
Rayment的博客
04-13 412
Problem BZOJ Solution 这道题的官方正解好像是bitset？ 据说数据很水的说，随机数据要么不连通要么就只有一条最短路 dp一下即可，记录每个点被几条最短路覆盖。 Code #include &lt;cstring&gt; #include &lt;cstdio&gt; #include &lt;queue&gt; using namespace std;...
BZOJ 5109: [CodePlus 2017]大吉大利，晚上吃鸡！ 最短路 拓扑 Dp
dilu0653的博客
09-18 309
title BZOJ 5109 LUOGU 4061 题意： 一张 \(n\) 个点 \(m\) 条无向边的图，节点编号为 \(1\) 到 \(n\) ，每条边具有一个正整数的长度。 假定大魔王都会从 \(S\) 点出发到达 \(T\) 点（ \(S\) 和 \(T\) 已知），并且只会走最短路，皮皮和毛毛会在 \(A\) 点和 \(B\) 点埋伏大魔王。 为了保证一定能埋伏到大魔王，...
[最短路][DP][传递闭包] BZOJ 5109 && LOJ #6252. 「CodePlus 2017 11 月赛」大吉大利，晚上吃鸡！
Vectorxj的博客
12-27 1129
SolutionSolution记fuf_u表示从SS到TT经过点uu的方案数。 这可以很简单的DP出来。 那么点对(A,B)(A,B)合法当且仅当 fA+fB=fTf_A+f_B=f_T AA，BB无法互相到达。 第一个用个表存，第二个按拓扑序做传递闭包就好了。 原来图会不连通 时间复杂度O(nlogn+nmw)O(n\log n+{nm\over w})#include <bits/st
【BZOJ5109】[CodePlus 2017]大吉大利，晚上吃鸡！ 最短路+拓扑排序+DP
aodanchui1057的博客
11-27 342
【BZOJ5109】[CodePlus 2017]大吉大利，晚上吃鸡！ Description 最近《绝地求生：大逃杀》风靡全球，皮皮和毛毛也迷上了这款游戏，他们经常组队玩这款游戏。在游戏中，皮皮和毛毛最喜欢做的事情就是堵桥，每每有一个好时机都能收到不少的快递。当然，有些时候并不能堵桥，皮皮和毛毛会选择在其他的必经之路上蹲点。K博士作为一个老年人，外加有心脏病，自然是不能玩这款游...
[传递闭包 BITSET] 「CodePlus 2017 11 月赛」大吉大利，晚上吃鸡！
CE玩家
11-27 1724
建出最短路图—（以下复制自官方题解） 定义 F(X)=F(X) = 从 SS 到 XX 的方案数 ×\times 从 XX 到 TT 的方案数 = 从 SS 经过 XX 到达 TT 的方案数，所以满足条件的点对 A,BA,B 为:F(A)+F(B)=F(T)F(A) + F(B) = F(T) AA 和 BB 不能相互到达对于条件 11 ，我们可以使用数据结构进行优化（使用std::map即可
[BZOJ5109]大吉大利，晚上吃鸡！
weixin_33910460的博客
08-22 136
[BZOJ5109]大吉大利，晚上吃鸡！ 题目大意： 一张\(n(n\le5\times10^4)\)个点\(m(m\le5\times10^4)\)条边的无向图，节点编号为\(1\)到\(n\)，边权为正整数。给定\(S\)和\(T\)，显然从\(S\)到\(T\)的最短路有一种或多种方案。 选择\(A,B\)两个点，约定\(A\)点和\(B\)点必须满足： 所有可能路径中，必定会经过\(A\...
bzoj5109(图论好题)
weixin_34297704的博客
01-23 283
我的参考题解：https://www.cnblogs.com/ccz181078/p/7907022.html； 不过我感觉题解的压位有问题，（1<<x）还不炸上天。不过这题数据水，好像怎么写都能对，这里放上我认为正确的写法 //正经解法是计算一个F(x)表示经过x的最短路条数，然后找F(A)+F(B)=F(T),用map来做。 //而且F(x)可能很大，还要取模，一个不够的...
[BZOJ5109][LOJ #6252][P4061][CodePlus 2017 11月赛]大吉大利，今晚吃鸡！(最短路+拓扑排序+传递闭包+map+bitset(hash+压位))...
weixin_30389003的博客
02-23 142
5109: [CodePlus 2017]大吉大利，晚上吃鸡！ Time Limit: 30 SecMemory Limit: 1024 MBSubmit: 107Solved: 57[Submit][Status][Discuss] Description 最近《绝地求生：大逃杀》风靡全球，皮皮和毛毛也迷上了这款游戏，他们经常组队玩这款游戏。在游戏中，皮皮 和毛毛最...
BZOJ2839:集合计数(容斥,组合数学)
最新发布
06-08
这个问题是一个比较经典的容斥原理的应用题目。我们可以考虑使用容斥原理来求解这个问题。具体来说，我们可以考虑求出所有集合的个数，然后减去至少包含一个特定元素的集合的个数，再加上至少包含两个特定元素的集合...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值