题目概述
给出
n
个点,
解题报告
在POJ做分数题简直就是作死,可以WA到飞起来,还是分数规划裸题(搞事啊,为什裸题我现在才做过啊),二分枚举答案
mid
,然后用新边权
wx−mid×t(x,y)
建边。只要存在一个环的边权和大于等于
0
,就说明验证成功。
等等……这怎么找啊……其实边权和大于等于
解题报告
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
typedef double DB;
const int maxn=1000,maxm=5000;
int n,m,v[maxn+5];
int E,lnk[maxn+5],nxt[maxm+5],son[maxm+5],w[maxm+5];
int que[maxn+5],ti,vis[maxn+5],tot[maxn+5];DB dis[maxn+5];
inline int fcmp(DB a,DB b) {if (fabs(a-b)<1e-10) return 0;if (a<b) return -1;return 1;}
#define AM(x) (((x)+1)%maxn)
inline bool check(DB mid)
{
for (int i=1;i<=n;i++) dis[i]=-1e100,tot[i]=0;ti++;
int Head=0,Tail=0;que[++Tail]=1;dis[1]=0;vis[1]=ti;
while (Head!=Tail)
{
int x=que[Head=AM(Head)];vis[x]=0;
for (int j=lnk[x];j;j=nxt[j])
if (dis[x]+v[x]-mid*w[j]>dis[son[j]])
{
dis[son[j]]=dis[x]+v[x]-mid*w[j];
if (vis[son[j]]!=ti)
{
que[Tail=AM(Tail)]=son[j];vis[son[j]]=ti;
if (++tot[son[j]]>n) return true;
}
}
}
return false;
}
#define Add(x,y,z) son[++E]=(y),w[E]=(z),nxt[E]=lnk[x],lnk[x]=E
int main()
{
freopen("program.in","r",stdin);
freopen("program.out","w",stdout);
while (~scanf("%d%d",&n,&m))
{
for (int i=1;i<=n;i++) scanf("%d",&v[i]);E=0;memset(lnk,0,sizeof(lnk));
for (int i=1,x,y,z;i<=m;i++) scanf("%d%d%d",&x,&y,&z),Add(x,y,z);
DB L=0,R=1e6;
for (DB mid=(L+R)/2;fcmp(L,R)<=0;mid=(L+R)/2)
if (check(mid)) L=mid+1e-4; else R=mid-1e-4;
printf("%.2f\n",R);
}
return 0;
}