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;
    }
    
  • 相关阅读:
    计算机程序的构造和解释
    StackOverflow之旅<1>------{去掉烦人的"!=null"判断}
    tomcat 启动显示指定的服务未安装
    Sql Server 表的复制
    Sql Server配置管理器与 Sql Server Management Studio
    “因为数据库正在使用,所以无法获得对数据库的独占访问权。”处理
    Win10磁贴 横向排列 增加多行 多列 磁贴横向展示
    解决使用Hyper-v 设置虚拟机网卡出现-从网络文件夹中隐藏
    Hyper-V网络设置(物理机+虚拟机)
    易语言等软件自动获取管理员权限,在64位Windows7系统非管理员帐户中执行
  • 原文地址:https://www.cnblogs.com/onionQAQ/p/10758228.html
Copyright © 2011-2022 走看看