傳送門

這道題是個Tarjan+DAG上求最長路的好題啊……來說說我的智障做法。

首先這題很明顯要縮點，這都沒啥，然後我發現要求最長路？

好啊，直接找昨天學的那種，深搜！過了樣例就交。

算法一，期望得分100，實際得分7.

？？？好的，發現這個算的是從一個點出發的最長路，我們要求的是終點還必須是酒館的一條路徑，於是乎有了算法二：從每個點倒著往回跑最長路，遇到出發點就停！

算法二，期望得分100，實際得分87.（好吧其實是沒卡我）

仔細想了一下發現，好像反著跑的話你可以一直不經過起點，這樣不符合要求……

然後想了一下，就又有了算法三：先在正圖上把所以起點能到的點都跑出來，這些之中的點再往回跑！

算法三：期望得分100，實際得分75.（？？？？？）

好吧後來發現，這樣反著往回跑的話，你還是有可能不經過起點。

於是乎最後只好直接用SPFA跑最長路了……這次終於得到了期望的100分……

看一下代碼。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#define rep(i,a,n) for(int i = a;i <= n;i++)
#define per(i,n,a) for(int i = n;i >= a;i--)
#define
 
 enter putchar(‘\n‘)
using namespace std;
typedef long long ll;
const int M = 5000005;
const int N = 10000005;
 
int read()
{
    int ans = 0,op = 1;
    char ch = getchar();
    while(ch < ‘0‘ || ch > ‘9‘)
    {
    if(ch == ‘-‘) op = -1;
    ch = getchar();
    }
    while(ch >=‘0‘ && ch <= ‘
 
9‘)
    {
    ans *= 10;
    ans += ch - ‘0‘;
    ch = getchar();
    }
    return ans * op;
}

struct edge
{
    int next,to,from,v;
}e[M],e1[M];

int head[M],dis[M],n,m,ecnt,ecnt1,low[M],dfn[M],idx,cnt,val[M],sval[M],maxn,x,y,s,p,belong[M];
int stack[M],top,head1[M];
bool bar[M],bar1[M],vis[M],pd[M];

void add(int x,int y)
{
    e[++ecnt].to = y;
    e[ecnt].from = x;
    e[ecnt].next = head[x];
    head[x] = ecnt;
}

void add1(int x,int y,int z)
{
    e1[++ecnt1].to = y;
    e1[ecnt1].from = x;
    e1[ecnt1].v = z;
    e1[ecnt1].next = head1[x];
    head1[x] = ecnt1;
}

void tarjan(int x)
{
    low[x] = dfn[x] = ++idx;
    vis[x] = 1,stack[++top] = x;
    for(int i = head[x];i;i = e[i].next)
    {
    if(!dfn[e[i].to]) tarjan(e[i].to),low[x] = min(low[x],low[e[i].to]);
    else if(vis[e[i].to]) low[x] = min(low[x],dfn[e[i].to]);
    }
    if(dfn[x] == low[x])
    {
    cnt++;
    int p;
    while((p = stack[top--]))
    {
        belong[p] = cnt,vis[p] = 0;
        sval[cnt] += val[p];
        if(bar[p]) bar1[cnt] = 1;
        if(p == x) break;
    }
    }
}

void rebuild()
{
    add1(0,belong[s],sval[belong[s]]);
    rep(i,1,ecnt)
    {
    int r1 = belong[e[i].from],r2 = belong[e[i].to];
    if(r1 != r2) add1(r1,r2,sval[r2]);
    }
}

void spfa()
{
    queue <int> q;
    dis[0] = 0,q.push(0);
    while(!q.empty())
    {
    int k = q.front();q.pop();vis[k] = 0;
    for(int i = head1[k];i;i = e1[i].next)
    {
        if(dis[e1[i].to] < dis[k] + e1[i].v)
        {
        dis[e1[i].to] = dis[k] + e1[i].v;
        if(!vis[e1[i].to]) vis[e1[i].to] = 1,q.push(e1[i].to);
        }
    }
    }
}

int main()
{
    n = read(),m = read();
    rep(i,1,m) x = read(),y = read(),add(x,y);
    rep(i,1,n) val[i] = read();
    s = read();
    rep(i,1,n) if(!dfn[i]) tarjan(i);
    rebuild();
    spfa();
    p = read();
    rep(i,1,p) x = read(),maxn = max(maxn,dis[belong[x]]);
    printf("%d\n",maxn);
    return 0;
}

APIO2009 搶掠計劃