「CodePlus 2017 11 月赛」大吉大利,晚上吃鸡!
n<=50000,m<=50000的图,给s和t,问有多少点对$(a,b)$满足
嗯。
不会。
首先最短路DAG造出来,然后两个条件转述一下:条件一,$N_a$表示从s到t经过a的路径,$N_a+N_b=N_t$;条件二,在最短路DAG上A不能到B,B不能到A。
条件一就迪杰斯特拉的时候算一下N,注意不在最短路DAG上的点$N_i=0$;然后对每个$N_t-N_b$的值存一个bitset,用以表示值为这么多的点的状态,枚举a查多少$N_t-N_b=N_a$即可。
条件二就正反拓扑序跑一下,然后传递闭包算出来即可知道最短路图上哪些点能到a和a能到哪些点,把这些点设为不可达点,取个交集即可算出每个a能和哪些b在条件二下配对。
然后就没了。
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<stdlib.h>
#include<bitset>
#include<queue>
#include<math.h>
//#include<time.h>
//#include<iostream>
using namespace std; int n,m,s,t;
#define maxn 50011
#define maxm 100011
#define LL long long struct Edge{int to,next,v;};
struct qnode
{
int id; LL v;
bool operator < (const qnode &b) const {return v<b.v;}
bool operator > (const qnode &b) const {return v>b.v;}
};
struct Graph
{
Edge edge[maxm<<]; int first[maxn],le;
Graph() {memset(first,,sizeof(first)); le=;}
void in(int x,int y,int v) {Edge &e=edge[le]; e.to=y; e.v=v; e.next=first[x]; first[x]=le++;}
void insert(int x,int y,int v) {in(x,y,v); in(y,x,v);}
priority_queue<qnode,vector<qnode>,greater<qnode> > q;
void dijkstra(int s,LL *dis,LL *f)
{
for (int i=;i<=n;i++) dis[i]=1e18,f[i]=;
dis[s]=; f[s]=;
q.push((qnode){s,});
while (!q.empty())
{
const int now=q.top().id; const LL d=q.top().v; q.pop();
if (d!=dis[now]) continue;
for (int i=first[now];i;i=edge[i].next)
{
const Edge &e=edge[i];
if (dis[e.to]>dis[now]+e.v)
{
dis[e.to]=dis[now]+e.v;
f[e.to]=f[now];
q.push((qnode){e.to,dis[e.to]});
}
else if (dis[e.to]==dis[now]+e.v) f[e.to]+=f[now];
}
}
}
}g; LL dis[][maxn],f[][maxn],val[maxn];
bitset<maxn> where[maxn],can[][maxn];
int indo[maxn],head,tail,que[maxn];
bool check(int x,int y,int v,int ty) {return dis[ty][x]+v+dis[ty^][y]==dis[][t];}
void toposort(int ty)
{
memset(indo,,sizeof(indo));
for (int i=;i<=n;i++)
for (int j=g.first[i];j;j=g.edge[j].next)
{
const Edge &e=g.edge[j];
if (check(i,e.to,e.v,ty)) indo[e.to]++;
}
head=tail=;
for (int i=;i<=n;i++) if (indo[i]==) que[tail++]=i;
while (head!=tail)
{
const int now=que[head++];
for (int i=g.first[now];i;i=g.edge[i].next)
{
const Edge &e=g.edge[i];
if (!check(now,e.to,e.v,ty)) continue;
indo[e.to]--; if (indo[e.to]==) que[tail++]=e.to;
can[ty][e.to]&=can[ty][now];
}
}
} qnode list[maxn];
int main()
{
scanf("%d%d%d%d",&n,&m,&s,&t);
for (int i=,x,y,v;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&v);
g.insert(x,y,v);
}
g.dijkstra(s,dis[],f[]); g.dijkstra(t,dis[],f[]); if (dis[][t]==1e18) {printf("%lld\n",1ll*n*(n-)/); return ;} for (int i=;i<=n;i++) can[][i].set(),can[][i].set(),can[][i][i]=can[][i][i]=;
toposort(); toposort(); for (int i=;i<=n;i++) if (dis[][i]+dis[][i]!=dis[][t]) val[i]=; else val[i]=f[][i]*f[][i];
for (int i=;i<=n;i++) list[i]=((qnode){i,val[i]});
sort(list+,list++n);
list[n+].v=1e18;
for (int i=,j=;i<=n+;i++) if (list[i].v!=list[i-].v)
{
int now=j;
for (;j<i;j++) where[now][list[j].id]=;
} LL ans=;
for (int i=;i<=n;i++)
{
int L=,R=n+; LL tmp=f[][t]-val[i];
while (L<R)
{
int mid=(L+R)>>;
if (list[mid].v>=tmp) R=mid;
else L=mid+;
}
if (L!=n+) ans+=(where[L]&can[][i]&can[][i]).count();
}
printf("%lld\n",ans/);
return ;
}
「CodePlus 2017 11 月赛」大吉大利,晚上吃鸡!的更多相关文章
- 「CodePlus 2017 11 月赛」大吉大利,晚上吃鸡!(dij+bitset)
从S出发跑dij,从T出发跑dij,顺便最短路计数. 令$F(x)$为$S$到$T$最短路经过$x$的方案数,显然这个是可以用$S$到$x$的方案数乘$T$到$x$的方案数来得到. 然后第一个条件就变 ...
- LOJ6252. 「CodePlus 2017 11 月赛」大吉大利,晚上吃鸡! 最短路+bitset
题目传送门 https://loj.ac/problem/6252 https://lydsy.com/JudgeOnline/problem.php?id=5109 题解 首先跑最短路,只保留 \( ...
- loj #6250. 「CodePlus 2017 11 月赛」找爸爸
#6250. 「CodePlus 2017 11 月赛」找爸爸 题目描述 小 A 最近一直在找自己的爸爸,用什么办法呢,就是 DNA 比对. 小 A 有一套自己的 DNA 序列比较方法,其最终目标是最 ...
- [LOJ 6249]「CodePlus 2017 11 月赛」汀博尔
Description 有 n 棵树,初始时每棵树的高度为 H_i,第 i 棵树每月都会长高 A_i.现在有个木料长度总量为 S 的订单,客户要求每块木料的长度不能小于 L,而且木料必须是整棵树(即不 ...
- [LOJ 6248]「CodePlus 2017 11 月赛」晨跑
Description “无体育,不清华”.“每天锻炼一小时,健康工作五十年,幸福生活一辈子” 在清华,体育运动绝对是同学们生活中不可或缺的一部分.为了响应学校的号召,模范好学生王队长决定坚持晨跑.不 ...
- 「CodePlus 2017 11 月赛」Yazid 的新生舞会(树状数组/线段树)
学习了新姿势..(一直看不懂大爷的代码卡了好久T T 首先数字范围那么小可以考虑枚举众数来计算答案,设当前枚举到$x$,$s_i$为前$i$个数中$x$的出现次数,则满足$2*s_r-r > 2 ...
- 「CodePlus 2017 11 月赛」可做题
这种题先二进制拆位,显然改的位置只有每一段确定的数的开头和结尾,只需要对于每一个可决策位置都尝试一下填1和0,然后取min即可. #include<iostream> #include&l ...
- 「CodePlus 2017 11 月赛」Yazid 的新生舞会
n<=500000的数字,问有多少个区间的众数出现次数严格大于区间长度的一半. 这么说来一个区间就一个众数了,所以第一反应是枚举数字,对下标进行处理.然后没有第二反应.很好. 在枚举一个数字的时 ...
- [LOJ#6259]「CodePlus 2017 12 月赛」白金元首与独舞
[LOJ#6259]「CodePlus 2017 12 月赛」白金元首与独舞 试题描述 到河北省 见斯大林 / 在月光下 你的背影 / 让我们一起跳舞吧 うそだよ~ 河北省怎么可能有 Stalin. ...
随机推荐
- Linux终端下对话
1.首先查看当前账户 [hadoop@weekend01 ~]$ who hadoop tty1 2016-11-14 09:31 (:0) hadoop pts/0 ...
- 函数的返回值return
'''1.什么是返回值 返回值是一个函数的处理结果 2.为什么要有返回值 如果我们需要在程序中拿到函数的处理结果做进一步的处理,则需要函数必须有返回值 3.函数的返回值的应用 函数的返回值用retur ...
- [SPOJ1812]Longest Common Substring II 后缀自动机 多个串的最长公共子串
题目链接:http://www.spoj.com/problems/LCS2/ 其实两个串的LCS会了,多个串的LCS也就差不多了. 我们先用一个串建立后缀自动机,然后其它的串在上面跑.跑的时候算出每 ...
- javaEE web 系统安装时自定义初始化
通常JavaWeb项目在第一次启动时我们需要做一些初始化工作,比如:初始化一个管理员的登录账户和密码,配置缓存.定时任务等,这些操作可以通过手工修改数据库完成,但是容易出错且繁琐,而且也很麻烦.如果这 ...
- vue HTTP请求(针对vue-resource)
//初始化页面需要做什么事情 //点击后需要做什么事情 //鼠标.键盘.冒泡.默认行为等事件 //前端调用接口就是按需展示数据//created和mounted里面都可以做数据处理,唯一不同的是cre ...
- vue1.0生命周期
<!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...
- [转载]iTOP-4418开发板Ubuntu系统烧写方法分享
本文转自迅为论坛:http://topeetboard.com 开发平台:iTOP-4418开发板系统:Ubuntu 1. TF卡读写速度测试烧写 Ubuntu 对于 TF 卡的要求比较高,很多老旧的 ...
- 理解 call, apply 的用法
callcall() 方法使用一个指定的 this 值和单独给出的一个或多个参数来调用一个函数. function list() { return Array.prototype.slice.call ...
- jQuery 开始动画,停止动画
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- 网络基础编程_5.4聊天室-IOCP服务器
聊天室-IOCP服务器 main 创建完成端口内核对象(CreateIoCompletionPort) 获取核心数并创建线程(GetSystemInfo + CreateThread) 创建套接字并绑 ...