青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

算法學社
記錄難忘的征途
posts - 141,comments - 220,trackbacks - 0
題目簡介:
    給你一個序列A,請你把序列A分成連續K個子段,每個子段的代價是 sum(A[i]*A[j]) 其中 i < j。請問如何分組使代價最小。
    數據范圍|A|,K <1000

吐槽:
    cppblog的editor可真是不咋的...  昨天弄了一下午LyX,想用XeTeX編譯,結果屢屢失敗。用vim寫的TeX文檔是可以用XeTeX編譯并正確顯示中文的。但是不知道為什么用LyX寫就不可以.... 有沒有TeX高手來指點一下啊...

解法分析:
    首先狀態轉移方程比較容易想到:DP[i][k]表示前i項分成了k組,那么有DP[i][k] = max(DP[j][k-1]+cost[j+1][i])對所有 j < i;
    這是一個2D/1D的方程,時間復雜度是O(n^3),對于n = 1000來說必然超時,需要優化...
    (A[j]+...+A[i])^2 = A[j]^2 + A[j+1]^2 + ... + A[i]^2 + 2* cost[j][i]; (沒有TeX寫這個就是不爽啊)
    于是乎 cost[j+1][i] = ((sum[i]-sum[j])^2 - suma[i] + suma[j])/2 = -sum[i] * sum[j] + (suma[j] - sum[j]^2)/2 + (sum[i]^2 - suma[i])/2;
    終于把這個轉成只和 i,j相關的表達式了, 可以斜率優化了:
    X(j) = sum[j]
    Y(j) = (suma[j] - sum[j]^2)/2 + dp[j][k-1];
    C(i) = (sum[i]^2 - suma[i])/2
    我們要求的是 dp[i][k] = -sum[j]*X(j) + Y(j) + C(i);
    移項得 Y(j) = X(j)*sum[j] + dp[i][k] - C(i) 轉成了線性規劃問題。 不難想到最優解一定在凸包上。
    而且X,Y是單調的,可以用棧維護凸包。斜率是單調的,可以用隊列維護最優決策。
代碼部分:
#include<iostream>
#include<cstdio>
#include<cassert>
using namespace std;
#define re(i,n) for(int i = 0; i< n; i++)
const int N = 1001;
typedef long long ll;
ll dp[N][N];
ll num[N],pos[N][2],sum[N],suma[N];
int Q[N];
bool cal(int p0, int p1 , ll k){
return pos[p1][1] - pos[p0][1] < k * ( pos[p1][0] - pos[p0][0]);
}
ll cross( int p0, int p1, int p2){
return (pos[p0][0] - pos[p1][0] ) * ( pos[p0][1] - pos[p2][1] ) - ( pos[p0][1] - pos[p1][1]) * ( pos[p0][0] - pos[p2][0] );
}
int main(){
int n,m ;
while(~scanf("%d%d",&n,&m) && n){
re(i,n)
scanf("%lld",&num[i]);
re(i,n) sum[i] = (i == 0 ? 0 : sum[i-1]) + num[i];
re(i,n) dp[i][0] = ( i == 0 ? 0 : dp[i-1][0] + num[i] * sum[i-1]);
re(i,n) suma[i] = (i == 0 ? 0 : suma[i-1])+ num[i] * num[i];
for(int j = 1; j<=m ; j++) {
int head = 0 ,tail =0 ;
for(int i = j; i < n ; i++){
pos[i][0] = sum[i-1]; pos[i][1] = dp[i-1][j-1] + (suma[i-1] + sum[i-1]*sum[i-1])/2;
while(tail - head > 1 && cross(Q[tail-2] , Q[tail-1], i) <= 0) tail --;
Q[tail ++ ] = i;
while(tail - head > 1 && cal(Q[head], Q[head +1 ] , sum[i]) ) head ++;
dp[i][j] = pos[Q[head]][1] - pos[Q[head]][0] * sum[i] + (sum[i]*sum[i] - suma[i])/2;
}
}
cout<< dp[n-1][m] <<endl;
}
}
posted on 2012-04-24 14:51 西月弦 閱讀(960) 評論(3)  編輯 收藏 引用 所屬分類: 解題報告

