[Luogu2656]采蘑菇

题目大意：
　　给你一个有向图，每条边有一个边权w以及恢复系数k，
　　你从s点出发乱走，经过某条边时会获得相应的收益w，而当第二次经过这条边时相应的收益为w*k下取整。
　　问你最大能获得的收益为多少？

思路：
　　缩点+DP。
　　首先跑一下Tarjan（只要从s开始跑，因为没跑到的地方肯定和答案没关系）。
　　对于每个强连通分量，我们算一下经过这个强联通分量能获得的总收益sum（就是拼命在这上面绕圈圈）。
　　把原图缩为一个DAG，然后就可以DP了。
　　设当前点为i，后继结点为j，边权为w，j的SCC的总收益为sum[j]，转移方程为f[j]=max{f[i]+w+sum[j]}。
　　当然也只要从s开始DP即可。

#include<cmath>
#include<stack>
#include<queue>
#include<cstdio>
#include<cctype>
#include<vector>
#include<cstring>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int N=80001,M=200000;
struct Edge {
    int from,to,w;
    double k;
    int next;
};
Edge e[M];
int head[N];
inline void add_edge(const int &u,const int &v,const int &w,const double &k,const int &i) {
    e[i]=(Edge){u,v,w,k,head[u]};
    head[u]=i;
}
int s;
int dfn[N],low[N],scc[N],cnt,id,sum[N],in[N];
bool ins[N];
std::stack<int> stack;
void tarjan(const int &x) {
    dfn[x]=low[x]=++cnt;
    stack.push(x);
    ins[x]=true;
    for(int i=head[x];~i;i=e[i].next) {
        const int &y=e[i].to;
        if(!dfn[y]) {
            tarjan(y);
            low[x]=std::min(low[x],low[y]);
        } else if(ins[y]) {
            low[x]=std::min(low[x],dfn[y]);
        }
    }
    if(dfn[x]==low[x]) {
        id++;
        int y=0;
        while(y!=x) {
            y=stack.top();
            stack.pop();
            ins[y]=false;
            scc[y]=id;
        }
    }
}
int ans,f[N];
std::queue<int> q;
inline void dp() {
    q.push(scc[s]);
    f[scc[s]]=sum[scc[s]];
    while(!q.empty()) {
        const int x=q.front();
        q.pop();
        ans=std::max(ans,f[x]);
        for(register int i=head[x];~i;i=e[i].next) {
            const int &y=e[i].to;
            f[y]=std::max(f[y],f[x]+e[i].w+sum[y]);
            if(!--in[y]) q.push(y);
        }
    }
}
int main() {
    const int n=getint(),m=getint();
    memset(head,-1,sizeof head);
    for(register int i=0;i<m;i++) {
        const int u=getint(),v=getint(),w=getint();
        double k;
        scanf("%lf",&k);
        add_edge(u,v,w,k,i);
    }
    s=getint();
    tarjan(s);
    memset(head,-1,sizeof head);
    for(register int i=0;i<m;i++) {
        const int &u=e[i].from,&v=e[i].to;
        if(!dfn[u]||!dfn[v]) continue; 
        if(scc[u]==scc[v]) {
            const double &k=e[i].k;
            int w=e[i].w;
            while(w>0) {
                sum[scc[u]]+=w;
                w=floor(w*k);
            }
        } else {
            in[scc[v]]++;
            add_edge(scc[u],scc[v],e[i].w,e[i].k,i);
        }
    }
    dp();
    printf("%d\n",ans);
    return 0;
}