Codeforces Round #189 (Div. 1) C - Kalila and Dimna in the Logging Industry 斜率优化dp

C - Kalila and Dimna in the Logging Industry

很容易能得到状态转移方程 dp[ i ] = min( dp[ j ] + b[ j ] * a[ i ] )， 然后斜率优化一下。

一直以为炸精度了， 忽然发现手贱把while 写成了if 。。。。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define PLL pair<LL, LL>
#define PLI pair<LL, int>
#define PII pair<int, int>
#define SZ(x) ((int)x.size())
#define ull unsigned long long
using namespace std;

const int N = 4e5 + ;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9 + ;
const double eps = 1e-;

int n, que[N], head, rear;
LL a[N], b[N], dp[N];

double calc(int k, int j) {
    return 1.0 * (dp[j] - dp[k]) * (b[k] - b[j]);
}

int main() {
    scanf("%d", &n);
    for(int i = ; i <= n; i++) scanf("%lld", &a[i]);
    for(int i = ; i <= n; i++) scanf("%lld", &b[i]);
    dp[] = ;
    head = , rear = ;
    que[++rear] = ;
    for(int i = ; i <= n; i++) {
        while(rear - head +  >=  && calc(que[head], que[head + ]) < a[i]) head++;
        dp[i] = dp[que[head]] + b[que[head]] * a[i];
        while(rear - head +  >=   && calc(que[rear-], que[rear]) >= calc(que[rear], i)) rear--;
        que[++rear] = i;
    }
    printf("%lld\n", dp[n]);
    return ;
}

/*
*/