博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
bzoj3931: [CQOI2015]网络吞吐量
阅读量:5217 次
发布时间:2019-06-14

本文共 3251 字,大约阅读时间需要 10 分钟。

将最短路图找出来,跑maxflow即可。有注意到数据范围。然后输出的时候%dWA了三次QAQ。。。

#include
#include
#include
#include
#include
using namespace std;#define rep(i,n) for(int i=1;i<=n;i++)#define ll long long#define clr(x,c) memset(x,c,sizeof(x))#define qwq(x) for(Edge *o=h[x];o;o=o->next)#define QWQ(x) for(edge *o=head[x];o;o=o->next)#define REP(i,s,t) for(int i=s;i<=t;i++)int read(){ int x=0;char c=getchar();bool f=true; while(!isdigit(c)) { if(c=='-') f=false; c=getchar(); } while(isdigit(c)) x=x*10+c-'0',c=getchar(); return f?x:-x;} const int nmax=1005;const int maxn=400005;const ll inf=(ll)1<<60;struct edge{ int to;ll cap;edge *next,*rev;};edge edges[maxn],*pt=edges,*head[nmax],*p[nmax],*cur[nmax];void add(int u,int v,ll w){ pt->to=v;pt->cap=w;pt->next=head[u];head[u]=pt++;}void adde(int u,int v,ll w){ add(u,v,w);add(v,u,0);head[u]->rev=head[v];head[v]->rev=head[u];}struct Edge{ int to;ll dist;Edge *next;};Edge Edges[maxn],*e=Edges,*h[nmax];void addedge(int u,int v,ll w){ e->to=v;e->dist=w;e->next=h[u];h[u]=e++;}int cnt[nmax],H[nmax];bool inq[nmax];ll d[nmax];void spfa(int s,int t,int n){ rep(i,n) d[i]=inf; clr(inq,0);inq[s]=1;d[s]=0; queue
q;q.push(s); while(!q.empty()){ int x=q.front();q.pop();inq[x]=0; qwq(x) if(d[o->to]>d[x]+o->dist){ d[o->to]=d[x]+o->dist; if(!inq[o->to]) q.push(o->to),inq[o->to]=1; } }}ll maxflow(int s,int t,int n){ clr(cnt,0);cnt[0]=n;clr(H,0); ll flow=0,a=inf;int x=s;edge *e; while(H[s]
next) if(e->cap>0&&H[x]==H[e->to]+1) break; if(e){ p[e->to]=cur[x]=e;a=min(a,e->cap);x=e->to; if(x==t){ while(x!=s) p[x]->cap-=a,p[x]->rev->cap+=a,x=p[x]->rev->to; flow+=a,a=inf; } }else{ if(!--cnt[H[x]]) break; H[x]=n; for(e=head[x];e;e=e->next) if(e->cap>0&&H[x]>H[e->to]+1) H[x]=H[e->to]+1,cur[x]=e; cnt[H[x]]++; if(x!=s) x=p[x]->rev->to; } } return flow;}int main(){ int n=read(),m=read(),u,v;ll w; rep(i,m) u=read(),v=read(),w=read(),addedge(u,v,w),addedge(v,u,w); /*rep(i,n){ qwq(i) printf("%d ",o->to);printf("\n"); }*/ spfa(1,n,n); //rep(i,n) printf("%d ",d[i]); rep(i,n) qwq(i) if(d[o->to]==d[i]+o->dist) if(i==1||i==n) adde(i,o->to,inf); else adde(i+n,o->to,inf); /*rep(i,n){ QWQ(i) printf("%d ",o->to);printf("\n"); }*/ rep(i,n) { w=read(); if(i!=1&&i!=n) adde(i,i+n,w); } printf("%lld\n",maxflow(1,n,n+n-2)); return 0;}

 

3931: [CQOI2015]网络吞吐量

Time Limit: 10 Sec  Memory Limit: 512 MB
Submit: 1328  Solved: 549
[][][]

Description

 路由是指通过计算机网络把信息从源地址传输到目的地址的活动,也是计算机网络设计中的重点和难点。网络中实现路由转发的硬件设备称为路由器。为了使数据包最快的到达目的地,路由器需要选择最优的路径转发数据包。例如在常用的路由算法OSPF(开放式最短路径优先)中,路由器会使用经典的Dijkstra算法计算最短路径,然后尽量沿最短路径转发数据包。现在,若已知一个计算机网络中各路由器间的连接情况,以及各个路由器的最大吞吐量(即每秒能转发的数据包数量),假设所有数据包一定沿最短路径转发,试计算从路由器1到路由器n的网络的最大吞吐量。计算中忽略转发及传输的时间开销,不考虑链路的带宽限制,即认为数据包可以瞬间通过网络。路由器1到路由器n作为起点和终点,自身的吞吐量不用考虑,网络上也不存在将1和n直接相连的链路。

 

Input

输入文件第一行包含两个空格分开的正整数n和m,分别表示路由器数量和链路的数量。网络中的路由器使用1到n编号。接下来m行,每行包含三个空格分开的正整数a、b和d,表示从路由器a到路由器b存在一条距离为d的双向链路。 接下来n行,每行包含一个正整数c,分别给出每一个路由器的吞吐量。

 

Output

输出一个整数,为题目所求吞吐量。

 

Sample Input

7 10
1 2 2
1 5 2
2 4 1
2 3 3
3 7 1
4 5 4
4 3 1
4 6 1
5 6 2
6 7 1
1
100
20
50
20
60
1

Sample Output

70

HINT

 

 对于100%的数据,n≤500,m≤100000,d,c≤10^9

 

Source

 
[ ][ ][ ]

转载于:https://www.cnblogs.com/fighting-to-the-end/p/5657138.html

你可能感兴趣的文章
关于MFC中窗口的销毁
查看>>
bzoj 2007: [Noi2010]海拔【最小割+dijskstra】
查看>>
BZOJ 1001--[BeiJing2006]狼抓兔子(最短路&对偶图)
查看>>
C# Dynamic通用反序列化Json类型并遍历属性比较
查看>>
128 Longest Consecutive Sequence 一个无序整数数组中找到最长连续序列
查看>>
定制jackson的自定义序列化(null值的处理)
查看>>
auth模块
查看>>
Java使用FileReader(file)、readLine()读取文件,以行为单位,一次读一行,一直读到null时结束,每读一行都显示行号。...
查看>>
Elipse安装Spring Tool Suite
查看>>
javascript keycode大全
查看>>
前台freemark获取后台的值
查看>>
log4j.properties的作用
查看>>
滴滴快车历史奖励政策:含工作日和周末的高峰奖励、翻倍奖励【历史政策】...
查看>>
文件操作类2
查看>>
游戏偶感
查看>>
Leetcode: Unique Binary Search Trees II
查看>>
转载------------Python多线程学习
查看>>
phpcurl类
查看>>
Hadoop伪分布式搭建
查看>>
C++ FFLIB 之FFDB: 使用 Mysql&Sqlite 实现CRUD
查看>>