题目描述:
给定一个 n个点m 条边的带权无向连通图 ,以及一个大小为k 的关键点集合S 。有个人要从点s走到点t,现在可以对所有边加上一个非负整数a,问最大的a,使得加上a后,满足:s到t的最短路长度=s到t且只能经过S中的点的最短路长度。
题目分析:
暴力
记x为只经过关键点的最短路长度,其路径条数为n
记y为可经过任意点的最短路长度,其路径条数为m
tip:路径条数意思这里指 覆盖最短路的边数
显然全部加上a之后最短路是不变的
也就是说我们要求这个东西
$$a*n+x=a*m+y$$
$$a=(\frac{y-x}{n-m})_{max}$$
所以我们怎么求这个东西呢
当然是D(bao)P(li)咯
假设f[i,j],g[i,j]表示到达i这个点,走过了j条边的最短路径
其中f[i,j]只经过关键点,g[i,j]经过任意一点,这两个数组暴力DP就能求出来
然后取出最大的a就好了
关于无解的情况:上面那个式子的值<0显然是不行的,m=n显然是取任意的(当然实现的话就不用这样了判断,这样讲只是方便理解)
代码:
1 //许多大佬貌似喜欢在代码上面写自己id,那我这个蒟蒻也就从此开始——跟啦! 2 //好吧Koko不是重点 3 //LevenKoko 4 #pragma GCC optimize("Ofast") 5 #include<bits/stdc++.h> 6 #define int long long 7 const int inf=1e18; 8 using namespace std; 9 inline int read(){ 10 int ans=0,f=1;char chr=getchar(); 11 while(!isdigit(chr)){if(chr=='-')f=-1;chr=getchar();} 12 while(isdigit(chr)) {ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();} 13 return ans*f; 14 }const int M=1e4+5,N=1e3+5; 15 int n,m,s,t,ans; 16 int head[N],nxt[M<<1],ver[M<<1],val[M<<1],tot,f[N][N],g[N][N],a[M],k; 17 inline void add(int x,int y,int z){ver[++tot]=y,nxt[tot]=head[x],val[tot]=z,head[x]=tot;} 18 inline void Clear_All(){ 19 memset(head,0,sizeof(head)); 20 memset(a,0,sizeof(a)); 21 for(int i=1;i<=n;i++)for(int j=0;j<=n;j++)f[i][j]=g[i][j]=inf; 22 tot=f[s][0]=g[s][0]=0;ans=-1; 23 } 24 inline void cmin(int &x,int y){return (void)(x=(x>y)?y:x);} 25 inline void cmax(int &x,int y){return (void)(x=(x>y)?x:y);} 26 inline void DP1(){ 27 for(int i=0;i<n;i++) 28 for(int j=1;j<=n;j++) 29 if(f[j][i]<inf&&a[j]) 30 for(int k=head[j];k;k=nxt[k]) 31 cmin(f[ver[k]][i+1],f[j][i]+val[k]); 32 } 33 inline void DP2(){ 34 for(int i=0;i<n;i++) 35 for(int j=1;j<=n;j++) 36 if(g[j][i]<inf) 37 for(int k=head[j];k;k=nxt[k]) 38 cmin(g[ver[k]][i+1],g[j][i]+val[k]); 39 } 40 inline void Check_Get(){ 41 for(int i=1;i<=n;i++){ 42 if(f[t][i]==inf) continue; 43 int j; 44 for(j=1;j<=i;j++) 45 if(g[t][j]<f[t][i]) break; 46 if(j<=i) continue; 47 for(j=1;j<i;j++) if(g[t][j]!=inf) break; 48 if(j>=i){ 49 ans=inf; break; 50 } 51 int cur=inf; 52 for(j=1;j<i;j++) 53 cur=min(cur,(-f[t][i]+g[t][j])/(i-j)); 54 for(j=i+1;j<=n;j++) 55 if(f[t][i]+i*cur>g[t][j]+j*cur) break; 56 if(j<=n) continue; 57 ans=max(ans,cur); 58 } 59 } 60 signed main(){ 61 // freopen("ernd.in","r",stdin); 62 // freopen("ernd.out","w",stdout); 63 int T=read(); 64 while(T--){ 65 n=read(),m=read(),s=read(),t=read();Clear_All(); 66 for(int i=1,x,y,z;i<=m;i++)x=read(),y=read(),z=read(),add(x,y,z),add(y,x,z); 67 k=read(); 68 for(int i=1,x;i<=k;i++) x=read(),a[x]=1; 69 DP1();DP2(); 70 Check_Get(); 71 if(ans==-1){puts("Impossible");continue;} 72 if(ans==inf){puts("Infinity");continue;} 73 printf("%lld\n",ans); 74 } 75 return 0; 76 }