BZOJ 1597: [Usaco2008 Mar]土地购买动态规划 + 斜率优化

最新推荐文章于 2019-08-05 15:01:00 发布

EM-LGH

最新推荐文章于 2019-08-05 15:01:00 发布

阅读量67

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/liyong1009s/article/details/94802011

本文深入探讨了使用 C++ 实现的一种复杂动态规划算法，通过对比和斜率优化，有效地解决了特定类型的问题。文章详细介绍了算法的实现过程，包括数据结构的设计和关键函数的定义，如比较函数和斜率计算函数。通过对输入数据的排序和动态调整，算法能够高效地找到最优解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Code:

#include<bits/stdc++.h> 
#define maxn 1000000 
#define ll long long 
#define x(i) (b[i+1]) 
#define y(i) (f[i]) 
#define setIO(s) freopen(s".in","r",stdin) 
using namespace std;  
int A[maxn],arr[maxn],q[maxn],tot,head,tail;  
ll a[maxn], b[maxn], f[maxn], u[maxn], v[maxn];         
bool cmp(int i,int j) 
{
    return u[i]==u[j]?v[i]<v[j]:u[i]<u[j]; 
}
double slope(int i,int j)
{
    return (double)(1.00*y(i)-y(j))/(double)(1.00*x(j)-x(i));       
}
int main()
{
    // setIO("input"); 
    int n,i,j; 
    scanf("%d",&n);   
    for(i=1;i<=n;++i) scanf("%lld%lld",&u[i],&v[i]),A[i]=i; 
    sort(A+1,A+1+n,cmp);   
    for(i=1;i<=n;++i) 
    { 
        while(v[A[i]]>=b[tot]&&tot) --tot; 
        ++tot, a[tot]=u[A[i]], b[tot]=v[A[i]];  
    }
    n=tot,head=tail=0;  
    for(i=1;i<=n;++i) 
    {           
        while(head<tail&&slope(q[head],q[head+1])<a[i]) ++head; 
        f[i]=f[q[head]]+a[i]*b[q[head]+1];               
        while(head<tail&&slope(q[tail], i)<slope(q[tail-1],i)) --tail; 
        q[++tail]=i; 
    }
    printf("%lld\n",f[n]);     
    return 0; 
}