首页 > 代码库 > 【CEOI2004】锯木厂选址

【CEOI2004】锯木厂选址

【题目描述】

从山顶上到山底下沿着一条直线种植了n棵老树。当地的政府决定把他们砍下来。为了不浪费任何一棵木材,树被砍倒后要运送到锯木厂。
木材只能按照一个方向运输:朝山下运。山脚下有一个锯木厂。另外两个锯木厂将新修建在山路上。你必须决定在哪里修建两个锯木厂,使得传输的费用总和最小。假定运输每公斤木材每米需要一分钱。

【输入描述】

输入的第一行为一个正整数n——树的个数(2≤n≤20 000)。树从山顶到山脚按照1,2……n标号。接下来n行,每行有两个正整数(用空格分开)。第i+1行含有:wi——第i棵树的重量(公斤为单位)和 di——第i棵树和第i+1棵树之间的距离,1≤wi ≤10 000,0≤di≤10 000。最后一个数dn,表示第n棵树到山脚的锯木厂的距离。保证所有树运到山脚的锯木厂所需要的费用小于2000 000 000分。

【输出描述】

输出只有一行一个数:最小的运输费用。

【分析】

斜率优化,看别人的,只提供代码,以后慢慢总结。(觉得自己实在是太水了==)

贴上链接:

http://wenku.baidu.com/linkurl=ptHSVhAevtr4BAOqM6U8MOzdJjFd_bjtEVcm34UULAEEnXXrltD6qERloCsy1CwgnjDY65pffDK3f7xWzvIztYntVK7z5c8iGinaP2EBaK7

 

 1 #include <cstdlib>
 2 #include <iostream>
 3 #include <cstring>
 4 #include <cstdio>
 5 #include <cmath>
 6 #include <algorithm>
 7 const int maxn=30005;
 8 const int INF=2000000005;
 9 using namespace std;
10 int sw[maxn],sd[maxn],cost[maxn],Q[2*maxn];
11 inline double g(int i,int j) {return ((sw[i]*sd[i]-sw[j]*sd[j]*1.0)/(sw[i]-sw[j]));}
12 inline int all(int i,int j) {return cost[j]-cost[i-1]-sw[i-1]*(sd[j]-sd[i-1]);}
13 int main()
14 {
15     int n,i,front=1,rear=1,ans=INF,w,d;
16     //文件操作
17     freopen("two.in","r",stdin);
18     freopen("two.out","w",stdout);
19     memset(sw,0,sizeof(sw));
20     memset(sd,0,sizeof(sd));
21     memset(cost,0,sizeof(cost));
22     scanf("%d",&n);
23     //从山顶往下输入 
24     for (i=1;i<=n;i++)
25     {
26         scanf("%d%d",&w,&d);
27         sw[i]=sw[i-1]+w;
28         sd[i+1]=sd[i]+d;
29         cost[i+1]=cost[i]+sw[i]*d;
30     }Q[1]=1;//初始化队列 
31     for (i=2;i<=n;i++)
32     {
33         //满足队列单调性 
34         while (front<rear && g(Q[front],Q[front+1])<=sd[i]) front++;
35         ans=min(ans,cost[Q[front]]+all(Q[front]+1,i)+all(i+1,n+1));
36         while (front<rear && g(Q[rear],i)<g(Q[rear-1],Q[rear]))rear--;//弹出队列 
37         Q[++rear]=i;//加入队列 
38     }
39     printf("%d\n",ans);
40     return 0;
41 }
View Code