[51nod2982] 大逃杀

Posted 2020-11-27 lslzf

问题描述

α?76星人入侵地球，尽管英勇的地球卫士们奋勇杀敌，但战火依然烧到了地球本球上，为了维护人类的利益，联合国决定建立m个避难所，第i个避难所能容纳ci个人，坐标为(pi,qi)。现在地球上还残余着n个国家，第i个国家有bi个人，坐标为(xi,yi)，现在联合国给出了一个初步的疏散方案，即一个n*m的矩阵a，其中ai,j表示第i个国家去第j个避难所的人数（(b_i=sum_{1le jle m}a_{i,j})且(sum_{1le ile 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 <queue>
#include <cstring>
#define N 5002
#define M 50002
using namespace std;
const int inf=1<<30;
int head[N],ver[M*2],nxt[M*2],cap[M*2],cost[M*2],l;
int n,m,s,t,i,j,pre[N],flow[N],dis[N],ans,maxflow,xn[N],yn[N],xm[N],ym[N],p[N];
bool in[N];
int abs(int x)
{
	return (x<0)?-x:x;
}
int read()
{
    char c=getchar();
    int w=0,f=1;
    while(c<‘0‘||c>‘9‘){
    	if(c==‘-‘) f=-1;
    	c=getchar();
	}
    while(c<=‘9‘&&c>=‘0‘){
        w=w*10+c-‘0‘;
        c=getchar();
    }
    return w*f;
}
void insert(int x,int y,int z,int w)
{
    ver[l]=y;
    cap[l]=z;
    cost[l]=w;
    nxt[l]=head[x];
    head[x]=l;
    l++;
    ver[l]=x;
    cost[l]=-w;
    nxt[l]=head[y];
    head[y]=l;
    l++;
}
bool SPFA()
{
    queue<int> q;
    memset(in,0,sizeof(in));
    for(int i=1;i<=t;i++) dis[i]=inf;
    q.push(s);
    dis[s]=0;in[s]=1;
    flow[s]=inf;
    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]>dis[x]+cost[i]&&cap[i]>0){
                dis[y]=dis[x]+cost[i];
                flow[y]=min(flow[x],cap[i]);
                pre[y]=i;
                if(!in[y]){
                    in[y]=1;
                    q.push(y);
                }
            }
        }
        in[x]=0;
    }
    if(dis[t]==inf) return 0;
    return 1;
}
void update()
{
    int x=t;
    while(x!=s){
        int i=pre[x];
        cap[i]-=flow[t];
        cap[i^1]+=flow[t];
        x=ver[i^1];
    }
    ans+=flow[t]*dis[t];
    maxflow+=flow[t];
}
int Dis(int i,int j)
{
	return abs(xn[i]-xm[j])+abs(yn[i]-ym[j])+1;
}
int main()
{
    memset(head,-1,sizeof(head));
    n=read();m=read();
    t=n+m+1;
    for(i=1;i<=n;i++) xn[i]=read(),yn[i]=read();
	for(i=1;i<=m;i++){
		xm[i]=read(),ym[i]=read(),p[i]=read();
		insert(n+i,t,p[i],0);
	}
	int tmp=0;
	for(i=1;i<=n;i++){
		for(j=1;j<=m;j++){
			int a=read();
			tmp+=a*Dis(i,j);
		}
		for(j=1;j<=m;j++) insert(i,n+j,inf,Dis(i,j));
	}
    cout<<tmp<<endl;
    return 0;
}