【BZOJ 1096】[ZJOI2007]仓库建设
【鏈接】 鏈接
【題意】
在這里輸入題意
【題解】
設(shè)f[i]表示在第i個(gè)地方設(shè)立一個(gè)倉(cāng)庫(kù),且前面符合要求的最小花費(fèi)。
則
\(f[i] = min(f[j] + c[i] + dis[i]*(sump[i] - sump[j])-(sumdp[i]-sumdp[j]));\)
其中
sump[]是p[]的前綴和,即\(sump[i] = p[1]+p[2]+..+p[i]\)
sumdp[]是x[]*p[]的前綴和,即\(sumdp[i] = x[1]*p[1]+x[2]*p[2]+...+x[i]*p[i]\)
這個(gè)表示j+1..i這一段的所有物品都運(yùn)輸?shù)絠這個(gè)倉(cāng)庫(kù)。
自己用筆寫(xiě)一下就知道了。
顯然是個(gè)N^2算法。
斜率優(yōu)化。
設(shè)x<y,且y優(yōu)于x.
balabala就能得到
\(\frac{f[y]+sumdp[y]-(f[x]+sumdp[x])}{sump[y]-sump[x]}<dis[i]\)
又dis是單調(diào)遞增的。
這是一個(gè)經(jīng)典的斜率優(yōu)化問(wèn)題了。
【錯(cuò)的次數(shù)】
在這里輸入錯(cuò)的次數(shù)
【反思】
在這里輸入反思
【代碼】
#include <bits/stdc++.h> #define ll long long using namespace std;const int N = 1e6;int n,dl[N+10],h,t; ll dis[N+10],p[N+10],c[N+10],sump[N+10],sumdp[N+10]; ll dp[N+10];double ju(int x,int y) {double fenzi = dp[y] + sumdp[y] - (dp[x]+sumdp[x]);double fenmu = sump[y] - sump[x];return fenzi/fenmu; }int main() {//freopen("F:\\rush.txt","r",stdin);scanf("%d",&n);for (int i = 1;i <= n;i++){scanf("%lld%lld%lld",&dis[i],&p[i],&c[i]);sump[i] = sump[i-1] + p[i];sumdp[i] = sumdp[i-1] + dis[i]*p[i];}for (int i = 1;i <= N;i++) dp[i] = 1e18;h = 1,t = 1;dl[1] = 0;for (int i = 1;i <= n;i++){while (h < t && ju(dl[h],dl[h+1])<dis[i]) h++;int j = dl[h];dp[i] = min(dp[i],dp[j] + c[i] +dis[i]*(sump[i]-sump[j])-(sumdp[i]-sumdp[j]));while (h < t && ju(dl[t-1],dl[t]) > ju(dl[t],i)) t--;dl[++t] = i;}printf("%lld\n",dp[n]);return 0; }轉(zhuǎn)載于:https://www.cnblogs.com/AWCXV/p/7637996.html
總結(jié)
以上是生活随笔為你收集整理的【BZOJ 1096】[ZJOI2007]仓库建设的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: php常用过滤htmlspecialch
- 下一篇: spring的jar包以及相关的API文