【bzoj1003】[ZJOI2006]物流运输trans 最短路+dp
2016-02-13 16:12
381 查看
非常容易的dp,一开始竟然还卡了一下。
dp[i]表示到第i天的最小花费
dp[i]=min{dp[j]+cost[j+1][i]*(i-j)+k} (0<=j<i)
cost[i][j]表示从第i天到第j天的最短路,预处理即可。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define inf 1000000000
using namespace std;
bool tag[21],a[21][110],vis[21];
int head[21],to[1010],next[1010],len[1010],dis[21],q[21];
int f[110][110],dp[110];
int n,m,num,k,Day,cnt;
void addedge(int x,int y,int z)
{
num++;to[num]=y;len[num]=z;next[num]=head[x];head[x]=num;
}
void spfa()
{
for (int i=2;i<=n;i++) dis[i]=inf;
dis[1]=0;
int l=0,r=1;
q[1]=1;vis[1]=1;
while (l!=r)
{
l++;if (l==21) l=0;
int x=q[l];
for (int p=head[x];p;p=next[p])
if (dis[x]+len[p]<dis[to[p]])
{
dis[to[p]]=dis[x]+len[p];
if (!vis[to[p]] && !tag[to[p]])
{
r++;if (r==21) r=0;
q[r]=to[p];vis[to[p]]=1;
}
}
vis[x]=0;
}
}
int main()
{
scanf("%d%d%d%d",&Day,&n,&k,&m);
for (int i=1;i<=m;i++)
{
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
addedge(x,y,z);addedge(y,x,z);
}
scanf("%d",&cnt);
for (int i=1;i<=cnt;i++)
{
int x,l,r;
scanf("%d%d%d",&x,&l,&r);
for (int j=l;j<=r;j++) a[x][j]=1;
}
for (int i=1;i<=Day;i++)
{
memset(tag,0,sizeof(tag));
for (int j=i;j<=Day;j++)
{
for (int k=1;k<=n;k++) tag[k]|=a[k][j];
spfa();
f[i][j]=dis
;
}
}
dp[0]=-k;
for (int i=1;i<=Day;i++)
{
dp[i]=inf;
for (int j=0;j<i;j++) if (f[j+1][i]!=inf) dp[i]=min(dp[i],dp[j]+f[j+1][i]*(i-j)+k);
}
printf("%d\n",dp[Day]);
return 0;
}
dp[i]表示到第i天的最小花费
dp[i]=min{dp[j]+cost[j+1][i]*(i-j)+k} (0<=j<i)
cost[i][j]表示从第i天到第j天的最短路,预处理即可。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define inf 1000000000
using namespace std;
bool tag[21],a[21][110],vis[21];
int head[21],to[1010],next[1010],len[1010],dis[21],q[21];
int f[110][110],dp[110];
int n,m,num,k,Day,cnt;
void addedge(int x,int y,int z)
{
num++;to[num]=y;len[num]=z;next[num]=head[x];head[x]=num;
}
void spfa()
{
for (int i=2;i<=n;i++) dis[i]=inf;
dis[1]=0;
int l=0,r=1;
q[1]=1;vis[1]=1;
while (l!=r)
{
l++;if (l==21) l=0;
int x=q[l];
for (int p=head[x];p;p=next[p])
if (dis[x]+len[p]<dis[to[p]])
{
dis[to[p]]=dis[x]+len[p];
if (!vis[to[p]] && !tag[to[p]])
{
r++;if (r==21) r=0;
q[r]=to[p];vis[to[p]]=1;
}
}
vis[x]=0;
}
}
int main()
{
scanf("%d%d%d%d",&Day,&n,&k,&m);
for (int i=1;i<=m;i++)
{
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
addedge(x,y,z);addedge(y,x,z);
}
scanf("%d",&cnt);
for (int i=1;i<=cnt;i++)
{
int x,l,r;
scanf("%d%d%d",&x,&l,&r);
for (int j=l;j<=r;j++) a[x][j]=1;
}
for (int i=1;i<=Day;i++)
{
memset(tag,0,sizeof(tag));
for (int j=i;j<=Day;j++)
{
for (int k=1;k<=n;k++) tag[k]|=a[k][j];
spfa();
f[i][j]=dis
;
}
}
dp[0]=-k;
for (int i=1;i<=Day;i++)
{
dp[i]=inf;
for (int j=0;j<i;j++) if (f[j+1][i]!=inf) dp[i]=min(dp[i],dp[j]+f[j+1][i]*(i-j)+k);
}
printf("%d\n",dp[Day]);
return 0;
}
相关文章推荐
- php session 写入数据库
- [iOS]APP代码实践:建立一个辅助的APP类,减少对AppDelegate的修改
- 英语语法最终珍藏版笔记-6“情态动词+have+ done”的含义
- 168. Excel Sheet Column Title LeetCode
- hash(哈希散列) 数据结构-9
- VMware Player 使用错误集锦
- 海马记忆训练
- usaco 1.5.3 sprime
- poj1001 2010.2.25
- hdu2824 2010.2.26
- 坚持读读的英语口语
- pig hive hbase之间的关系和差异
- 函数指针示例
- 2016/2/13 Head First Java 第7节第二小结 多态
- 【Hibernate系列】(四):基本映射
- poj1125 2010.2.22
- 英语语法最终珍藏版笔记-5情态动词
- Java程序设计的DOS命令基础
- 将博客搬至CSDN
- 英语语法最终珍藏版笔记-4助动词和情态动词