• 博客园logo
  • 会员
  • 众包
  • 新闻
  • 博问
  • 闪存
  • 赞助商
  • HarmonyOS
  • Chat2DB
    • 搜索
      所有博客
    • 搜索
      当前博客
  • 写随笔 我的博客 短消息 简洁模式
    用户头像
    我的博客 我的园子 账号设置 会员中心 简洁模式 ... 退出登录
    注册 登录
james1207

博客园    首页    新随笔    联系   管理    订阅  订阅

poj 1180 斜率优化dp

这个题目要是顺着dp的话很难做,但是倒着推就很容易退出比较简单的关系式了。

dp[i]=min(dp[u]+(sum[u-1]-sum[i-1]+s)*f[i]);dp[i]代表从i到结尾需要花费的代价,sum[i]表示1到i的时间和,f[i]代表i到n的代价和。

然后对于i状态来说,j由于k等价于 (dp[j]-dp[k])/(sum[k-1]-sum[j-1])<f[i] 

然后f[i]随着i递减而递增,所以就可以利用斜率优化的办法来搞了。

 

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=1e4+9;
long long sum[maxn],f[maxn];
long long dp[maxn];
int que[maxn];

bool chk1(int k,int j,int i)
{
    return (dp[j]-dp[k])<f[i]*(sum[k-1]-sum[j-1]);
}

bool chk2(int k,int j,int i)
{
    return ((dp[j]-dp[k])*(sum[j-1]-sum[i-1]))>((dp[i]-dp[j])*(sum[k-1]-sum[j-1]));
}

int main()
{
    int n,s;
    while(scanf("%d %d",&n,&s)!=EOF)
    {
        sum[0]=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%lld %lld",&sum[i],&f[i]);
            sum[i]+=sum[i-1];
        }
        for(int i=n-1;i>=1;i--) f[i]+=f[i+1];

        int front=1,end=0;
        que[++end]=n+1;
        dp[n+1]=0;
        for(int i=n;i>=1;i--)
        {
            while(front<end&&chk1(que[front],que[front+1],i))
            front++;
            int u=que[front];
            dp[i]=dp[u]+(sum[u-1]-sum[i-1]+s)*f[i];
            while(front<end&&chk2(que[end-1],que[end],i))
            end--;
            que[++end]=i;
        }
        printf("%lld\n",dp[1]);
    }
    return 0;
}


 


 

posted @ 2013-08-16 18:22  Class Xman  阅读(173)  评论(0)    收藏  举报
刷新页面返回顶部
博客园  ©  2004-2025
浙公网安备 33010602011771号 浙ICP备2021040463号-3