16.7.17 C 1782. Travel

　给出一个有 个顶点 条边的有向图，对于一条边长度为len的边有两种走法。
　　1、如果a和b可以互达，则走过这条边的时间为len
　　2、如果a和b不可以互达，则走过这条边的时间为2*len
　　现在给出一个k，问，从顶点1到顶点n，满足第二种走法不超过k次的最短时间是多少。

Input

　　第一行有3个整数n,m,k(1<=n<=100,1<=m<=10000,0<=k<=10),表示有n个顶点，m条边。
　　接下来有m行，每行有3个整数xi,yi,leni(1<=xi,yi<=n,1<=leni<=10000)，表示长度为leni的有向边。
　　注意，两个点可能有多条边连接。

Output

　　一行一个整数，表示最短时间。
　　如果没有满足题目条件的路径，则输出-1

Sample Input

7 7 3

1 2 2

1 3 2

2 4 3

4 7 5

3 5 4

5 6 1

6 4 2

Sample Output

20

Data Constraint

Hint

【数据约定】
　　对于30%的数据n<=10,m<=10,
　　对于100%的数据，如题目描述
解题思路：用floyd求出最短路，然后分层sbfa,f[i][j] 表示1到i用了j次所用的时间，这样答案就是f[n][0..k]
代码如下：

#include<cstdio>
#include<iostream>
#include<cstring>
#define F(i,x,y) for(int i=x;i<=y;i++)
#define d(i,x,y) for(int i=x;i>=y;i--)
using namespace std;
const int maxn=120;
int n,m,k;
int a[maxn][maxn],c[maxn][maxn];
bool b[maxn][17];
long long  f[maxn][17];
int d[1000005][3];
int x,y,t;
int p;
long long ans;
void init()
{
    scanf("%d %d %d",&n,&m,&k);
    memset(a,5,sizeof(a));
    memset(c,0,sizeof(c));
    ans=0x7fffffff;
    p=a[1][1];
    F(i,1,m)
    {
        scanf("%d %d %d",&x,&y,&t);
        if(a[x][y]==p)
        {
            a[x][y]=t;
            c[x][0]++;
            c[x][c[x][0]]=y;
        }
        a[x][y]=min(a[x][y],t);
    }
}
void floyed()
{
    F(i,1,n)
     F(j,1,n)
      F(k,1,n) 
        if(i != j && i != k && j != k && a[i][k] != p && a[k][j] != p)
          a[i][j]=min(a[i][j],a[i][k] + a[k][j]);

}
void spfa()
{
    memset(b,true,sizeof(b));

    F(i,1,n)
     F(j,0,10) f[i][j]=0x7fffffff;

    d[1][1]=1;
    d[1][2]=0;
    f[1][0]=0;
    b[1][0]=false;
    int t = 0,w = 1;
    int l,q,r;
    while(t < w)
    {
        t++;
        l = d[t][1];
        q = d[t][2];

        F(i,1,c[l][0])
        {
            r = c[l][i];
            if(a[r][l] != p)
            {
                if(f[l][q] + a[l][r]<f[r][q])
                {
                    f[r][q] = a[l][r] + f[l][q];
                    if(b[r][q])
                    {
                        w++;
                        d[w][1] = r;
                        d[w][2] = q;
                        b[r][q] = false;
                    }
                }
            }
            else 
            {
                if(q == k) continue;
                if(f[l][q]+a[l][r] * 2 < f[r][q+1])
                {
                   f[r][q+1] = a[l][r] * 2 + f[l][q];
                   if(b[r][q+1])
                   {
                      w++;
                      d[w][1] = r;
                      d[w][2] = q + 1;
                      b[r][q+1] = false;
                   }
                }

            }
        }
        b[l][q]=true;
    }
    F(i,0,k) ans=min(ans,f[n][i]);
}
void work()
{
    floyed();
    spfa();
}
void getout()
{
    if(ans==0x7fffffff) ans=-1;

    printf("%lld\n",ans);
}
int main()
{
    freopen("4.in","r",stdin);

    init();
    work();
    getout();
    return 0;
}