codevs1403 新三国争霸

题目描述 Description

PP 特别喜欢玩即时战略类游戏,但他觉得那些游戏都有美中不足的地方。灾害总不降临道路,而只降临城市,而且道路不能被占领,没有保护粮草的真实性。于是他就研发了《新三国争霸》。
在这款游戏中,加入灾害对道路的影响(也就是一旦道路W[i,j]受到了灾害的影响,那么在一定时间内,这条路将不能通过)和道路的占领权(对于一条道路W[i,j],至少需要K[i,j]个士兵才能守住)。
PP可真是高手,不一会,就攻下了N-1座城市,加上原来的就有N座城市了,但他忽略了一点……那就是防守同样重要,不过现在还来的及。因为才打完仗所以很多城市都需要建设,PP估算了一下,大概需要T天。他现在无暇分身进攻了,只好在这T天内好好的搞建设了。所以他秒要派士兵占领一些道路,以确保任何两个城市之间都有路(不然敌人就要分而攻之了,是很危险的)。士兵可不是白干活的,每个士兵每天都要吃掉V的军粮。因为有灾害,所以方案可能有变化(每改变一次就需要K的军粮,初始方案也需要K的军粮)。
因为游戏是PP编的,所以他知道什么时候有灾害。PP可是一个很节约的人,他希望这T天在道路的防守上花最少的军粮。
N<=300,M<=5000 ,T<=50;

输入描述 Input Description

第一行有5个整数N,M,T,V,K。N表示有城市数,M表示道路数,T表示需要修养的天数,V表示每个士兵每天吃掉的军粮数,K表示修改一次花掉的军粮数。
以下M行,每行3个数A,B,C。表示A与B有一条路(路是双向的)需要C个士兵才能守住。
第M+2行是一个数P,表示有P个灾害。
以下P行,每行4个数,X,Y,T1,T2。表示X到Y的这条路,在T1到T2这几天都会受灾害。

输出描述 Output Description

T天在道路的防守上花费最少的军粮。

样例输入 Sample Input

3 3 5 10 30
1 2 1
2 3 2
1 3 4
1
1 3 2 5

样例输出 Sample Output

180

数据范围及提示 Data Size & Hint

各个测试点1s

/*
这道题和物流运输这道题非常像,还是枚举i~j天的最小生成树,然后 搞一搞区间dp,一开始并查集都打错了- -!智障
*/
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define ll long long
#define fo(i,l,r) for(int i = l;i <= r;i++)
#define fd(i,l,r) for(int i = r;i >= l;i--)
using namespace std;
const int N = 405,M=12005;
ll read(){
    ll x=0,f=1;
    char ch=getchar();
    while(!(ch>='0'&&ch<='9')){if(ch=='-')f=-1;ch=getchar();};
    while(ch>='0'&&ch<='9'){x=x*10+(ch-'0');ch=getchar();};
    return x*f;
}
struct edge{
    int u;
    int v;
    ll w;
    friend bool operator < (edge a,edge b){
        return a.w < b.w;
    }
}e[M];
ll n,m,t,k,val,p,sumv;
ll rec[55][55],dp[55];
bool zh[305][305][55],vis[M];
int cnt,fa[N],tot;
int findf(int x){
    return x == fa[x] ? fa[x] : fa[x] = findf(fa[x]);
}
void input(){
    n=read();m=read();t=read();val=read();k=read();
    fo(i,1,m){
        e[i].u = read();
        e[i].v = read();
        if(e[i].u > e[i].v) swap(e[i].u,e[i].v);
        e[i].w = read();
        sumv += e[i].w;
    }
    p =read();
    int u,v,t1,t2;
    fo(i,1,p){
        u = read();v = read();t1 = read();t2 = read();
        if(u > v) swap(u,v);
        if(t1 > t2) swap(t1,t2);
        fo(j,t1,t2) zh[u][v][j] = true;
    }
}
void mst(int lp,int rp){
    int u,v;
    int cnt_n = 0;
    fo(i,1,m){
        fo(j,lp,rp){
            if(zh[e[i].u][e[i].v][j]){
                vis[i] = true;
                break;
            }
        }
    }
    fo(i,1,m){
        u = findf(e[i].u);
        v = findf(e[i].v);
        if(fa[u] == v || vis[i]) continue;
        fa[u] = v;
        tot += e[i].w * val;
        cnt_n++;
        if(cnt_n == n-1) break;
    }
    if(cnt_n != n-1) rec[lp][rp] = 987654321012345LL;
    else rec[lp][rp] = tot*(rp-lp+1);
}
void get_rec(){
    sort(e+1,e+1+m);
    fo(i,1,t){
        fo(k,i,t){
            fo(j,1,N) fa[j] = j;
            fo(j,1,m) vis[j] = false;
            tot = 0;
            mst(i,k);
        }    
    }
}
void get_ans(){
    fo(i,1,51) dp[i] = 987654321012345LL;
    fo(i,1,t){
        fd(j,1,i){
            dp[i] = min(dp[i],rec[j][i] + dp[j-1] + k);
        }
    }
    cout<<dp[t]<<endl;
}
int main(){
    input();
    get_rec();
    get_ans();
    return 0;
} 

 

posted @ 2016-11-06 15:19  ACforever  阅读(148)  评论(0编辑  收藏  举报