【CodeChef-FLYDIST】Flight Distance（线性规划）（对偶原理）（单纯形）

传送门

---

题解：

要求两点之间的连边不大于最短路？？？

也就是说这条边的权值就是最短路的权值。

除此以外还要求不能够继续松弛，否则就不是最短路。

设$w_{uv}$表示$u,v$之间连边的原边权，$t_{uv}^{+}$表示修改边权的增加量，$t_{uv}^{-}$表示修改边权的减少量，$d_{uv}$表示新图中的$u,v$最短路大小，则我们得到如下线性规划：

$$\begin{array}{rlrlrlrl}limits:& & & & & & d_{uv}-t_{uv}^{+}+t_{uv}^{-}& \ge w_{uv}\\ & & & & & & -d_{uv}+t_{uv}^{+}-t_{uv}^{-}& \ge -w_{uv}\\ & & & & & & t_{uv}^{+}-t_{uv}^{-}+d_{uj}-d_{vj}& \ge -w_{uv}\\ minimize:& & & & & & \sum t_{uv}^{+}+t_{uv}^{-}& \end{array}$$

直接做需要初始化找解，对偶一下就可以直接上单纯形了。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

typedef long double ld;
cs int N=13,M=47;

namespace LP{
	int n,m;
	cs ld eps=1e-9,INF=1e17;
	ld a[3*M][2*M+2*M*N];
	
	inline void pivot(int l,int e){
		ld t=a[l][e];a[l][e]=1;
		for(int re i=0;i<=n;++i)a[l][i]/=t;
		for(int re i=0;i<=m;++i)if(i!=l&&fabs(a[i][e])>eps){
			t=a[i][e],a[i][e]=0;
			for(int re j=0;j<=n;++j)a[i][j]-=t*a[l][j];
		}
	}
	inline void simplex(){
		while(true){
			int l=0,e=0;ld mn=INF;
			for(int re j=1;j<=n;++j)if(a[0][j]>eps){e=j;break;}
			if(!e)break;
			for(int re i=1;i<=m;++i)
			if(a[i][e]>eps&&a[i][0]/a[i][e]<mn)mn=a[i][0]/a[i][e],l=i;
			pivot(l,e);
		}
	}
	inline void solve(){
		simplex();
		ld val=-a[0][0];
		for(int re i=1;;++i)if(val*i-(ll)(val*i+eps)<eps){
			printf("%d/%d",int(val*i+eps),i);
			return ;
		}
	}
}
using LP::a; 

int n,m,&mm=LP::n,&nn=LP::m;
inline int d(int i,int j){
	if(i<j)std::swap(i,j);
	static int id[N][N];
	return id[i][j]?id[i][j]:id[i][j]=++nn;
}
inline int tp(int i,int j){
	if(i<j)std::swap(i,j);
	static int id[N][N];
	return id[i][j]?id[i][j]:id[i][j]=++nn;
}
inline int tm(int i,int j){
	if(i<j)std::swap(i,j);
	static int id[N][N];
	return id[i][j]?id[i][j]:id[i][j]=++nn;
}
signed main(){
#ifdef zxyoi
	freopen("flydist.in","r",stdin);
#endif
	scanf("%d%d",&n,&m);
	for(int re i=1;i<=m;++i){
		int u,v,w;scanf("%d%d%d",&u,&v,&w);++u,++v;
		int i1=tp(u,v),i2=tm(u,v),i3=d(u,v);
		a[i1][0]=a[i2][0]=1;
		a[0][++mm]=w,a[i3][mm]=1,a[i1][mm]=-1,a[i2][mm]=1;
		a[0][++mm]=-w,a[i3][mm]=-1,a[i1][mm]=1,a[i2][mm]=-1;
		for(int re j=1;j<=n;++j)if(j!=u&&j!=v){
			a[0][++mm]=-w,a[i1][mm]=1,a[i2][mm]=-1,a[d(v,j)][mm]=1,a[d(u,j)][mm]=-1;
			a[0][++mm]=-w,a[i1][mm]=1,a[i2][mm]=-1,a[d(u,j)][mm]=1,a[d(v,j)][mm]=-1;
		}
	}
	LP::solve();
	return 0;
}