【GDOI 2016 Day2】第一题 SigemaGO

题目

分析

拆点连边+spfa。
首先把图分成2lim+1层，也就是每个点拆成2lim+1个点。
如果a和b之间、b和c有一条有向边，那么连边(k,a)-->(k+1,b)，(k+1,b)-->(k+2,c)（k=1、3、5、···、2lim+1，是当前点所在的层数），这两条边的权值和是l。也就是说当走了(k,a)-->(k+1,b)-->(k+2,c)这条路线时，就是抄了一次近道。
注意：在spfa中，当走到的点在第2、4、6、···、2lim层时，就只能往上走，因为当前点一定在抄近道的过程中。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
using namespace std;
int next[800000],last[800000],to[800000],dis[13][200000],n,m,lim,l,tot,v[800000],ans,d[10000000][2];
bool bz[200000];
int bj(int x,int y,int z)
{
    next[++tot]=last[x];
    last[x]=tot;
    to[tot]=y;
    v[tot]=z;
}
int spfa()
{
    int i,j,head=0,tail=1,k;
    d[1][0]=0;
    d[1][1]=1;
    dis[0][1]=0;
    while(head<tail)
    {
        k=++head;
        bz[d[k][0]*n+d[k][1]]=true;
        for(i=last[d[k][1]+n*d[k][0]];i;i=next[i])
        {
            j=to[i];
            int x=(j-1)/n;
            if(dis[x][(j-1)%n+1]>dis[d[k][0]][d[k][1]]+v[i])
            {
                dis[x][(j-1)%n+1]=dis[d[k][0]][d[k][1]]+v[i];
                if(bz[j])
                {
                    bz[j]=false;
                    d[++tail][0]=x;
                    d[tail][1]=(j-1)%n+1;
                }
            }
        }
    }
}
int main()
{
    freopen("sigemago.in","r",stdin);
    freopen("sigemago.out","w",stdout);
    scanf("%d%d%d%d",&n,&m,&l,&lim);
    for(int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        for(int j=1;j<=lim+1;j++)
        {
            bj(x+n*(j*2-2),y+n*(j*2-2),z);
        }
        for(int j=1;j<=lim*2;j++)
        {
            bj(x+n*(j-1),y+n*j,l*(j%2));
        }
    }
    memset(dis,60,sizeof(dis));
    memset(bz,true,sizeof(bz));
    spfa();
    ans=maxlongint;
    for(int i=1;i<=lim+1;i++)
    {
        if(ans>dis[i*2-2][n]) ans=dis[i*2-2][n];
    }
    if(ans>=dis[0][0])  printf("-1\n");
    else
    printf("%d\n",ans);
}