洛谷——P1828 香甜的黄油 Sweet Butter

香甜的黄油 Sweet Butter

题目描述

农夫John发现做出全威斯康辛州最甜的黄油的方法：糖。把糖放在一片牧场上，他知道N（1<=N<=500）只奶牛会过来舔它，这样就能做出能卖好价钱的超甜黄油。当然，他将付出额外的费用在奶牛上。

农夫John很狡猾。像以前的Pavlov，他知道他可以训练这些奶牛，让它们在听到铃声时去一个特定的牧场。他打算将糖放在那里然后下午发出铃声，以至他可以在晚上挤奶。

农夫John知道每只奶牛都在各自喜欢的牧场（一个牧场不一定只有一头牛）。给出各头牛在的牧场和牧场间的路线，找出使所有牛到达的路程和最短的牧场（他将把糖放在那）

输入输出格式

输入格式：

 

第一行: 三个数：奶牛数N，牧场数（2<=P<=800），牧场间道路数C(1<=C<=1450)

第二行到第N+1行: 1到N头奶牛所在的牧场号

第N+2行到第N+C+1行： 每行有三个数：相连的牧场A、B，两牧场间距离D（1<=D<=255），当然,连接是双向的

 

输出格式：

 

一行 输出奶牛必须行走的最小的距离和

 

输入输出样例

输入样例#1：

3 4 5
2
3
4
1 2 1
1 3 5
2 3 7
2 4 3
3 4 5

输出样例#1：

8

说明

{样例图形

          P2  
P1 @--1--@ C1
         |
         | 
       5  7  3
         |   
         |     C3
       C2 @--5--@
          P3    P4

} {说明：

放在4号牧场最优

}

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

const long MAXN=9999;
const long maxn=0x7fffffff;

bool visit[MAXN];
long MIN=maxn;
long dis[MAXN],N,P,C,u,v,len,mc;
long head[MAXN],num_edge,mcs[MAXN],temp;
queue<long>que;

struct Edge{
    long pre;
    long to;
    long len;
}edge[MAXN];

void in()
{
    num_edge=0;
    for(long i=1;i<=C;i++){
        cin>>u>>v>>len;
        edge[++num_edge].pre=head[u];
        edge[num_edge].to=v;
        edge[num_edge].len=len;
        head[u]=num_edge;
        swap(u,v);
        edge[++num_edge].pre=head[u];
        edge[num_edge].to=v;
        edge[num_edge].len=len;
        head[u]=num_edge;
    }
}

void spfa(int a){
    que.push(a);
    visit[a]=true;
    while(!que.empty()){    
        long k=que.front();
        for(long j=head[k];j!=-1;j=edge[j].pre){
            long w=edge[j].len;
            if(dis[edge[j].to]>dis[k]+w){
                dis[edge[j].to]=dis[k]+w;
                if(!visit[edge[j].to]){
                    visit[edge[j].to]=true;
                    que.push(edge[j].to);
                }
            }
        }
        visit[k]=false;    
        que.pop();
    }
}

int main()
{
    memset(visit,false,sizeof(visit));
    memset(head,-1,sizeof(head));
    
    while(!que.empty())
    que.pop();
    cin>>N>>P>>C;
    for(int i=1;i<=N;i++){
          cin>>mc;
          ++mcs[mc];
    }
    in();
    for(int i=1;i<=P;i++)
    {
        fill(dis,dis+MAXN,maxn);
        temp=0;
        dis[i]=0;
        spfa(i);
        for(int i=1;i<=P;i++)
            temp+=dis[i]*mcs[i];
        if(temp<MIN)
        MIN=temp;
    }
    cout<<MIN<<endl;
    return 0;
}