求助精度问题
  • 板块P2656 采蘑菇
  • 楼主junee
  • 当前回复2
  • 已保存回复2
  • 发布时间2024/11/13 14:50
  • 上次更新2024/11/13 18:20:37
查看原帖
求助精度问题
757101
junee楼主2024/11/13 14:50
#include<iostream>
#include<iomanip>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<random>
#include<chrono>
#include<queue>
using namespace std;
mt19937_64 rnd(chrono::steady_clock::now().time_since_epoch().count());
typedef long long LL;
const int N=2e5+10;
int n,m,s;
int h[N],hh[N],e[N*2],ne[N*2],w[N*2],idx;
long double t[N*2];
int dfn[N],low[N],timestamp;
int stk[N],tt,instk[N];
int id[N],sz[N],din[N],cnt;
int f[N];
queue<int>q;
void add(int h[],int a,int b,int c,long double d){
	e[idx]=b,ne[idx]=h[a],w[idx]=c,t[idx]=d,h[a]=idx++;
}
void tarjan(int u){
	dfn[u]=low[u]=++timestamp;
	stk[++tt]=u,instk[u]=1;
	for(int i=h[u];~i;i=ne[i]){
		int j=e[i];
		if(!dfn[j]){
			tarjan(j);
			low[u]=min(low[u],low[j]);
		}
		else if(instk[j])low[u]=min(low[u],dfn[j]);
	}
	if(dfn[u]==low[u]){
		int y;cnt++;
		do{
			y=stk[tt--];
			instk[y]=0;
			id[y]=cnt;
		}while(y!=u);
	}
}
int get(int x,long double t){
	int res=0;
	while(x){
		res+=x;
		x*=t;
	}
	return res;
}
void dp(){
	f[s]=sz[s];
	q.push(s);
	while(q.size()){
		int ver=q.front();q.pop();
		for(int i=hh[ver];~i;i=ne[i]){
			int j=e[i];
			f[j]=max(f[j],f[ver]+w[i]+sz[j]);
			q.push(j);
		}
	}
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	memset(h,-1,sizeof h);
	memset(hh,-1,sizeof hh);
	cin>>n>>m;
	for(int i=1;i<=m;i++){
		int a,b,c;long double d;
		cin>>a>>b>>c>>d;
		add(h,a,b,c,d);
	}
	cin>>s;
	for(int i=1;i<=n;i++){
		if(!dfn[i])tarjan(i);
	}
	for(int u=1;u<=n;u++){
		for(int i=h[u];~i;i=ne[i]){
			int j=e[i];
			int a=id[u],b=id[j];
			if(a==b)sz[a]+=get(w[i],t[i]);
			else add(hh,a,b,w[i],0);
		}
	}
//	for(int i=1;i<=n;i++)cout<<id[i]<<'\n';
	s=id[s];
	dp();
	int ans=0;
	for(int i=1;i<=cnt;i++)ans=max(ans,f[i]);
	cout<<ans;
	return 0;
}

这份代码开 double 过不了:(,有人知道为啥要开 long double 精度才够吗。

2024/11/13 14:50
加载中...