【CCF】最优配餐

试题名称： 最优配餐

时间限制： 1.0s

内存限制： 256.0MB

问题描述： 问题描述

　　栋栋最近开了一家餐饮连锁店，提供外卖服务。随着连锁店越来越多，怎么合理的给客户送餐成为了一个急需解决的问题。

　　栋栋的连锁店所在的区域可以看成是一个n×n的方格图（如下图所示），方格的格点上的位置上可能包含栋栋的分店（绿色标注）或者客户（蓝色标注），有一些格点是不能经过的（红色标注）。

　　方格图中的线表示可以行走的道路，相邻两个格点的距离为1。栋栋要送餐必须走可以行走的道路，而且不能经过红色标注的点。



　　送餐的主要成本体现在路上所花的时间，每一份餐每走一个单位的距离需要花费1块钱。每个客户的需求都可以由栋栋的任意分店配送，每个分店没有配送总量的限制。

　　现在你得到了栋栋的客户的需求，请问在最优的送餐方式下，送这些餐需要花费多大的成本。

输入格式

　　输入的第一行包含四个整数n, m, k, d，分别表示方格图的大小、栋栋的分店数量、客户的数量，以及不能经过的点的数量。

　　接下来m行，每行两个整数xi, yi，表示栋栋的一个分店在方格图中的横坐标和纵坐标。

　　接下来k行，每行三个整数xi, yi, ci，分别表示每个客户在方格图中的横坐标、纵坐标和订餐的量。（注意，可能有多个客户在方格图中的同一个位置）

　　接下来d行，每行两个整数，分别表示每个不能经过的点的横坐标和纵坐标。

输出格式

　　输出一个整数，表示最优送餐方式下所需要花费的成本。

样例输入

10 2 3 3

1 1

8 8

1 5 1

2 3 3

6 7 2

1 2

2 2

6 8

样例输出

29

评测用例规模与约定

　　前30%的评测用例满足：1<=n <=20。

　　前60%的评测用例满足：1<=n<=100。

　　所有评测用例都满足：1<=n<=1000，1<=m, k, d<=n^2。可能有多个客户在同一个格点上。每个客户的订餐量不超过1000，每个客户所需要的餐都能被送到。

答题以来第一次出现了运行超时的错误。原因是用广度优先搜索的时候读取一个计算一个，而不是读取一个计算四个。而且一旦所有的客户都送到餐以后就可以退出广度搜索了。

最后结果没有用long long(_int64)存储导致只有80分。改成long long以后通过。

至此201312 - 201503的CCF前四道题都解决了。估计目前自己还没达到第五题的水平。

#include<iostream>
#include<queue>
using namespace std;

class P{
public:
int x;
int y;
P(){
x = -1;
y = -1;
}
P(int x,int y){
this->x = x;
this->y = y;
}
};

class T:public P{
public:
int z;
T(){
P();
z = -1;
}
T(int x,int y,int z){
this->x =x;
this->y =y;
this->z = z;
}
};

int n,m,k,r;//方格图的大小、栋栋的分店数量、客户的数量，以及不能经过的点的数量
int vis[1001][1001] ={0};
queue<T> q;
long long cost = 0;

//(x+1,y) (x,y+1) (x,y-1) (x-1,y)

void Z(){
T p;
int cnt = 0;
int x,y,z;
while(!q.empty()){
p = q.front();
q.pop();
x = p.x + 1;
y = p.y;
z = p.z + 1;
if(vis[x][y]!=-1 && x>=1 && x<=n && y>=1 && y<=n){
if(vis[x][y]>0 ){
cost += vis[x][y] * z;
cnt++;
if(cnt==k){
return;
}
}
vis[x][y] = -1;
q.push(T(x,y,z));
}
x = p.x;
y = p.y + 1;
if(vis[x][y]!=-1 && x>=1 && x<=n && y>=1 && y<=n){
if(vis[x][y]>0 ){
cost += vis[x][y] * z;
cnt++;
if(cnt==k){
return;
}
}
vis[x][y] = -1;
q.push(T(x,y,z));
}
x = p.x - 1;
y = p.y;
if(vis[x][y]!=-1 && x>=1 && x<=n && y>=1 && y<=n){
if(vis[x][y]>0 ){
cost += vis[x][y] * z;
cnt++;
if(cnt==k){
return;
}
}
vis[x][y] = -1;
q.push(T(x,y,z));
}
x = p.x;
y = p.y - 1;
if(vis[x][y]!=-1 && x>=1 && x<=n && y>=1 && y<=n){
if(vis[x][y]>0 ){
cost += vis[x][y] * z;
cnt++;
if(cnt==k){
return;
}
}
vis[x][y] = -1;
q.push(T(x,y,z));
}
}
}

int main(){
cin>>n>>m>>k>>r;
int i,x,y,z;
for(i=-1;++i<m;){
cin>>x>>y;
q.push(T(x,y,0));
}
for(i=-1;++i<k;){
cin>>x>>y>>z;
if( vis[x][y] != 0 ){
vis[x][y] += z;
}
else{
vis[x][y] = z;
}
}
for(i=-1;++i<r;){
cin>>x>>y;
vis[x][y] = -1;
}
Z();
cout<<cost<<endl;
//	system("pause");
return 0;
}