bzoj3118: Orz the MST（线性规划+单纯形法）

传送门

不难发现，对于每一条树边肯定要减小它的权值，对于每一条非树边要增加它的权值

对于每一条非树边\(j\)，他肯定与某些树边构成了一个环，那么它的边权必须大于等于这个环上的所有边

设其中一条边为\(i\)，变化量为\(x\)，那么就要满足\(w_i-x_i\leq w_j+x_j\)，即\(x_i+x_j\geq w_i-w_j\)

然后这就是个线性规划了。因为这线性规划的目标函数要取最小，所以我们把它对偶一下就可以了

//minamoto
#include<bits/stdc++.h>
#define R register
#define Loli true
#define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
    R int res,f=1;R char ch;
    while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
    for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
    return res*f;
}
const int N=1005,M=10005;const double eps=1e-8,inf=1e18;
struct eg{int v,nx,id;}e[N<<1];int head[N],tot;
inline void add(R int u,R int v,R int id){e[++tot]={v,head[u],id},head[u]=tot;}
int n,m,nn,mm,u,v,x,dep[N],fa[N],id[N],U[N],V[N],W[N],F[N],A[N],B[N];
double a[N][M];
void dfs(int u){go(u)if(v!=fa[u])fa[v]=u,id[v]=e[i].id,dep[v]=dep[u]+1,dfs(v);}
void pivot(int l,int e){
	double t=a[l][e];a[l][e]=1;fp(i,0,m)a[l][i]/=t;
	fp(i,0,n)if(i!=l&&fabs(a[i][e])>eps){
		t=a[i][e],a[i][e]=0;
		fp(j,0,m)a[i][j]-=t*a[l][j];
	}
}
void simplex(){
	while(Loli){
		int l=0,e=0;double mn=inf;
		fp(i,1,m)if(a[0][i]>eps){e=i;break;}if(!e)return;
		fp(i,1,n)if(a[i][e]>eps&&a[i][0]/a[i][e]<mn)mn=a[i][0]/a[i][e],l=i;
		pivot(l,e);
	}
}
int main(){
//	freopen("testdata.in","r",stdin);
	nn=read(),mm=read();
	fp(i,1,mm){
		U[i]=read(),V[i]=read(),W[i]=read(),F[i]=read(),A[i]=read(),B[i]=read();
		if(F[i])add(U[i],V[i],i),add(V[i],U[i],i);
	}dfs(1),n=mm;
	fp(i,1,mm)if(F[i])a[i][0]=B[i];
	else{
		a[i][0]=A[i],u=U[i],v=V[i];
		while(u!=v){
			if(dep[u]<dep[v])swap(u,v);
			x=id[u];if(W[x]>W[i])++m,a[x][m]=a[i][m]=1,a[0][m]=W[x]-W[i];
			u=fa[u];
		}
	}simplex();printf("%.0lf\n",-a[0][0]);return 0;
}