FeedBack:
# re: hdu 2829 動態規劃+斜率優化
2012-08-23 21:18 | ZeroClock
學習了。思路和代碼都很飄逸...可以加你好友嗎?  回復  更多評論
  
# re: hdu 2829 動態規劃+斜率優化
2012-08-24 12:43 | 西月弦
@ZeroClock
可以哦,我的QQ和姓名在首頁上.
我在CF和TC的ID是hanfei19910905  回復  更多評論
  
# re: hdu 2829 動態規劃+斜率優化[未登錄]
2013-05-04 15:57 | 123
于是乎 cost[j+1][i] = ((sum[i]-sum[j])^2 - suma[i] + suma[j])/2 = -sum[i] * sum[j] + (suma[j] - sum[j]^2)/2 + (sum[i]^2 - suma[i])/2;
這里是不是化錯了。。sum[j]是正的,到了后面就變成負的了。這題我也差不多是類似做法,不過我沒有把cost化開來。不知道哪里有問題,所以便來google了。。  回復  更多評論
  
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            亚洲网站在线观看| 一区二区三区四区国产| 久久国产精品一区二区三区四区 | 欧美精品1区| 国语自产精品视频在线看一大j8 | 亚洲成色精品| 中日韩美女免费视频网址在线观看| 久久国产视频网| 国产精品一区二区久久国产| 在线亚洲电影| 亚洲乱码精品一二三四区日韩在线| 久久中文字幕一区| 亚洲二区视频在线| 嫩草影视亚洲| 欧美www视频| 亚洲人午夜精品| 亚洲国产精品va在线看黑人动漫 | 欧美精品一区二区三区在线看午夜 | 欧美日韩小视频| 亚洲资源av| 亚洲欧美三级在线| 国产一区二区黄色| 久久一区二区三区四区五区| 久久久久久综合| 亚洲经典视频在线观看| 亚洲人体大胆视频| 国产精品乱看| 久久综合狠狠综合久久综合88| 久久久久国产精品人| 亚洲国产黄色| 99re66热这里只有精品4| 国产精品yjizz| 久久精品一区二区三区四区| 久久国产欧美| 一本一本久久a久久精品综合麻豆| 日韩一区二区免费看| 国产精品一区二区在线| 久久亚洲欧洲| 欧美日韩国产成人在线观看| 香蕉久久国产| 免费一级欧美片在线观看| 亚洲性视频网址| 久久精品一区二区三区不卡| 最新日韩在线| 亚洲一级片在线观看| 在线观看国产精品网站| 99riav久久精品riav| 国产美女精品在线| 欧美激情麻豆| 国产精品午夜在线观看| 欧美激情视频一区二区三区在线播放| 欧美日韩精品免费观看视一区二区 | 国产精品午夜视频| 国产日本精品| 麻豆精品视频在线观看| 欧美日韩一区在线播放| 久久五月天婷婷| 欧美午夜a级限制福利片| 久久女同互慰一区二区三区| 欧美日韩精品免费观看| 欧美成人激情视频| 国产精品一区二区久久精品| 亚洲电影av在线| 国产日本亚洲高清| 亚洲精品中文字幕女同| 精品成人免费| 午夜精品久久久久久久99黑人| 亚洲毛片视频| 久久天堂国产精品| 久久激情五月激情| 国产精品国产三级欧美二区 | 欧美激情一区三区| 久久这里只有精品视频首页| 国产精品高潮呻吟久久av黑人| 亚洲国产欧美国产综合一区 | 久久久久久久国产| 国产精品久久久久久久第一福利| 亚洲丰满少妇videoshd| 伊人久久综合| 久久精品中文字幕一区| 欧美中文在线视频| 国产精品亚洲精品| 亚洲网站啪啪| 亚洲欧美视频在线观看视频| 欧美福利专区| 亚洲区欧美区| 日韩视频在线观看国产| 欧美/亚洲一区| 亚洲东热激情| 一本久道久久综合狠狠爱| 免费一区二区三区| 亚洲国产精品久久久久秋霞影院| 在线成人欧美| 久久亚洲精品一区| 欧美成人久久| 亚洲精品一区二区三区四区高清| 免费在线欧美视频| 亚洲黄色天堂| 一本色道综合亚洲| 国产精品美女www爽爽爽视频| 亚洲图片在区色| 久久精品国产99| 在线观看欧美日韩国产| 免费不卡在线观看| 亚洲欧洲精品一区二区三区波多野1战4 | 亚洲精品免费在线播放| 欧美电影电视剧在线观看| 亚洲欧洲精品一区二区| 欧美成年人在线观看| 亚洲第一主播视频| 中文国产一区| 国产精品毛片a∨一区二区三区| 亚洲欧美一区二区三区极速播放 | 国产一区二区三区免费观看 | 欧美视频一二三区| av不卡在线| 欧美日韩久久精品| 亚洲一区二区精品| 久久精品中文字幕一区二区三区 | 久热精品视频在线观看一区| 亚洲高清123| 欧美日韩第一区| 亚洲一区二区高清视频| 老司机精品久久| 一区二区三区四区五区视频| 国产日韩欧美日韩大片| 久久久久久综合网天天| 日韩视频一区二区在线观看| 亚洲免费人成在线视频观看| 国内精品一区二区三区| 欧美韩日一区| 性感少妇一区| 亚洲精品免费网站| 香蕉乱码成人久久天堂爱免费 | 欧美成人精品高清在线播放| 一区二区三区 在线观看视| 久久狠狠一本精品综合网| 亚洲精品色婷婷福利天堂| 国产精品入口夜色视频大尺度| 久久久免费av| 亚洲永久精品国产| 欧美激情精品久久久久久变态| 亚洲欧美中文在线视频| 91久久精品久久国产性色也91| 国产精品美女在线观看| 欧美高清视频www夜色资源网| 欧美一区二区观看视频| 99国产精品99久久久久久粉嫩| 欧美成人精品h版在线观看| 亚洲欧美日韩一区在线| 99re这里只有精品6| 国内精品模特av私拍在线观看| 日韩视频不卡中文| 亚洲国产经典视频| 久久九九国产| 亚洲男人的天堂在线观看| 日韩视频第一页| 在线免费日韩片| 国产亚洲激情在线| 国产乱子伦一区二区三区国色天香| 欧美国产乱视频| 久久亚洲国产成人| 久久狠狠一本精品综合网| 亚洲永久免费精品| 亚洲一区亚洲二区| 中日韩男男gay无套| 99国内精品| 99国产精品久久久久久久久久| 亚洲激情另类| 欧美激情一区二区三区全黄| 久久久噜久噜久久综合| 欧美一区视频在线| 欧美一区二区三区在线| 午夜精品久久久久久久久久久| 亚洲午夜精品福利| 中文在线资源观看视频网站免费不卡| 亚洲精品日韩欧美| 亚洲免费观看高清完整版在线观看| 亚洲国产精品t66y| 国产日韩欧美在线一区| 久久精品一区二区国产| 欧美日韩国语| 免费试看一区| 国产欧美日韩一区二区三区在线观看 | 久久人人爽人人爽| 久久精品女人天堂| 久久婷婷国产综合精品青草| 久久久久久久久久久久久9999| 久久久不卡网国产精品一区| 欧美在线视频一区二区三区| 久久er精品视频| 看片网站欧美日韩| 亚洲国产精品福利| 亚洲激情电影在线| 亚洲麻豆国产自偷在线| 亚洲视频1区| 欧美一区亚洲二区| 久久久综合网| 亚洲黄色在线看| 一区二区三区|亚洲午夜|