AW361 观光奶牛

本文探讨了在特定图上将0/1规划问题转化为乘积形式,并使用SPFA算法进行合理性判断的方法。通过实例代码详细解析了算法实现过程,包括二分查找、图的构建和最短路径寻找等关键步骤。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目地址


易错点:

  • 二分时,l,r,mid的类型都要设为double.
  • 0/1规划问题在该类图上可转化为乘积的形式并使用spfa判断当前解是否合理.
  • 注意区分i和i所代表的点a[i].

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
#define eps 1e-3
using namespace std;
const int MAXN=2e3,MAXM=1e4;
struct Edge{
	int from,to,nxt;
	double w;
}e[MAXM];
int head[MAXN],edgeCnt=1;
void addEdge(int u,int v,double w){
	e[++edgeCnt].from=u;
	e[edgeCnt].to=v;
	e[edgeCnt].w=w;
	e[edgeCnt].nxt=head[u];
	head[u]=edgeCnt;
}
double dis[MAXN];
int n,enterNum[MAXN];
bool inQueue[MAXN];
bool spfa(){
	memset(dis,0x3f,sizeof(dis));
	memset(enterNum,0,sizeof(enterNum));
	queue<int> q;
	for(int i=1;i<=n;i++){
		q.push(i);
		inQueue[i]=1;
	}
	while(!q.empty()){
		if(dis[q.front()]>dis[q.back()])swap(q.front(),q.back());
		int nowPoint=q.front();q.pop();
		inQueue[nowPoint]=0;
		for(int i=head[nowPoint];i;i=e[i].nxt){
			int nowV=e[i].to;
			if(dis[nowV]>dis[nowPoint]+e[i].w){
				dis[nowV]=dis[nowPoint]+e[i].w;
				enterNum[nowV]++;
				if(enterNum[nowV]>n)return 1;
				if(inQueue[nowV])continue;
				q.push(nowV);
				inQueue[nowV]=1;
			}
		}
	}
	return 0;
}
int a[MAXM],b[MAXM],t[MAXM],f[MAXN];
int m;
bool judge(double nowValue){
	memset(head,0,sizeof(head));
	edgeCnt=1;
	for(int i=1;i<=m;i++){
		addEdge(a[i],b[i],nowValue*t[i]-f[a[i]]);
	}
	return spfa();
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&f[i]);
	}
	for(int i=1;i<=m;i++){
		scanf("%d%d%d",&a[i],&b[i],&t[i]);
	}
	double l=0,r=1e3;
	while(r-l>eps){
		double mid=(l+r)/2;
		if(judge(mid))l=mid;
		else r=mid;
	}
	printf("%.2f\n",l);
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值