码迷,mamicode.com
首页 > 其他好文 > 详细

[51nod2982] 大逃杀

时间:2020-04-08 22:26:43      阅读:98      评论:0      收藏:0      [点我收藏+]

标签:using   def   abs   pre   代码   stream   方案   name   efi   

问题描述

α?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 <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;
}

[51nod2982] 大逃杀

标签:using   def   abs   pre   代码   stream   方案   name   efi   

原文地址:https://www.cnblogs.com/LSlzf/p/12663067.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!