BZOJ 1003: [ZJOI2006]物流运输(spfa+dp)

传送门：http://www.lydsy.com/JudgeOnline/problem.php?id=1003
这题真是一个想法题啊！
从这题中又获得一些思考。算法只是用来解决问题的工具，而不是用来套题目的解题模板，看到一个题目以后，得到一些想法，要用算法来实现想法来解决题目。而不是看到题目以后，从大脑里枚举算法来套。以前大部分时间，我都在干这种蠢事，缺少思考。
看到这题以后，还是完全没有思路，不知道怎么做，看了别人的题解才看懂怎么做。
有一个地方很值得注意，就是点数<=20，这个地方要是没注意到的话，就很难想到这题怎么做。因为点数少，所以求一次最短路非常快，所以就可以n^2处理出每个时间区间中，起点到终点的距离。（以前都以为最短路不快，而没有意识到，点数少的时候，最短路也是飞快的）。处理出这些距离后，就可以dp了。
状态转移方程：dp[i]=min(dp[i],dp[j]+cost[j+1][i]*(i-j)+k);

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <iostream>
#include <vector>
#include <map>
#include <queue>
#include <ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
#define PB push_back
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define calm (l+r)>>1
const ll INF=(ll)1e9+7;

struct EE{
    int to,v,next;
    EE(){}
    EE(int to,int v,int next):to(to),v(v),next(next){}
}edge[100000];
struct PP{
    int id,l,r;
    PP(){}
    PP(int id,int l,int r):id(id),l(l),r(r){}
}block[100000];
int n,m,k,e,q,Ecnt;
int head[22],d[22];
ll cost[111][111],dp[111];
bool vis[22];

inline void addedge(int from,int to,int v){
    edge[Ecnt]=EE(to,v,head[from]);
    head[from]=Ecnt++;
    edge[Ecnt]=EE(from,v,head[to]);
    head[to]=Ecnt++;
}
bool judge(int s,int t,int l,int r){
    for(int i=0;i<q;i++)if(block[i].id==s||block[i].id==t){
        if(l<=block[i].r&&r>=block[i].l)return false;
    }
    return true;
}
int spfa(int l,int r){
    queue<int> Q;
    fill(d,d+m+1,INF);
    memset(vis,false,sizeof vis);
    Q.push(1);d[1]=0;vis[1]=true;
    while(!Q.empty()){
        int s=Q.front();Q.pop();
        vis[s]=false;
        for(int i=head[s];~i;i=edge[i].next){
            int t=edge[i].to;
            if(judge(s,t,l,r)){
                if(d[t]>d[s]+edge[i].v){
                    d[t]=d[s]+edge[i].v;
                    if(!vis[t]){
                        vis[t]=true;Q.push(t);
                    }
                }
            }
        }
    }
    return d[m];
}
void cal(){
    for(int i=1;i<=n;i++){
        for(int j=i;j<=n;j++){
            cost[i][j]=spfa(i,j);
        }
    }
    dp[0]=0;
    for(int i=1;i<=n;i++){
        dp[i]=cost[1][i]*i;
        for(int j=1;j<=i;j++){
            ll c=cost[j+1][i];
            dp[i]=min(dp[i],dp[j]+c*(i-j)+k);
        }
    }
    printf("%lld\n",dp[n]);
}
int main()
{
    //freopen("/home/xt/code/acm/input.txt","r",stdin);
    scanf("%d%d%d%d",&n,&m,&k,&e);
    Ecnt=0;memset(head,-1,sizeof head);
    for(int i=0;i<e;i++){
        int a,b,c;scanf("%d%d%d",&a,&b,&c);
        addedge(a,b,c);
    }
    scanf("%d",&q);//数组开大点
    for(int i=0;i<q;i++){
        scanf("%d%d%d",&block[i].id,&block[i].l,&block[i].r);
    }
    cal();
    //printf("[Run in %.1fs]\n",(double)clock()/CLOCKS_PER_SEC);
    return 0;
}