题解:
其实就是一个简单的最小割判断是否唯一解。。。
可是我写了一上午还没过。。。T_T
把1-n的最短路上的边提出来做最小割。
然后从s,t分别bfs判断必须在某个割的点。如果有的点没有被bfs到,那么最小割方案不为1。
因为s到它的边满流,它到t的边也满流,哪条边都可以作为割边。
但还是有很多坑点啊!!!一条路两端的权值相同。。。
现在还没过。。。
代码:
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cmath> 4 #include<cstring> 5 #include<algorithm> 6 #include<iostream> 7 #include<vector> 8 #include<map> 9 #include<set> 10 #include<queue> 11 #include<string> 12 #define inf 1000000000 13 #define maxn 100000 14 #define maxm 100000 15 #define eps 1e-10 16 #define ll long long 17 #define pa pair<int,int> 18 #define for0(i,n) for(int i=0;i<=(n);i++) 19 #define for1(i,n) for(int i=1;i<=(n);i++) 20 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 21 #define for3(i,x,y) for(int i=(x);i>=(y);i--) 22 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go) 23 #define mod 1000000007 24 using namespace std; 25 inline int read() 26 { 27 int x=0,f=1;char ch=getchar(); 28 while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} 29 while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();} 30 return x*f; 31 } 32 int n,m,s,t,maxflow,tot=1,a[maxn],u[maxn],v[maxn],w[maxn],head[maxn],cur[maxn],h[maxn]; 33 queue<int>q; 34 ll d[2][maxn]; 35 bool vv[maxn],can[maxn]; 36 struct edge{int go,next,v;}e[maxm]; 37 void add(int x,int y,int v) 38 { 39 cout<<x<<' '<<y<<' '<<v<<endl; 40 e[++tot]=(edge){y,head[x],v};head[x]=tot; 41 e[++tot]=(edge){x,head[y],0};head[y]=tot; 42 } 43 void add2(int x,int y,int v) 44 { 45 e[++tot]=(edge){y,head[x],v};head[x]=tot; 46 e[++tot]=(edge){x,head[y],v};head[y]=tot; 47 } 48 bool bfs() 49 { 50 for(int i=s;i<=t;i++)h[i]=-1; 51 q.push(s);h[s]=0; 52 while(!q.empty()) 53 { 54 int x=q.front();q.pop(); 55 for(int i=head[x];i;i=e[i].next) 56 if(e[i].v&&h[e[i].go]==-1) 57 { 58 h[e[i].go]=h[x]+1;q.push(e[i].go); 59 } 60 } 61 return h[t]!=-1; 62 } 63 int dfs(int x,int f) 64 { 65 if(x==t) return f; 66 int tmp,used=0; 67 for(int i=cur[x];i;i=e[i].next) 68 if(e[i].v&&h[e[i].go]==h[x]+1) 69 { 70 tmp=dfs(e[i].go,min(e[i].v,f-used)); 71 e[i].v-=tmp;if(e[i].v)cur[x]=i; 72 e[i^1].v+=tmp;used+=tmp; 73 if(used==f)return f; 74 } 75 if(!used) h[x]=-1; 76 return used; 77 } 78 void dinic() 79 { 80 maxflow=0; 81 while(bfs()) 82 { 83 for (int i=s;i<=t;i++)cur[i]=head[i];maxflow+=dfs(s,inf); 84 } 85 } 86 inline void dfss(int k,int x) 87 { 88 can[x]=1; 89 cout<<"AAAAA "<<x<<endl; 90 for4(i,x)if(e[i^k].v&&!can[y])dfss(k,y); 91 } 92 void spfa(int k,int s) 93 { 94 for (int i=1;i<=n;i++){vv[i]=0;d[k][i]=inf;} 95 q.push(s);d[k][s]=0;vv[s]=1; 96 while(!q.empty()) 97 { 98 int x=q.front();q.pop();vv[x]=0; 99 for (int i=head[x],y;i;i=e[i].next) 100 if(e[i].v&&d[k][x]+(ll)e[i].v<d[k][y=e[i].go]) 101 { 102 d[k][y]=d[k][x]+(ll)e[i].v; 103 if(!vv[y]){vv[y]=1;q.push(y);} 104 } 105 } 106 } 107 int main() 108 { 109 freopen("input.txt","r",stdin); 110 freopen("output.txt","w",stdout); 111 int T=read(); 112 while(T--) 113 { 114 n=read();m=read(); 115 for1(i,n-1)a[i]=read();a[n]=inf; 116 memset(head,0,sizeof(head));tot=0; 117 for1(i,m){u[i]=read();v[i]=read();w[i]=read();add2(u[i],v[i],w[i]);} 118 spfa(0,1);spfa(1,n); 119 memset(head,0,sizeof(head));tot=1; 120 for1(i,m) 121 { 122 if(d[0][u[i]]+w[i]+d[1][v[i]]==d[0][n])add(u[i],v[i],min(a[u[i]],a[v[i]])); 123 if(d[0][v[i]]+w[i]+d[1][u[i]]==d[0][n])add(v[i],u[i],min(a[u[i]],a[v[i]])); 124 } 125 s=1;t=n; 126 dinic(); 127 memset(can,0,sizeof(can)); 128 dfss(0,s);dfss(1,t); 129 bool flag=0; 130 for1(i,n)for4(j,i) 131 { 132 cout<<i<<' '<<e[j].v<<' '<<e[j].go<<' '<<can[e[j].go]<<' '<<a[i]<<' '<<a[e[j].go]<<endl; 133 if((j&1)==0&&e[j].v==0&&a[i]==a[e[j].go])flag=1; 134 if(!can[e[j].go])flag=1; 135 } 136 printf("%s %d\n",flag?"No":"Yes",maxflow); 137 } 138 return 0; 139 }