BZOJ1003: [ZJOI2006]物流运输

Portal

数据中n,m非常小，所以可以预处理出dis[i][]]表示从第i天到第j天走同一条路线的最短距离。然后dp就行了。

$$f[i]=min(f[j]+K+dis[j+1][i]*(i-j))$$

【代码】

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#define N 25
#define M 805
#define INF 1000000001
using namespace std;
typedef long long ll;

int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n,m,K,E,D,cnt;
int b[M],p[N],nextedge[M],w[M];
int Dis[N],dis[105][105];
bool Flag[N],mark[N];
int c[N][105],sum[N][105],f[N];

void Add(int x,int y,int z){
    cnt++;
    b[cnt]=y;
    nextedge[cnt]=p[x];
    p[x]=cnt;
    w[cnt]=z;
}

void Anode(int x,int y,int z){
    Add(x,y,z);Add(y,x,z);
}

void Spfa()
{
    queue<int>q;
    q.push(1);
    for(int i=2;i<=n;i++) Dis[i]=INF,Flag[i]=0;
    while(!q.empty())
    {
        int k=q.front();q.pop();Flag[k]=0;
        for(int i=p[k];i;i=nextedge[i])
        {
            int v=b[i];
            if(mark[v]) continue;
            if(Dis[v]>Dis[k]+w[i])
            {
                Dis[v]=Dis[k]+w[i];
                if(!Flag[v]) Flag[v]=1,q.push(v);
            }
        }
    }
}

void Input_Init()
{
    m=read(),n=read();K=read(),E=read();
    for(int i=1;i<=E;i++)
    {
        static int x,y,z;
        x=read(),y=read(),z=read();
        Anode(x,y,z);
    }
    D=read();
    for(int i=1;i<=D;i++) 
    {
        static int x;
        x=read();c[x][read()]++,c[x][read()+1]--;
    }
    for(int i=1;i<=n;i++)
    {
        int x=0;
        for(int j=1;j<=m;j++)
            x+=c[i][j],sum[i][j]=sum[i][j-1]+(x!=0);
    }
}

void Pre(int st,int ed)
{
    for(int i=1;i<=n;i++)
        mark[i]=(sum[i][ed]-sum[i][st-1]==0)?0:1;
}

void Work()
{
    for(int i=1;i<=m;i++)
    for(int j=i;j<=m;j++)
    {
        Pre(i,j);
        Spfa();dis[i][j]=Dis[n];
    }
}

void DP()
{
    for(int i=1;i<=m;i++)
    {
        f[i]=INF;if(dis[1][i]!=INF) f[i]=dis[1][i]*i;
        for(int j=2;j<i;j++) if(dis[j+1][i]!=INF)
            f[i]=min(f[i],f[j]+K+dis[j+1][i]*(i-j));
    }
    printf("%d\n",f[m]);
}

int main()
{
    Input_Init();
    Work();
    DP();
    return 0;
}