zoukankan      html  css  js  c++  java
  • POJ 1821 Fence (算竞进阶习题)

    单调队列优化dp

    我们把状态定位F[i][j]表示前i个工人涂了前j块木板的最大报酬(中间可以有不涂的木板)。

    第i个工人不涂的话有两种情况: 那么F[i - 1][j], F[i][j - 1]就成为了转移状态的候选。

    那如果第i个工人要涂的话,我们可以假设我们是从k+1涂到j的,根据题意可以求出k的取值范围,然后状态转移的条件限制了j的取值范围。

    我们考虑每j从小到大增加的过程,j对应的k的取值是一个上界不变下届变大的区间,是一个滑动窗口,那我们可以用单调队列来维护决策k的最优候选。

    #include <bits/stdc++.h>
    #define INF 0x3f3f3f3f
    #define full(a, b) memset(a, b, sizeof a)
    using namespace std;
    typedef long long ll;
    inline int lowbit(int x){ return x & (-x); }
    inline int read(){
        int X = 0, w = 0; char ch = 0;
        while(!isdigit(ch)) { w |= ch == '-'; ch = getchar(); }
        while(isdigit(ch)) X = (X << 3) + (X << 1) + (ch ^ 48), ch = getchar();
        return w ? -X : X;
    }
    inline int gcd(int a, int b){ return a % b ? gcd(b, a % b) : b; }
    inline int lcm(int a, int b){ return a / gcd(a, b) * b; }
    template<typename T>
    inline T max(T x, T y, T z){ return max(max(x, y), z); }
    template<typename T>
    inline T min(T x, T y, T z){ return min(min(x, y), z); }
    template<typename A, typename B, typename C>
    inline A fpow(A x, B p, C lyd){
        A ans = 1;
        for(; p; p >>= 1, x = 1LL * x * x % lyd)if(p & 1)ans = 1LL * x * ans % lyd;
        return ans;
    }
    const int N = 20000;
    int n, m, dp[200][N], q[N];
    struct Node {
        int l, s, p;
        bool operator < (const Node &rhs) const {
            return s < rhs.s;
        }
    }w[N];
    
    int calc(int i, int k){
        return dp[i - 1][k] - w[i].p * k;
    }
    
    int main(){
    
        while(~scanf("%d%d", &n, &m)){
            full(q, 0), full(dp, 0);
            for(int i = 1; i <= m; i ++)
                w[i].l = read(), w[i].p = read(), w[i].s = read();
            sort(w + 1, w + m + 1);
            dp[0][0] = 0;
            for(int i = 1; i <= m; i ++){
                int l = 1, r = 0;
                for(int k = max(0, w[i].s - w[i].l); k <= w[i].s - 1; k ++){
                    while(l <= r && calc(i, q[r]) <= calc(i, k)) r --;
                    q[++r] = k;
                }
                for(int j = 1; j <= n; j ++){
                    dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
                    if(j >= w[i].s){
                        while(l <= r && j - q[l] > w[i].l) l ++;
                        if(l <= r) dp[i][j] = max(dp[i][j], w[i].p * j + calc(i, q[l]));
                    }
                }
            }
            printf("%d
    ", dp[m][n]);
        }
        return 0;
    }
    
  • 相关阅读:
    创造别人想要的东西(2)----创业15问
    创造别人想要的东西(1)----创业的本质
    有点迷茫
    普林斯顿大学算法公开课(1)----介绍
    工业控制系统信息安全与信息安全、功能安全的区别【转自工控工程网】
    十条概念区分信息安全与隐私保护【转自安全牛】
    NBT(NetBIOS Over TCP)名称解析概述
    SMB协议概述
    PCI DSS合规建设ASV扫描介绍
    信息安全等级测评师(初级 技术)
  • 原文地址:https://www.cnblogs.com/onionQAQ/p/10758228.html
Copyright © 2011-2022 走看看