【Luogu】P2445动物园(最大流)
2024-08-24 02:53:38
题目本身还是比较水的吧……容易发现是最大流套上dinic跑一遍就好了,并不会超时。
比较不偷税的一点是关于某动物的所有目击报告都符合才能连边……qwqqwqqwq
#include<cstdio>
#include<algorithm>
#include<cctype>
#include<cstring>
#include<cstdlib>
#include<queue>
#define maxn 1000
using namespace std;
inline long long read(){
long long num=,f=;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') f=-;
ch=getchar();
}
while(isdigit(ch)){
num=num*+ch-'';
ch=getchar();
}
return num*f;
} int u[]={,,,,-};
int w[]={,,,-,}; struct Edge{
int next,to,val;
}edge[maxn*];
int head[maxn],num;
inline void addedge(int from,int to,int val){
edge[++num]=(Edge){head[from],to,val};
head[from]=num;
}
inline void add(int from,int to,int val){
addedge(from,to,val);
addedge(to,from,);
} inline int count(int i){ return i&?i+:i-; } int pre[maxn];
int Start,End;
bool vis[maxn];
int dfn[maxn];
int list[maxn]; bool bfs(){
memset(vis,,sizeof(vis)); dfn[Start]=; vis[Start]=;
queue<int>q; q.push(Start);
while(!q.empty()){
int from=q.front(); q.pop();
for(int i=head[from];i;i=edge[i].next){
int to=edge[i].to;
if(vis[to]||edge[i].val<=) continue;
vis[to]=; dfn[to]=dfn[from]+;
q.push(to);
}
}
return vis[End];
} int dfs(int x,int val){
if(val==||x==End) return val;
vis[x]=; int flow=;
for(int &i=list[x];i;i=edge[i].next){
int to=edge[i].to;
if(vis[to]||dfn[to]!=dfn[x]+||edge[i].val<=) continue;
int now=dfs(to,min(val,edge[i].val));
if(x!=End) pre[to]=x; edge[i].val-=now; edge[count(i)].val+=now; val-=now; flow+=now;
if(val<=) break;
}
if(val!=flow) dfn[x]=-;
return flow;
} inline int maxflow(){
int ans=;
while(bfs()){
memset(vis,,sizeof(vis));
for(int i=Start;i<=End;++i) list[i]=head[i];
int now=dfs(Start,0x7fffffff);
if(now==) break;
ans+=now;
}
return ans;
} int s[maxn][maxn];
char c[maxn];
bool ext[maxn][maxn];
int dis[][][];
int v[maxn];
bool exa[maxn];
int f[maxn][maxn];
int sum[maxn]; struct Node{
int x,y;
}que[maxn]; int main(){
memset(dis,/,sizeof(dis)); int inf=dis[][][];
int n=read();
for(int i=;i<=n;++i){
scanf("%s",c+);
for(int j=;j<=n;++j)
if(c[j]=='*') s[i][j]=;
else s[i][j]=;
}
int p=read(); End=p*+;
for(int i=;i<=p;++i) que[i]=(Node){read(),read()};
for(int i=;i<=p;++i){
add(Start,i,); add(i+p,End,);
queue<Node>q; q.push(que[i]); dis[i][que[i].x][que[i].y]=;
while(!q.empty()){
Node from=q.front(); q.pop();
//printf("%d %d %d<<<<<<<<\n",i,from.x,from.y);
for(int j=;j<;++j){
int nx=from.x+u[j];
int ny=from.y+w[j];
//printf("%d %d %d %d><>\n",nx,ny,s[nx][ny],dis[i][from.x][from.y]);
if(nx<||nx>n||ny<||ny>n||s[nx][ny]||dis[i][nx][ny]!=inf) continue;
dis[i][nx][ny]=dis[i][from.x][from.y]+;
q.push((Node){nx,ny});
}
}
}
for(int i=;i<=p;++i) v[i]=read();
int r=read();
for(int i=;i<=r;++i){
int t=read(),x=read(),y=read(),d=read(); sum[d]++;
for(int j=;j<=p;++j){
//printf("%d %d>>>\n",j,dis[j][x][y]);
int ret=v[d]*t;
if(ret>=dis[j][x][y]) f[j][d]++;
}
}
for(int i=;i<=p;++i)
for(int j=;j<=p;++j)
if(sum[i]==f[j][i]) add(j,i+p,);
maxflow();
for(int i=+p;i<End;++i){
for(int j=head[i];j;j=edge[j].next){
int to=edge[j].to;
if(to==End||edge[j].val==) continue;
printf("%d %d %d\n",i-p,que[to].x,que[to].y);
break;
}
}
return ;
}
最新文章
- 设置tomcat远程debug
- Windows批处理:自动检查网络连通性
- eclipse配置gradle
- OSX10.11 删除系统自带的软件
- foundation框架之反射机制
- hive 操作(转)
- Mac双系统切换
- 社区发现(Community Detection)算法 [转]
- 如何重载ComboBox 使其下拉按钮(带下箭头的)和下拉列表的垂直滚动条的宽度改变?(自绘ComboBox) [转]
- 实现一个servlet的helloworld小程序(不适用Eclipse)
- web端及时通讯原理
- 关于const限定符
- [转]JAVA的动态代理机制及Spring的实现方式
- python3 第四章 - 输入与输出
- react children
- 20165223 《JAVA程序设计》第五周学习总结
- tcpdump抓包常用参数和用法
- hadoop动态添加删除节点datanode及恢复
- svg 配合cesium使用
- UVM:8.4.3 用factory 机制创建实例的接口
热门文章
- 如何在Sierra运行 Specials K 的patch
- 使用com.sun.imageio.plugins.png.PNGMetadata读取图片的元数据
- MyEclipse Update Progress Error解决方法
- mybatis(一):思维导图
- C#数组排序方法
- Unicode与ASCiI之间有什么区别?java当中的转义字符 Character类的使用 String类的使用
- 初学redis,redis基本数据类型
- 201621123080 《Java程序设计》 第7周学习总结
- ubuntu 设置定时任务
- SQL语句小练习