noip2012 文化之旅（深搜，最优性剪枝）

P2070 [NOIP2012P4]文化之旅

时间: 1000ms / 空间: 131072KiB / Java类名: Main

背景

NOIP 2012 普及组题4

描述

有一位使者要游历各国，他每到一个国家，都能学到一种文化，但他不愿意学习任何一种文化超过一次（即如果他学习了某种文化，则他就不能到达其他有这种文化的国家）。不同的国家可能有相同的文化。不同文化的国家对其他文化的看法不同，有些文化会排斥外来文化（即如果他学习了某种文化，则他不能到达排斥这种文化的其他国家）。
现给定各个国家间的地理关系，各个国家的文化，每种文化对其他文化的看法，以及这位使者游历的起点和终点（在起点和终点也会学习当地的文化），国家间的道路距离，试求从起点到终点最少需走多少路。

输入格式

第一行为五个整数 N，K，M，S，T，每两个整数之间用一个空格隔开，依次代表国家个数（国家编号为 1 到 N），文化种数（文化编号为 1 到 K），道路的条数，以及起点和终点的编号（保证 S 不等于 T）；
第二行为 N 个整数，每两个整数之间用一个空格隔开，其中第 i 个数 Ci，表示国家 i 的文化为 Ci。

接下来的 K 行，每行 K 个整数，每两个整数之间用一个空格隔开，记第 i 行的第 j 个数为 aij，aij= 1 表示文化 i 排斥外来文化 j（i 等于 j 时表示排斥相同文化的外来人），aij= 0 表示不排斥（注意 i 排斥 j 并不保证 j 一定也排斥 i）。

接下来的 M 行，每行三个整数 u，v，d，每两个整数之间用一个空格隔开，表示国家 u

与国家 v 有一条距离为 d 的可双向通行的道路（保证 u 不等于 v，两个国家之间可能有多条道路）。

输出格式

输出只有一行，一个整数，表示使者从起点国家到达终点国家最少需要走的距离数（如果无解则输出-1）。

测试样例1

输入

样例1：
2 2 1 1 2
1 2
0 1
1 0
1 2 10
样例2：
2 2 1 1 2
1 2
0 1
0 0
1 2 10

输出

样例1：
-1
样例2：
10

备注

对于样例1：

由于到国家2必须要经过国家1，而国家2的文明却排斥国家1的文明，所以不可能到国家2。

对于样例2：

路线为1->2。

【数据范围】

对于 20%的数据，有 2≤N≤8，K≤5；

对于 30%的数据，有 2≤N≤10，K≤5；

对于 50%的数据，有 2≤N≤20，K≤8；

对于 70%的数据，有 2≤N≤100，K≤10；

对于 100%的数据，有 2≤N≤100，1≤K≤100，1≤M≤N2，1≤ki≤K，1≤u, v≤N，1≤d≤1000， S≠T，1 ≤S, T≤N。

解析：最优性剪枝。用floyd计算出任意两点间的可能最短距离（不考虑同色或者排斥），假设搜索到点 i ，已走len，若len+map[i][t]>=ans,那么就没有继续搜索的必要了。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 100
using namespace std;int ans,n,p,m,s,t;
int c[maxn+10],map[maxn+10][maxn+10],q[maxn+10];
bool b[maxn+10][maxn+10],flag[maxn+10];void readdata()
{int i,j,k,u,v,d;scanf("%d%d%d%d%d",&n,&p,&m,&s,&t);for(i=1;i<=n;i++)scanf("%d",&c[i]);for(i=1;i<=p;i++)for(j=1;j<=p;j++)scanf("%d",&k),b[i][j]=k;for(i=1;i<=p;i++)b[i][i]=1;if(b[c[t]][c[s]]){printf("-1\n");exit(0);}memset(map,90,sizeof(map));    for(i=1;i<=m;i++){scanf("%d%d%d",&u,&v,&d);map[u][v]=map[v][u]=d;}for(i=1;i<=n;i++)map[i][i]=0;
}void dfs(int x,int len)
{if(x==t){ans=min(ans,len);return;}int i,j,k,l,r;for(i=1;i<=n;i++)if(!flag[i] && ans-len-map[x][i]>map[i][t]){l=q[0]+1;q[++q[0]]=i,flag[i]=1;for(j=1;j<=n;j++)if(!flag[j] && b[c[j]][c[i]])q[++q[0]]=j,flag[j]=1;r=q[0];dfs(i,len+map[x][i]);for(j=l;j<=r;j++)flag[q[j]]=0;q[0]=l-1;    }
}void work()
{int i,j,k;for(k=1;k<=n;k++)for(i=1;i<=n;i++)if(!b[c[k]][c[i]])for(j=1;j<=n;j++)if(!b[c[j]][c[k]] && !b[c[j]][c[i]])if(map[i][j]-map[i][k]>map[k][j])map[i][j]=map[i][k]+map[k][j];for(i=1;i<=n;i++)if(i!=s && i!=t)if(b[c[i]][c[s]] || b[c[t]][c[i]])q[++q[0]]=i,flag[i]=1;ans=map[maxn+5][maxn+5];q[++q[0]]=s,flag[s]=1,dfs(s,0);  if(ans==map[maxn+5][maxn+5])printf("-1\n");else printf("%d\n",ans);
}int main()
{readdata();work();return 0;
}