codevs 1403 新三国争霸

用dp[i]表示第i天的最少花费
cost[i][j]表示从i到j天方案不变时的士兵数
用双层循环+kruskal预处理出cost[i][j]

P[u][v][t]表示u到v的路在t时间是否可用 (即灾害)
按照时间轴进行dp，外层循环i，内层循环j表示它的方案从第j天转移而来
dp[i]=min(dp[i],dp[j]+cost[j+1][i]*V*(i-j));
即从第j天改变方案，选择最小合法方案后，一直到第i天
初始化dp[i]=cost[1][i]*i*V; //注意cost[1][i]可能等于inf
最后dp[i]+=K;表示在某一天改变了一次方案

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#define re register
#define llint long long
#define wl putchar('
')
#define read1(a) scanf("%d",&a)
#define De(a) cout<<a<<endl
using namespace std;
template <typename ty> inline void read(ty &x)
{
    x=0;int f=1;re char c=getchar();
    for(;c<'0'||c>'9';c=getchar()) if(c=='-') f=-1;
    for(;c>='0'&&c<='9';c=getchar()) x=(x<<1)+(x<<3)+(c^48);
    x*=f;
}
template <typename ty> inline void write(ty x)
{
    if(x<0) x=-x,putchar('-');
    if(x>9) write(x/10);
    putchar(x%10+48);
}
const int inf = 0x3f3f3f3f;
struct edge{
    int u,v,c;
}e[5005<<1];
int fa[5000];
int cost[500][500],dp[500];
inline bool cmp(edge a,edge b)
{
    return a.c<b.c;
}
bool P[305][305][60];
int N,M,T,V,K;
int q;

inline int find(int x)
{
    return x==fa[x]?x:fa[x]=find(fa[x]);
}
bool check(int u,int v,int t1,int t2)
{
    for(int i=t1;i<=t2;i++)
        if(P[u][v][i]) return false;
    return true;
}
inline int kruskal(int t1,int t2)
{
    int cnt=1,ans=0;
    for(int i=1;i<=N;i++)
        fa[i]=i;
    
    for(int i=1;i<=M;i++)
    {
        if(e[i].u==e[i].v) continue;
        if(check(e[i].u,e[i].v,t1,t2)==0) continue;
        int fx=find(e[i].u),fy=find(e[i].v);
        
        if(fx==fy) continue;
        
        ans+=e[i].c;
        fa[fy]=fx;
        cnt++;//在这里数是否把所有的点都加入mst 
        if(cnt == N) break;
    }
    if(cnt==N) return ans;
    return inf;
}

int main()
{
    read(N),read(M),read(T),read(V),read(K);
    for(re int i=1;i<=M;i++)
    {
        read(e[i].u),read(e[i].v),read(e[i].c);
    }
    sort(e+1,e+M+1,cmp);
    read(q);
    for(re int i=1;i<=q;i++)
    {
        int u,v,t1,t2;
        read(u),read(v),read(t1),read(t2);
        for(int j=t1;j<=t2;j++)
        {
            P[u][v][j]=P[v][u][j]=1;//uv路径不可用 
        }
    }
    for(int i=1;i<=T;i++)
    {
        for(int j=i;j<=T;j++)
        {
            cost[i][j] = kruskal(i,j);//预处理从i到j天的mst 
            printf("%d ",cost[i][j]);
        }
        puts("");
    }
    for(int i=1;i<=T;i++)
    {
        if(cost[1][i]!=inf) dp[i]=cost[1][i]*i*V;//表示初始化为一直用第一天开始的方案
        else dp[i]=cost[1][i];//防止爆掉
        for(int j=1;j<i;j++)
        {
            if(cost[j+1][i]!=inf)//若从j+1到i天存在mst 
                dp[i]=min(dp[i],dp[j]+cost[j+1][i]*V*(i-j));//cost[j+1][i]*V*(i-j)表示从j+1天到i天的花费 
        }
        dp[i]+=K;//这里加k，以后转移的时候也会包含这个k 
    }
    write(dp[T]);
}