BZOJ1003

1003: [ZJOI2006]物流运输trans

思路：最短路+DP
每一天都走最短路显然是不现实的
不过我们发现天数很小
设cost[i][j]——i到j天内都可以走的最短路距离(若不存在则为inf)
然后我们在进行一个DP
设f[i]——前i天运输的最小费用
显然，状态转移方程为：
f[i]=min(f[j]+cost[j+1][i]*(i-j)+K) (1<=j< i)
为了方便，初始状态为f[0]=-K
感觉做了好多事啊！

代码：
这应该是我改代码风格前的作品
难怪dijkstra都打错了……

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <iostream>
#define MIN(a,b) ((a)<(b)?(a):(b))
using namespace std;
const int MAXDAY=110,MAXN=25,INF=1000000;
struct Front
{
    int x,y;
    int num;
}front[2*MAXN*MAXN];
int orderx[2*MAXN*MAXN],stdx[MAXN],num[MAXN];
bool use[MAXDAY][MAXN];
bool check[MAXN];
int cost[MAXDAY][MAXDAY];
int f[MAXDAY];
int day,n,K,e,m,d,maxroad;
inline void add(int x,int y,int num)
{
    front[++m].x=x;
    front[m].y=y;
    front[m].num=num;
    return ;
}
inline void sortx()
{
    for (int i=1;i<=m;i++)
        num[front[i].x]++;
    for (int i=1;i<=n;i++)
        {
        num[i]+=num[i-1];
        stdx[i]=num[i];
        }
    for (int i=1;i<=m;i++)
        orderx[++num[front[i].x-1]]=i;
    return ;
}
inline int dijstra(int l,int r)
{
    int ans[MAXN];
    memset(check,false,sizeof(check));
    for (int i=0;i<=n;i++)
        ans[i]=INF;
    ans[1]=0;
    for (int i=1;i<n;i++)
        {
        int expend=0;
        for (int j=1;j<=n;j++)
            if (!check[j] && ans[j]<ans[expend])
                expend=j;
        check[expend]=true;
        for (int j=stdx[expend-1]+1;j<=stdx[expend];j++)
            {
            int next=front[orderx[j]].y,COST=front[orderx[j]].num;
            bool pass=true;
            for (int k=l;k<=r;k++)
                if (!use[k][next])
                    {
                    pass=false;
                    continue;
                    }
            if (!pass)
                continue;
            ans[next]=MIN(ans[expend]+COST,ans[next]);
            }
        }
    return ans[n];
}
inline void init()
{
    cin>>day>>n>>K>>e;
    for (int i=1;i<=e;i++)
        {
        int x,y,num;
        scanf("%d%d%d",&x,&y,&num);
        add(x,y,num);
        add(y,x,num);
        }
    memset(use,true,sizeof(use));
    sortx();
    int d;
    cin>>d;
    while (d--)
        {
        int x,l,r;
        scanf("%d%d%d",&x,&l,&r);
        for (int i=l;i<=r;i++)
            use[i][x]=false;
        }
    for (int i=1;i<=day;i++)
        for (int j=i;j<=day;j++)
            cost[i][j]=dijstra(i,j);
    for (int i=1;i<=day;i++)
        f[i]=INF;
    f[0]=-K;
    for (int i=1;i<=day;i++)
        for (int j=0;j<i;j++)
            f[i]=MIN(f[i],f[j]+cost[j+1][i]*(i-j)+K);
    cout<<f[day];
    return ;
}
int main()
{
    init();
    return 0;
}