CODE[VS] 1403 新三国争霸

题目描述 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

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

DP+最小生成树~

DP就是线性的，主要是用最小生成树预处理出d[i][j]，即从i天到j天不改变路径的最小花费很神奇~

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;int n,m,t,v,k,x,y,p,k1,k2,fa[301],f[51],d[51][51];bool b[51][301][301];struct node{int x,y,v;}a[5001];bool cmp(node u,node v){return u.v<v.v;}int findd(int u){return fa[u]==u ? u:fa[u]=findd(fa[u]);}bool che(int u,int v,int xx,int yy){for(int i=u;i<=v;i++)  if(b[i][xx][yy]) return 0;return 1;}int cal(int u,int v){for(int i=1;i<=n;i++) fa[i]=i;int tot=0,kx,ky;for(int i=1;i<=m;i++)  if((kx=findd(a[i].x))!=(ky=findd(a[i].y)) && che(u,v,a[i].x,a[i].y))  {  fa[kx]=ky;tot+=a[i].v;  }for(int i=2;i<=n;i++)  if(findd(i)!=findd(1)) return 999999999;return tot;}int main(){scanf("%d%d%d%d%d",&n,&m,&t,&v,&k);for(int i=1;i<=m;i++){scanf("%d%d%d",&x,&y,&a[i].v);if(x>y) swap(x,y);a[i].x=x;a[i].y=y;}sort(a+1,a+m+1,cmp);scanf("%d",&p);while(p--){scanf("%d%d%d%d",&x,&y,&k1,&k2);if(x>y) swap(x,y);for(int i=k1;i<=k2;i++) b[i][x][y]=1;}for(int i=1;i<=t;i++)  for(int j=i;j<=t;j++) d[i][j]=cal(i,j);for(int i=1;i<=t;i++){f[i]=k+d[1][i]*v*i;for(int j=1;j<i;j++)  if(d[j+1][i]!=999999999 && f[i]>d[j+1][i]*(i-j)*v+k+f[j])    f[i]=d[j+1][i]*(i-j)*v+k+f[j];}printf("%d\n",f[t]);return 0;}