洛谷 P2656 采蘑菇（Tarjan缩点，dp）

传送门

---

解题思路

对于每一条路，很显然只有两种情况，一种是走一次，一种是走无限次，而第二种情况的条件是这条边在一个环中。

所以先用Tarjan缩点，每个点的权值更新为这个环上所有边一直摘的蘑菇数的和。

然后跑一边既有边权也有点权的DAG上的dp即可。

//数据有锅，后四个点答案有误，用double存实际上错了，但是标程是double，为了过题只能改成double了。

具体原因请看这篇博客：关于float、double小项问题

AC代码

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<stack>
using namespace std;
const int maxn=80008;
const int maxm=200005;
int n,m,scc_cnt,cnt,p[maxn][2],dp[maxn],w[maxn];
int sscc[maxn],dfn[maxn],low[maxn],times,num[maxn],ss;
stack<int> s;
struct node{
    int v,next,value,sum;
}e[maxm][2];
void insert(int u,int v,int value,int sum,int id){
    cnt++;
    e[cnt][id].v=v;
    e[cnt][id].value=value;
    e[cnt][id].sum=sum;
    e[cnt][id].next=p[u][id];
    p[u][id]=cnt;
}
void dfs(int u){
    dfn[u]=low[u]=++times;
    s.push(u);
    for(int i=p[u][0];i!=-1;i=e[i][0].next){
        int v=e[i][0].v;
        if(!dfn[v]){
            dfs(v);
            low[u]=min(low[u],low[v]);
        }else if(num[v]==0){
            low[u]=min(low[u],dfn[v]);
        }
    }
    if(dfn[u]==low[u]){
        ++scc_cnt;
        while(1){
            int x=s.top();
            s.pop();
            num[x]=scc_cnt;
            if(x==u) break;
        }
    }
}
void dfs2(int u){
    if(dp[u]) return;
    dp[u]=sscc[u];
    for(int i=p[u][1];i!=-1;i=e[i][1].next){
        int v=e[i][1].v;
        dfs2(v);
        dp[u]=max(dp[u],dp[v]+e[i][1].value+sscc[u]);
    }
}
int main()
{
    memset(p,-1,sizeof(p));
    cin>>n>>m;
    for(int i=1;i<=m;i++){
        int u,v,value,sum=0;
        double p;
        scanf("%d%d%d%lf",&u,&v,&value,&p);
        int vv=value;
        while(value>0){
            sum+=value;
            value=floor((double)value*p); 
        }
        insert(u,v,vv,sum,0);
    }
    cin>>ss;
    for(int i=1;i<=n;i++) if(!dfn[i]) dfs(i);
    cnt=0;
    for(int u=1;u<=n;u++){
        for(int i=p[u][0];i!=-1;i=e[i][0].next){
            int v=e[i][0].v;
            if(num[u]==num[v]){
                sscc[num[u]]+=e[i][0].sum;
            }else{
                insert(num[u],num[v],e[i][0].value,0,1);
            }
        }
    }
    dfs2(num[ss]);
    cout<<dp[num[ss]];
    return 0;
}