ccf-201409-4(广度优先+精度陷阱)

问题描述

试题编号：	201409-4
试题名称：	最优配餐
时间限制：	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，每个客户所需要的餐都能被送到。

本题的解题思路时广度优先搜索，同时在广度优先上进行优化

比如：不能一个一个用户单独送，而是同时从分店出发，来查找哪一个分店的线路先到，若找到客户了，就把消耗加到结果上

还有一个就是用户可能在同一个位置，在赋值的时候不是简单的赋值，而是采用+=的形式

注意精度，最后的结果可能int型保存不下，要采用longlong型来解决

以下是代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <vector>
#include <queue>
#include <cstring>
#include <algorithm>
#include <fstream>
#include <map>
using namespace std;

///201409

int n,m,k,d;
int addr[1005][1005],notto[1005][1005],a[1005][1005];
int flag[1005][1005];///flag数组和notto数组可以采用同一个数组
int dx[]={1,0,-1,0};
int dy[]={0,1,0,-1};

long long sum_price=0;

struct Node
{
    int x,y,len;
    Node(int x1,int y1,int len1)
    {
        len=len1;
        x=x1;
        y=y1;
    }
};

bool islegal(int x,int y)
{
    return x>=1&&x<=n&&y>=1&&y<=n&&!notto[x][y];
}

queue <Node> node;

int bfs()
{
    while(!node.empty())
    {
        Node q=node.front();
        node.pop();
        for(int i=0;i<4;i++)
        {
            int x=q.x+dx[i];
            int y=q.y+dy[i];
            int len=q.len;
            if(islegal(x,y)&&!flag[x][y])
            {
                sum_price+=a[x][y]*(len+1);
                flag[x][y]=1;
                node.push(Node(x,y,len+1));
            }
        }
    }
}

int main()
{
    freopen("in.txt","r",stdin);
    memset(addr,0,sizeof(addr));
    memset(notto,0,sizeof(notto));
    memset(a,0,sizeof(a));
    memset(flag,0,sizeof(flag));
    cin>>n>>m>>k>>d;
    for(int i=0;i<m;i++)
    {
        int x,y;
        cin>>x>>y;
        addr[x][y]=1;
        node.push(Node(x,y,0));
    }
    for(int i=0;i<k;i++)
    {
        int x,y,z;
        cin>>x>>y>>z;
        a[x][y]+=z;
    }
    for(int i=0;i<d;i++)
    {
        int x,y;
        cin>>x>>y;
        notto[x][y]=1;
    }
    bfs();
    cout<<sum_price<<endl;
    return 0;
}