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;
}