比赛 2024暑假C班集训A 评测结果 AAAAAAAAAA
题目名称 制作人偶 最终得分 100
用户昵称 op_组撒头屯 运行时间 0.606 s
代码语言 C++ 内存使用 4.07 MiB
提交时间 2024-07-10 08:07:23
显示代码纯文本
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
#define ld long double
#define pii pair<int,int>
#define fi first
#define se second
#define pb push_back
const int N=500+5;
const int M=5000+5;
const ld eps=1e-8;
int n,m;
int c[N],tot=0;
int w[M],sum[N];
pii p[M];

struct edge{int to,next;ld c;};
edge eg[M*8];int head[N],ne=1;
inline void adde(int x,int y,ld c){
	eg[++ne]={y,head[x],c},head[x]=ne;
	eg[++ne]={x,head[y],0},head[y]=ne;
}

int S,T;
int dep[N];
queue<int>Q;
bool bfs(){
	while(!Q.empty())Q.pop();
	memset(dep,0,sizeof(dep));
	Q.push(S),dep[S]=1;
	while(!Q.empty()){
		int u=Q.front();Q.pop();
		for (int i=head[u];i;i=eg[i].next){
			int v=eg[i].to;ld c=eg[i].c;
			if (dep[v]||fabs(c)<eps)continue;
			dep[v]=dep[u]+1,Q.push(v);
			if (v==T)return 1;
		}
	}return 0;
}
ld dinic(int u,ld flow){
	if (u==T||fabs(flow)<eps)return flow;
	ld res=flow;
	for (int i=head[u];i;i=eg[i].next){
		int v=eg[i].to;ld c=eg[i].c;
		if (dep[v]!=dep[u]+1||fabs(c)<eps)continue;
		ld cur=dinic(v,min(res,c));
		if(fabs(cur)<eps){dep[v]=0;continue;}
		eg[i].c-=cur,eg[i^1].c+=cur;
		res-=cur;
		if (fabs(res)<eps)break;
	}
	return flow-res;
}
bool check(ld mid){
	memset(head,0,sizeof(head)),ne=1;
	S=n+1,T=n+2;
	for (int i=1;i<=n;i++)adde(S,i,0.5*sum[i]),adde(i,T,mid*c[i]);
	for (int i=1;i<=m;i++)adde(p[i].fi,p[i].se,0.5*w[i]),adde(p[i].se,p[i].fi,0.5*w[i]);
	ld ans=0,flow;
	while(bfs())while(flow=dinic(S,1e8))ans+=flow;
	return ans+eps<tot;
}
int main(){
	freopen ("asiram.in","r",stdin);
	freopen ("asiram.out","w",stdout);
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++)scanf("%d",&c[i]);
	for (int i=1;i<=m;i++){
		scanf("%d%d%d",&p[i].fi,&p[i].se,&w[i]);
		sum[p[i].fi]+=w[i],sum[p[i].se]+=w[i];
		tot+=w[i];
	}
	ld l=0,r=5e5;
	while(fabs(r-l)>eps){
		ld mid=(l+r)/2;
		if (check(mid))l=mid;
		else r=mid;
	}
	printf("%.8Lf\n",l);
	return 0;
}