黄色网页视频 I 影音先锋日日狠狠久久 I 秋霞午夜毛片 I 秋霞一二三区 I 国产成人片无码视频 I 国产 精品 自在自线 I av免费观看网站 I 日本精品久久久久中文字幕5 I 91看视频 I 看全色黄大色黄女片18 I 精品不卡一区 I 亚洲最新精品 I 欧美 激情 在线 I 人妻少妇精品久久 I 国产99视频精品免费专区 I 欧美影院 I 欧美精品在欧美一区二区少妇 I av大片网站 I 国产精品黄色片 I 888久久 I 狠狠干最新 I 看看黄色一级片 I 黄色精品久久 I 三级av在线 I 69色综合 I 国产日韩欧美91 I 亚洲精品偷拍 I 激情小说亚洲图片 I 久久国产视频精品 I 国产综合精品一区二区三区 I 色婷婷国产 I 最新成人av在线 I 国产私拍精品 I 日韩成人影音 I 日日夜夜天天综合

bzoj 1096: [ZJOI2007]倉庫建設(shè)

系統(tǒng) 2237 0

? dp是很好想的了,關(guān)鍵是數(shù)據(jù)太大,普通dp肯定超時,所以一定有用某種優(yōu)化,dp優(yōu)化也就那么幾種,這道題用的是 斜率優(yōu)化 ,先寫出普通的狀態(tài)轉(zhuǎn)移方程: dp[i] = min{ ?dp[j] + Σ ( p[k] * (x[i] - x[k] ) )?, ?j+1 <=k <= i , ?0 <= j <= i-1}

這個式子應(yīng)該是很好理解的。接下來,就要進行優(yōu)化。dp[j] 無法改變, 所以只好放眼于第二項, 即sigma那一項

Σ ( p[k] * (x[i] - x[k] ) =?Σ (p[k] * x[i] - p[k] * x[k]) = p[ j+1 ~ i] * x[i] - p[ j+1~i] * x[ j+1 ~i]

我們發(fā)現(xiàn),這個式子中,x[i] 為當(dāng)前點的量,而p[j+1 ~i] 和 p[j+1 ~i] * x[j+1 ~i] 很容易預(yù)處理得到。

于是,我們把 a[i] 定義為 到 i 為止所有貨物的個數(shù) 即 sum( p[1~i] ) ; 把b[i] 定義為到 i 為止 所有 p[j] * x[j] 之和

即 sum( p[ j+1~i] * x[ j+1 ~i]?) ;

我們有了這兩個預(yù)處理,就可以轉(zhuǎn)化成斜率優(yōu)化來做了,一開始的式子,我們轉(zhuǎn)化為

dp[i] = min { dp[j] + x[i] * (a[i] - a[j]) - (b[i] - b[j]) }

剩下的就是斜率優(yōu)化的內(nèi)容了,這里不再贅述,注意一點, a[i] - a[j] 是負(fù)數(shù) 除過去要變號。

代碼如下

      #include <cstdio>
      
        

#include 
      
      <cstring>
      
        

#include 
      
      <cstdlib>
      
        

#include 
      
      <iostream>
      
        

#include 
      
      <algorithm>
      
        

#include 
      
      <queue>


      
        #define
      
       N 1000100


      
        using
      
      
        namespace
      
      
         std;

 


      
      
        int
      
      
         n;

deque
      
      <
      
        int
      
      >
      
         q;

deque
      
      <
      
        int
      
      >
      
        ::iterator x, y;


      
      
        long
      
      
        long
      
       a[N]={
      
        0
      
      }, b[N]={
      
        0
      
      
        };


      
      
        long
      
      
        long
      
      
         f[N], dis[N], c[N];

 


      
      
        long
      
      
        long
      
       getup(
      
        int
      
       j, 
      
        int
      
       k) { 
      
        return
      
       f[j] - f[k] + b[j] -
      
         b[k]; }


      
      
        long
      
      
        long
      
       getdown(
      
        int
      
       j, 
      
        int
      
       k) { 
      
        return
      
       a[j] -
      
         a[k]; }


      
      
        long
      
      
        long
      
       getans(
      
        int
      
       j, 
      
        int
      
       now) { 
      
        return
      
       f[j] + (a[now] - a[j]) * dis[now] -b[now] + b[j] +
      
        c[now]; }


      
      
        bool
      
       ketan(
      
        int
      
       j, 
      
        int
      
       k, 
      
        int
      
       now) { 
      
        return
      
       getup(j, k) < dis[now] *
      
         getdown(j, k); }


      
      
        bool
      
       keya(
      
        int
      
       i, 
      
        int
      
       j, 
      
        int
      
       k) { 
      
        return
      
       getup(i, j) * getdown (j, k) > getup(j, k) *
      
         getdown(i, j); }

 


      
      
        int
      
      
         main()

{

    scanf(
      
      
        "
      
      
        %d
      
      
        "
      
      ,&
      
        n);

    
      
      
        for
      
       (
      
        int
      
       i = 
      
        1
      
      ; i <= n; ++
      
        i)

    {

        
      
      
        long
      
      
        long
      
      
         z;

        scanf(
      
      
        "
      
      
        %lld%lld%lld
      
      
        "
      
      , &dis[i], &z, &
      
        c[i]);

        a[i] 
      
      = a[i-
      
        1
      
      ] + z; b[i] = b[i-
      
        1
      
      ] + z*
      
        dis[i];

    }

    f[
      
      
        0
      
      ] = 
      
        0
      
      ; q.push_front(
      
        0
      
      
        );

    
      
      
        for
      
       (
      
        int
      
       i = 
      
        1
      
      ; i <= n; ++
      
        i)

    {

        
      
      
        while
      
       (q.size() > 
      
        1
      
      
        )

        {

            x 
      
      = q.begin(); y = x; y++
      
        ;

            
      
      
        if
      
       (!ketan(*y, *x, i)) 
      
        break
      
      
        ;

            q.pop_front();

        }

        x 
      
      =
      
         q.begin();

        f[i] 
      
      = getans(*
      
        x, i);

        
      
      
        while
      
       (q.size() > 
      
        1
      
      
        )

        {

            x 
      
      = q.end(); x--; y = x; y--
      
        ;

            
      
      
        if
      
       (!keya(*y, *x, i)) 
      
        break
      
      
        ;

            q.pop_back();

        }

        q.push_back(i);

    }

    printf(
      
      
        "
      
      
        %lld\n
      
      
        "
      
      
        , f[n]);

}
      
    

?

bzoj 1096: [ZJOI2007]倉庫建設(shè)


更多文章、技術(shù)交流、商務(wù)合作、聯(lián)系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯(lián)系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦?。?!

發(fā)表我的評論
最新評論 總共0條評論