问题描述
α−76星人入侵地球，尽管英勇的地球卫士们奋勇杀敌，但战火依然烧到了地球本球上，为了维护人类的利益，联合国决定建立m个避难所，第i个避难所能容纳ci个人，坐标为(pi,qi)。现在地球上还残余着n个国家，第i个国家有bi个人，坐标为(xi,yi)，现在联合国给出了一个初步的疏散方案，即一个n*m的矩阵a，其中ai,j表示第i个国家去第j个避难所的人数（\(b_i=\sum_{1\le j\le m}a_{i,j}\)且\(\sum_{1\le i\le n}a_{i,j}\le c_j\)）
不过，运送这些人是要花钱的，运送一个人从坐标(x1,y1)到(x2,y2)要花费|x1−x2|+|y2−y1|+1美元，现在要求你找到一种避难方案使得总代价最小，输出这个总代价。
输入格式
第一行两个整数n,m表示国家数和避难所数
接下来n行每行三个整数xi,yi表示每个国家的坐标
接下来m行每行三个整数pi,qi,ci表示每个避难所的坐标和能容纳的人数
接下来n行，每行m个整数表示联合国给出的避难方案，a[i][j]表示第i个国家去第j个避难所避难的人数。
输出格式
一行一个整数表示最小的总代价。
样例输入
输入样例1:
1 1
615 -707
-904 -91 24
14
输入样例2:
2 5
750 593
734 364
-244 476 39
946 287 46
743 -975 2
309 265 20
941 822 42
5 4 2 20 32
34 15 0 0 10
输入样例3:
6 1
615 -985
173 932
-133 668
790 -755
-428 524
-308 904
-559 -81 31
28
3
0
0
0
0
样例输出
输出样例1:
14520
输出样例2:
87700
输出样例3:
63450
数据范围
30%的数据： 1<=n,m<=10
60%的数据： 1<=n,m<=50
100%的数据： 1<=n,m<=100,坐标在[-1000,1000]，1<=c[i],a[i][j]<=50
解析
这……由于题目的问题，只要按照联合国方案计算代价输出即可通过。
没错，样例都过不了，但他就是能过。
代码本来是写费用流的，后来发现了，也没有改代码。
代码
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#define N 1000002
#define M 1000002
using namespace std;
const int inf=1<<30;
int head[N],ver[M*2],nxt[M*2],cap[M*2],l;
int n,m,s,t,i,dis[N],ans;
int read()
{
char c=getchar();
int w=0;
while(c<'0'||c>'9') c=getchar();
while(c<='9'&&c>='0'){
w=w*10+c-'0';
c=getchar();
}
return w;
}
void insert(int x,int y,int z)
{
ver[l]=y;
cap[l]=z;
nxt[l]=head[x];
head[x]=l;
l++;
ver[l]=x;
nxt[l]=head[y];
head[y]=l;
l++;
}
bool bfs()
{
queue<int> q;
memset(dis,-1,sizeof(dis));
q.push(s);
dis[s]=0;
while(!q.empty()){
int x=q.front();
q.pop();
for(int i=head[x];i!=-1;i=nxt[i]){
int y=ver[i];
if(dis[y]==-1&&cap[i]>0){
dis[y]=dis[x]+1;
q.push(y);
}
}
}
return (dis[t]!=-1);
}
int dfs(int x,int flow)
{
if(x==t||flow==0) return flow;
int ans=0;
for(int i=head[x];i!=-1;i=nxt[i]){
int y=ver[i];
if(dis[y]==dis[x]+1&&cap[i]>0){
int a=dfs(y,min(flow,cap[i]));
flow-=a;
ans+=a;
cap[i]-=a;
cap[i^1]+=a;
}
if(flow==0) break;
}
if(flow) dis[x]=-1;
return ans;
}
int main()
{
memset(head,-1,sizeof(head));
n=read();m=read();
t=n+2*m+1;
for(i=2;i<n;i++){
int va=read();
insert(s,i,va);
ans+=va;
}
for(i=2;i<n;i++){
int vb=read();
insert(i,t,vb);
ans+=vb;
}
insert(s,1,inf);
insert(n,t,inf);
for(i=1;i<=m;i++){
int u=read(),v=read(),ea=read(),eb=read(),ec=read();
insert(u,v,ec);insert(v,u,ec);
insert(s,n+i,ea);insert(n+i,u,inf);insert(n+i,v,inf);
insert(n+m+i,t,eb);insert(u,n+m+i,inf);insert(v,n+m+i,inf);
ans+=ea+eb;
}
while(bfs()) ans-=dfs(s,inf);
printf("%d\n",ans);
return 0;
}