[JSOI2019]精准预测
这么明显的限制条件显然是\(\text{2-sat}\)
考虑按照时间拆点,\((0/1,x,t)\)表示\(x\)个人在时间\(t\)是生/死
有一些显然的连边
\[(0,x,t+1)->(0,x,t)\]
就是如果想在\(t+1\)时刻还是活着那必须在\(t\)时刻还活着
\[(1,x,t)->(1,x,t+1)\]
就是如果\(t\)时刻死了,那么\(t+1\)时刻也得是死的
对于第一种限制,显然是连\((1,x,t)->(1,y,t+1)\),别忘了连对称边,就是\(y\)在\(t+1\)时刻活着,那么\(x\)在\(t\)时刻可不能死,于是\((0,y,t+1)->(0,x,t)\)
对于第二种限制,连\((0,x,t)->(1,y,t)\),对称边连\((0,y,t)->(1,x,t)\)
发现按照上面的方式建图我们仅是点数就是\(nT\)级别,不能接受
发现每个人的生/死都连成了一条链,我们把上面每一其他边相连的点缩掉,这样图的点数就达到了\(n+m\)级别
简单分析我们发现这是一个\(\rm DAG\)
要对每个人求出有多少个人能和他一起活到最后
首先如果第\(i\)个人,\((0,i,T+1)\)能走到\((1,i,T+1)\),即\(i\)活着能推出\(i\)死了,于是\(i\)显然不能活着,自然也就不可能有点能和\(i\)一起活到最后。对于这样的点直接输出\(0\)即可。
否则,我们就看看如果这个点活着,那么能推出多少个点死了,即\((0,i,T+1)\)能到达多少种不同的\((1,y,T+1)\),这是经典的有向图传递闭包问题,只能做到\(\rm bitset\)复杂度。
那\(\rm bitset\)复杂度就可以了,\(O(\frac{m(n+m)}{w})\)大概就能过了。
但是这样\((0,i,T+1)\)可能没有办法走到那些必死点,我们把必死点做一个\(\rm bitset\),和\(i\)点的\(\rm bitset\)取一个或就可了。
代码
#include<bits/stdc++.h>
#define re register
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#pragma GCC optimize(3)
#pragma GCC optimize("-fcse-skip-blocks")
#define getchar() (S==TT&&(TT=(S=BB)+fread(BB,1,1<<15,stdin),S==TT)?EOF:*S++)
char BB[1<<18],*S=BB,*TT=BB;
inline int read() {
char c=getchar();int x=0;while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
const int maxn=4.6e5+1;
const int W=10000;
struct E{int v,nxt;}e[600005];
int head[maxn],vis[maxn];
std::bitset<10000> dp[maxn],g;
std::vector<int> id[2][50001],ti[50001];
int xx[100001],t[100001],yy[100001],op[100001];
int n,m,T,cnt,num,chk[50001],ans[50001];
inline void add(int x,int y) {e[++num].v=y;e[num].nxt=head[x];head[x]=num;}
inline int find(int x,int c) {
int l=0,r=ti[x].size()-1;
while(l<=r) {
int mid=l+r>>1;
if(ti[x][mid]==c) return mid;
if(ti[x][mid]<c) l=mid+1;else r=mid-1;
}
return 0;
}
void dfs(int x) {
if(vis[x]) return;vis[x]=1;
for(re int i=head[x];i;i=e[i].nxt)
dfs(e[i].v),dp[x]|=dp[e[i].v];
}
void calc(int l,int r) {
if(l>n) return;r=min(n,r);g.reset();
for(re int i=1;i<=cnt;i++) vis[i]=0,dp[i].reset();
for(re int i=l;i<=r;i++)
for(re int j=0;j<ti[i].size();++j)
dp[id[1][i][j]][i-l]=1;
for(re int i=1;i<=n;i++) dfs(id[0][i][ti[i].size()-1]);
for(re int i=l;i<=r;i++) if(dp[id[0][i][ti[i].size()-1]][i-l]) g[i-l]=1,chk[i]=1;
for(re int i=1;i<=n;i++) {
if(chk[i]) continue;
int x=id[0][i][ti[i].size()-1];;
ans[i]+=(dp[x]|g).count();
}
}
int main() {
T=read(),n=read(),m=read();
for(re int i=1;i<=m;i++) {
op[i]=read(),t[i]=read(),xx[i]=read(),yy[i]=read();
ti[xx[i]].push_back(t[i]);
ti[yy[i]].push_back(t[i]+(op[i]==1?0:1));
}
for(re int i=1;i<=n;i++) ti[i].push_back(T+1);
for(re int i=1;i<=n;i++) {
std::sort(ti[i].begin(),ti[i].end());
ti[i].erase(std::unique(ti[i].begin(),ti[i].end()),ti[i].end());
for(re int j=0;j<ti[i].size();++j)
id[0][i].push_back(++cnt),id[1][i].push_back(++cnt);
for(re int j=0;j<ti[i].size()-1;++j)
add(id[1][i][j],id[1][i][j+1]);
for(re int j=ti[i].size()-1;j;--j)
add(id[0][i][j],id[0][i][j-1]);
}
for(re int i=1;i<=m;i++) {
int x=find(xx[i],t[i]),y;
if(op[i]) {
y=find(yy[i],t[i]);
add(id[0][xx[i]][x],id[1][yy[i]][y]);
add(id[0][yy[i]][y],id[1][xx[i]][x]);
}
else {
y=find(yy[i],t[i]+1);
add(id[1][xx[i]][x],id[1][yy[i]][y]);
add(id[0][yy[i]][y],id[0][xx[i]][x]);
}
}
calc(W+1,2*W);calc(4*W+1,5*W);calc(1,W);calc(3*W+1,4*W);calc(2*W+1,3*W);
for(re int i=1;i<=n;i++,putchar(' ')) if(chk[i]) putchar('0');else printf("%d",n-1-ans[i]);
return 0;
}
原文地址:https://www.cnblogs.com/asuldb/p/11564842.html
- Python 异常处理完整指南
- 有趣的安全实验:利用多线程资源竞争技术上传shell
- NotSerializableException解决方法NotSerializableException
- 谁蹭了我的WiFi?浅谈家用无线路由器攻防
- 当输入流和输出流同时作用一个文件
- 揭秘:充电宝是如何盗取你的个人隐私的?
- 微软:暴力破解面前,增强密码复杂性基本没用
- 创建被图像填充的组件解释几处做法解释几点
- Android原生嵌入React Native
- React Native控件之Listview
- React Native控件只TextInput
- Qzone React Native改造
- android 特卖列表倒计时卡顿问题
- ios zxing扫码问题
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- springboot开发之配置自定义的错误界面和错误信息
- springboot开发之删除员工
- LeetCode | 1.两数之和
- springboot开发之配置嵌入式Servlet容器两种方式
- jvaa之初始化块
- java之匿名内部类
- springboot配置之获取配置文件中属性的第二种方法(@Value)不同于@ConfigurationProperties
- mybatis动态sql之内置参数_parameter和_databaseId
- Redis | Redis 字符串相关命令
- java之抽象类
- springboot配置之yaml
- properties和yaml配置文件
- java之模板方法设计模式
- java之接口
- PHP 代码混淆处理思路