zoukankan      html  css  js  c++  java
  • BZOJ 1705: [Usaco2007 Nov]Telephone Wire 架设电话线 DP + 优化 + 推导

    Description

    最近,Farmer John的奶牛们越来越不满于牛棚里一塌糊涂的电话服务 于是,她们要求FJ把那些老旧的电话线换成性能更好的新电话线。 新的电话线架设在已有的N(2 <= N <= 100,000)根电话线杆上, 第i根电话线杆的高度为height_i米(1 <= height_i <= 100)。 电话线总是从一根电话线杆的顶端被引到相邻的那根的顶端 如果这两根电话线杆的高度不同,那么FJ就必须为此支付 C*电话线杆高度差(1 <= C <= 100)的费用。当然,你不能移动电话线杆, 只能按原有的顺序在相邻杆间架设电话线。Farmer John认为 加高某些电话线杆能减少架设电话线的总花费,尽管这项工作也需要支出一定的费用。 更准确地,如果他把一根电话线杆加高X米的话,他得为此付出X^2的费用。 请你帮Farmer John计算一下,如果合理地进行这两种工作,他最少要在这个电话线改造工程上花多少钱。

    Input

    * 第1行: 2个用空格隔开的整数:N和C

    * 第2..N+1行: 第i+1行仅有一个整数:height_i

    Output

    * 第1行: 输出Farmer John完成电话线改造工程所需要的最小花费

    题解: 

    首先,加高后的最高高度一定不大于当前最高高度. 

    $f[i][j]$ 表示第 $i$ 个柱子高度为 $j$ 时的最小花费.

    $f[i][j]=(j-h[i])^{2}+min(f[i-1][k]+c imes left | k-j ight |)$

    绝对值符号有些不好处理,再把转移方程拆一下:

    $f[i][j]=(j-h[i])^{2}+min(f[i-1][k]-c imes k)$ ,$k<=j$

    $f[i][j]=(j-h[i])^{2}+min(f[i-1][k]+c imes k)$, $k>=j$

    不难看出,对于 $f[i][j]$ 来说, $(j-h[i])^{2}$ 是固定的,直接算即可.

    我们设 $min1[i][j]$ 为 $k<=j$ 时 $f[i][k]-c imes k$ 的最小值.

    不难得出 : $min1[i][j]=min(min1[i][j-1], f[i][j]-c imes j)$

    而 $f[i][j]=(j-h[i])^{2}+min1[i-1][j]$.

    于是,我们用下标为 $i-1$ 的答案来更新当前下标为 $i$ 的答案,再用当前的 $f[i][j]$ 来更新下标为 $i$ 的 $min1[i][j]$ 即可,每次转移时
    $O(1)$ 的。
      
    #include<bits/stdc++.h>
    #define setIO(s) freopen(s".in","r",stdin) 
    #define Max(a,b) (a=a>b?a:b) 
    #define Min(a,b) (a=a<b?a:b) 
    using namespace std;
    int f[103],h[100003],min1[103],min2[103];
    int main()
    {
    	// setIO("input");  
    	int n,c,maxv=0; 
    	scanf("%d%d",&n,&c);
    	for(int i=1;i<=n;++i) 
    		scanf("%d",&h[i]),Max(maxv,h[i]); 
    	memset(f,0x3f,sizeof(f));  
    	for(int i=h[1];i<=maxv;++i) f[i]=(i-h[1])*(i-h[1]); 
    	for(int i=2;i<=n;++i)
    	{
    		memset(min1,0x3f,sizeof(min1)); 
    		memset(min2,0x3f,sizeof(min2)); 
    		for(int j=h[i-1];j<=maxv;++j) 
    			min1[j]=min(min1[j-1],f[j]-c*j); 
    		for(int j=maxv;j>=1;--j) 
    			min2[j]=min(min2[j+1],j>=h[i-1]?f[j]+c*j:0x3f3f3f3f); 
    		for(int j=h[i];j<=maxv;++j) 
    			f[j]=(j-h[i])*(j-h[i])+min(min1[j]+c*j,min2[j]-c*j); 
    	}
    	int ans=0x3f3f3f3f; 
    	for(int i=h[n];i<=maxv;++i) ans=min(ans,f[i]); 
    	printf("%d
    ",ans); 
    	return 0; 
    }
    

      

  • 相关阅读:
    C#调用C++ ---参数传递
    Retained Mode Versus Immediate Mode
    mind map in latex
    vk example
    基本环境
    [转]ld 和 ld.gold 和 ld.bfd
    [转] c++11 int&& 右值引用
    [转] c++11列表初始化
    [转] c++ const, volatile, mutable用法
    [转] c++11 模板元编程
  • 原文地址:https://www.cnblogs.com/guangheli/p/10974185.html
Copyright © 2011-2022 走看看