http://codeforces.com/contest/319/problem/C
代码太挫,不忍细看...
想着想着,开始胡思乱想...
想着先按顺序砍,例如 1 2 3 4 5
如果发生变化,就会有数字跳到前面来,比如 1 3 2 4 5
这时发现 2 已无卵用,反正不需要b[2]了,干脆放最后好了
然后再变化的时候,若跳到3前面,3就是白跳了
问题就转化成 从序列中选取若干点,按顺序砍,使到n的花费最小
转移方程就是 dp[i] = min(dp[j] + B[j] * A[i]) (j < i)
然而我并不会斜率优化,顺道学了一下,随便找了一篇
http://www.cnblogs.com/ka200812/archive/2012/08/03/2621345.html
#include
using namespace std;
typedef long long ll;
#define rep(i,n) for(int i = 0;i < n;i ++)
#define repe(i,n) for(int i=1;i<=n;i++)
#define mst(A,k) memset(A,k,sizeof(A))
const ll INF = 1e18;
int n;
ll A[100005],B[100005],dp[100005];
int que[200005],st,ed;
int main() {
scanf("%d",&n);
for(int i=1; i<=n; i++) scanf("%d",A+i);
for(int i=1; i<=n; i++) scanf("%d",B+i);
dp[1] = 0;
que[st = 0] = 1;
ed = 1;
for(int i=2; i<=n; i++)
{
while(ed - st > 1 && (dp[que[st+1]] - dp[que[st]]) * 1.0 <= A[i] * 1.0 * (B[que[st]] - B[que[st+1]]))
{
st++;
}
dp[i] = dp[que[st]] + B[que[st]] * M * A[i];
while(ed - st > 1 && (dp[que[ed-1]] - dp[que[ed-2]]) * 1.0 * (B[i] - B[que[ed-1]]) < (dp[i] - dp[que[ed-1]]) * 1.0 * (B[que[ed-1]] - B[que[ed-2]]))
{
ed--;
}
que[ed++] = i;
}
cout<
<<"\n";
return 0;
}