[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^'';
     while(isdigit(ch=getchar())) x=(((x<<)+x)<<)+(ch^'');
     return x;
 }
 const int N=,M=;
 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=;
         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,-,sizeof head);
     for(register int i=;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,-,sizeof head);
     for(register int i=;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>) {
                 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 ;
 }