#include<stdio.h>
/*思路:显然对于每个地方,
只有一种供水方式就足够了,这样也能保证花费最小,
而每个地方都可以自己挖井,所以是不可能出现无解的情况的,
为了方便思考,我们引入一个虚拟点,把所有自己挖井的都连到这个点,
边权为挖井的花费,而如果i能从j处引水,则从j向i连边,边权为引水的花费,
然后对这个有向图,以虚拟点为根,求最小树形图即可(最小树形图即为有向图的最小生成树)。*/
#include<string.h>
#include<math.h>
#include<stdlib.h>
#define N 1100
#define inf 999999999
struct node {
int x,y,z;
}f[N];
int n;
struct nodee{
int u,v,w;
}edge[N*N];
int manhadun(int x,int y) {
return abs(f[x].x-f[y].x)+abs(f[x].y-f[y].y)+abs(f[x].z-f[y].z);
}
int yong,visit[N],pre[N],fffma[N],id[N];
void addedge(int u,int v,int w) {
edge[yong].u=u;
edge[yong].v=v;
edge[yong++].w=w;
}
int zhuliu(int root) {
int sum=0,i;
while(1) {
for(i=0;i<=n;i++)
fffma[i]=inf;
memset(visit,-1,sizeof(visit));
memset(id,-1,sizeof(id));
for(i=0;i<yong;i++){
int u=edge[i].u,v=edge[i].v;
if(edge[i].w<fffma[v]&&u!=v) {//求所有除根节点外的点,选择一天最小权值的边
pre[v]=u;//记录前驱
fffma[v]=edge[i].w;
}
}
fffma[root]=0;
pre[root]=root;
for(i=0;i<=n;i++){//若有孤立点则无解
if(fffma[i]==inf)return -1;
sum+=fffma[i];
}
int res=0;
for(i=0;i<=n;i++)//所选择的边有没有环和环的个数
if(visit[i]==-1) {
int v=i;
while(visit[v]!=i&&id[v]==-1) {//这个地方需要注意
visit[v]=i;
v=pre[v];
}
if(visit[v]!=i||v==root)continue;
int u;
for(u=pre[v];u!=v;u=pre[u])
id[u]=res;
id[v]=res++;
}
if(res==0) return sum;//如果无环,则直接返回解
for(i=0;i<=n;i++)//将环外的点都加入
if(id[i]==-1)
id[i]=res++;
for(i=0;i<yong;i++) {//缩点出边不变,入边变
edge[i].w-=fffma[edge[i].v];
edge[i].u=id[edge[i].u];
edge[i].v=id[edge[i].v];
}
n=res-1;root=id[root];//相当于建立了一个新图,root也会改变
}
return sum;
}
int main() {
int i,j,k,x,y,z;
while(scanf("%d%d%d%d",&n,&x,&y,&z),n||x||y||z) {
yong=0;
for(i=1;i<=n;i++) {
scanf("%d%d%d",&f[i].x,&f[i].y,&f[i].z);
addedge(0,i,f[i].z*x);//构造的虚拟点0为根节点,权值为每家建井的费用
}
for(i=1;i<=n;i++) {
scanf("%d",&k);
while(k--) {
scanf("%d",&j);
if(j==i)continue;
if(f[j].z>f[i].z)
addedge(i,j,manhadun(i,j)*y+z);
else
addedge(i,j,manhadun(i,j)*y);
}
}
printf("%d\n",zhuliu(0));
}
return 0;
}
转载于:https://www.cnblogs.com/thefirstfeeling/p/4410706.html