ZOJ 2630 Plane Partition(轮廓线 状态压缩DP)

最新推荐文章于 2024-09-08 10:19:43 发布
最新推荐文章于 2024-09-08 10:19:43 发布
阅读量1.7k 收藏 1
点赞数 2
分类专栏： ACM_ZOJ ACM_DP
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/jayye1994/article/details/12357429
版权

题意：

给你一个a*b的矩阵，往里面填充非负整数且满足

x(i, j) <= x(i-1, j)

x(i, j) <= x(i, j-1)

问最大的数字不超过c的填充方式有多少种。 0 < a, b, c <= 6


解题思路：

构造一个像插头DP那样的轮廓线，逐格递推DP。如下图



红色的就是轮廓线，蓝色斜线的格子就是当前决策的格子，轮廓线以上的填充数字都已经确定了，而对轮廓线一下的数字的选择填充有影响的只有轮廓线以上的黄色斜线的格子，所以有这样的DP，

dp[ x ][ y ][ zt ] 表示从格子x，y开始填充数字且轮廓线以上所有数字的压缩状态为zt的填充方式数。

我是直接把x, y, zt (轮廓线上每个格子的的数字)压缩成一个7进制的数，复杂度为 O (7^8) ，需要打表输出。


#include <stdio.h>
#include <string.h>
typedef long long ll;

const int maxn = 6000000;

ll dp[maxn];
int a, b, c, q[11];
ll ww[11][11][11];

int get(int x, int y, int q[]) {
    int ret = x;
    ret = ret*7 + y;
    for(int i = 0;i < b; i++)
        ret = ret*7 + q[i];
    return ret;
}

ll dfs(int zt) {
    int prezt = zt;
    for(int i = b-1;i >= 0; i--) {
        q[i] = zt % 7;
        zt /= 7;
    }
    int x, y;
    y = zt % 7;
    x = zt / 7;
    if(y == b) {
        x++;
        y = 0;
        return dfs( get(x, y, q) );
    }
    if(x == a)
        return 1;

    if(dp[prezt] != -1) return dp[prezt];
    int tmp = q[y];
    if(y > 0 && tmp > q[y-1])   tmp = q[y-1];
    ll ret = 0;
    int ss[11];
    for(int i = 0;i < b; i++)
        ss[i] = q[i];
    for(int i = tmp;i >= 0; i--) {
        ss[y] = i;
        ret += dfs( get(x, y+1, ss) );
    }

    return dp[prezt] = ret;
}

void init() {
ww[1][1][1] = (ll)2;
ww[1][1][2] = (ll)3;
ww[1][1][3] = (ll)4;
ww[1][1][4] = (ll)5;
ww[1][1][5] = (ll)6;
ww[1][1][6] = (ll)7;
ww[1][2][1] = (ll)3;
ww[1][2][2] = (ll)6;
ww[1][2][3] = (ll)10;
ww[1][2][4] = (ll)15;
ww[1][2][5] = (ll)21;
ww[1][2][6] = (ll)28;
ww[1][3][1] = (ll)4;
ww[1][3][2] = (ll)10;
ww[1][3][3] = (ll)20;
ww[1][3][4] = (ll)35;
ww[1][3][5] = (ll)56;
ww[1][3][6] = (ll)84;
ww[1][4][1] = (ll)5;
ww[1][4][2] = (ll)15;
ww[1][4][3] = (ll)35;
ww[1][4][4] = (ll)70;
ww[1][4][5] = (ll)126;
ww[1][4][6] = (ll)210;
ww[1][5][1] = (ll)6;
ww[1][5][2] = (ll)21;
ww[1][5][3] = (ll)56;
ww[1][5][4] = (ll)126;
ww[1][5][5] = (ll)252;
ww[1][5][6] = (ll)462;
ww[1][6][1] = (ll)7;
ww[1][6][2] = (ll)28;
ww[1][6][3] = (ll)84;
ww[1][6][4] = (ll)210;
ww[1][6][5] = (ll)462;
ww[1][6][6] = (ll)924;
ww[2][1][1] = (ll)3;
ww[2][1][2] = (ll)6;
ww[2][1][3] = (ll)10;
ww[2][1][4] = (ll)15;
ww[2][1][5] = (ll)21;
ww[2][1][6] = (ll)28;
ww[2][2][1] = (ll)6;
ww[2][2][2] = (ll)20;
ww[2][2][3] = (ll)50;
ww[2][2][4] = (ll)105;
ww[2][2][5] = (ll)196;
ww[2][2][6] = (ll)336;
ww[2][3][1] = (ll)10;
ww[2][3][2] = (ll)50;
ww[2][3][3] = (ll)175;
ww[2][3][4] = (ll)490;
ww[2][3][5] = (ll)1176;
ww[2][3][6] = (ll)2520;
ww[2][4][1] = (ll)15;
ww[2][4][2] = (ll)105;
ww[2][4][3] = (ll)490;
ww[2][4][4] = (ll)1764;
ww[2][4][5] = (ll)5292;
ww[2][4][6] = (ll)13860;
ww[2][5][1] = (ll)21;
ww[2][5][2] = (ll)196;
ww[2][5][3] = (ll)1176;
ww[2][5][4] = (ll)5292;
ww[2][5][5] = (ll)19404;
ww[2][5][6] = (ll)60984;
ww[2][6][1] = (ll)28;
ww[2][6][2] = (ll)336;
ww[2][6][3] = (ll)2520;
ww[2][6][4] = (ll)13860;
ww[2][6][5] = (ll)60984;
ww[2][6][6] = (ll)226512;
ww[3][1][1] = (ll)4;
ww[3][1][2] = (ll)10;
ww[3][1][3] = (ll)20;
ww[3][1][4] = (ll)35;
ww[3][1][5] = (ll)56;
ww[3][1][6] = (ll)84;
ww[3][2][1] = (ll)10;
ww[3][2][2] = (ll)50;
ww[3][2][3] = (ll)175;
ww[3][2][4] = (ll)490;
ww[3][2][5] = (ll)1176;
ww[3][2][6] = (ll)2520;
ww[3][3][1] = (ll)20;
ww[3][3][2] = (ll)175;
ww[3][3][3] = (ll)980;
ww[3][3][4] = (ll)4116;
ww[3][3][5] = (ll)14112;
ww[3][3][6] = (ll)41580;
ww[3][4][1] = (ll)35;
ww[3][4][2] = (ll)490;
ww[3][4][3] = (ll)4116;
ww[3][4][4] = (ll)24696;
ww[3][4][5] = (ll)116424;
ww[3][4][6] = (ll)457380;
ww[3][5][1] = (ll)56;
ww[3][5][2] = (ll)1176;
ww[3][5][3] = (ll)14112;
ww[3][5][4] = (ll)116424;
ww[3][5][5] = (ll)731808;
ww[3][5][6] = (ll)3737448;
ww[3][6][1] = (ll)84;
ww[3][6][2] = (ll)2520;
ww[3][6][3] = (ll)41580;
ww[3][6][4] = (ll)457380;
ww[3][6][5] = (ll)3737448;
ww[3][6][6] = (ll)24293412;
ww[4][1][1] = (ll)5;
ww[4][1][2] = (ll)15;
ww[4][1][3] = (ll)35;
ww[4][1][4] = (ll)70;
ww[4][1][5] = (ll)126;
ww[4][1][6] = (ll)210;
ww[4][2][1] = (ll)15;
ww[4][2][2] = (ll)105;
ww[4][2][3] = (ll)490;
ww[4][2][4] = (ll)1764;
ww[4][2][5] = (ll)5292;
ww[4][2][6] = (ll)13860;
ww[4][3][1] = (ll)35;
ww[4][3][2] = (ll)490;
ww[4][3][3] = (ll)4116;
ww[4][3][4] = (ll)24696;
ww[4][3][5] = (ll)116424;
ww[4][3][6] = (ll)457380;
ww[4][4][1] = (ll)70;
ww[4][4][2] = (ll)1764;
ww[4][4][3] = (ll)24696;
ww[4][4][4] = (ll)232848;
ww[4][4][5] = (ll)1646568;
ww[4][4][6] = (ll)9343620;
ww[4][5][1] = (ll)126;
ww[4][5][2] = (ll)5292;
ww[4][5][3] = (ll)116424;
ww[4][5][4] = (ll)1646568;
ww[4][5][5] = (ll)16818516;
ww[4][5][6] = (ll)133613766;
ww[4][6][1] = (ll)210;
ww[4][6][2] = (ll)13860;
ww[4][6][3] = (ll)457380;
ww[4][6][4] = (ll)9343620;
ww[4][6][5] = (ll)133613766;
ww[4][6][6] = (ll)1447482465;
ww[5][1][1] = (ll)6;
ww[5][1][2] = (ll)21;
ww[5][1][3] = (ll)56;
ww[5][1][4] = (ll)126;
ww[5][1][5] = (ll)252;
ww[5][1][6] = (ll)462;
ww[5][2][1] = (ll)21;
ww[5][2][2] = (ll)196;
ww[5][2][3] = (ll)1176;
ww[5][2][4] = (ll)5292;
ww[5][2][5] = (ll)19404;
ww[5][2][6] = (ll)60984;
ww[5][3][1] = (ll)56;
ww[5][3][2] = (ll)1176;
ww[5][3][3] = (ll)14112;
ww[5][3][4] = (ll)116424;
ww[5][3][5] = (ll)731808;
ww[5][3][6] = (ll)3737448;
ww[5][4][1] = (ll)126;
ww[5][4][2] = (ll)5292;
ww[5][4][3] = (ll)116424;
ww[5][4][4] = (ll)1646568;
ww[5][4][5] = (ll)16818516;
ww[5][4][6] = (ll)133613766;
ww[5][5][1] = (ll)252;
ww[5][5][2] = (ll)19404;
ww[5][5][3] = (ll)731808;
ww[5][5][4] = (ll)16818516;
ww[5][5][5] = (ll)267227532;
ww[5][5][6] = (ll)3184461423;
ww[5][6][1] = (ll)462;
ww[5][6][2] = (ll)60984;
ww[5][6][3] = (ll)3737448;
ww[5][6][4] = (ll)133613766;
ww[5][6][5] = (ll)3184461423;
ww[5][6][6] = (ll)55197331332;
ww[6][1][1] = (ll)7;
ww[6][1][2] = (ll)28;
ww[6][1][3] = (ll)84;
ww[6][1][4] = (ll)210;
ww[6][1][5] = (ll)462;
ww[6][1][6] = (ll)924;
ww[6][2][1] = (ll)28;
ww[6][2][2] = (ll)336;
ww[6][2][3] = (ll)2520;
ww[6][2][4] = (ll)13860;
ww[6][2][5] = (ll)60984;
ww[6][2][6] = (ll)226512;
ww[6][3][1] = (ll)84;
ww[6][3][2] = (ll)2520;
ww[6][3][3] = (ll)41580;
ww[6][3][4] = (ll)457380;
ww[6][3][5] = (ll)3737448;
ww[6][3][6] = (ll)24293412;
ww[6][4][1] = (ll)210;
ww[6][4][2] = (ll)13860;
ww[6][4][3] = (ll)457380;
ww[6][4][4] = (ll)9343620;
ww[6][4][5] = (ll)133613766;
ww[6][4][6] = (ll)1447482465;
ww[6][5][1] = (ll)462;
ww[6][5][2] = (ll)60984;
ww[6][5][3] = (ll)3737448;
ww[6][5][4] = (ll)133613766;
ww[6][5][5] = (ll)3184461423;
ww[6][5][6] = (ll)55197331332;
ww[6][6][1] = (ll)924;
ww[6][6][2] = (ll)226512;
ww[6][6][3] = (ll)24293412;
ww[6][6][4] = (ll)1447482465;
ww[6][6][5] = (ll)55197331332;
ww[6][6][6] = (ll)1478619421136;
}

int main() {

/*
    for(int i = 1;i <= 6; i++) {
        for(int j = 1;j <= 6; j++) {
            for(int k = 1;k <= 6; k++) {
                a = i; b = j; c = k;
                memset(dp, -1, sizeof(dp));
                for(int l = 0;l < b; l++) {
                    q[l] = c;
                }
                ww[i][j][k] = dfs( get(0, 0, q) );
            }
        }
    }
    for(int i = 1;i <= 6; i++) {
        for(int j = 1;j <= 6; j++) {
            for(int k = 1;k <= 6; k++) {
                printf("ww[%d][%d][%d] = (ll)%I64d;\n", i, j, k, ww[i][j][k]);
            }
        }
    }
*/
    init();
    while(scanf("%d%d%d", &a, &b, &c) != -1) {
        printf("%lld\n", ww[a][b][c]);
    }

    return 0;
}


JayYe
关注
专栏目录
状态压缩DP总结
yhjpku的博客
06-27 1831
参考博客：https://blog.csdn.net/accry/article/details/6607703状态压缩主要指的是用位运算代替枚举压缩DP的时间，如果某一个状态和之前状态的顺序没有关系，那么就可以将之前的选或者不选压缩到一个二进制数中。在选择第i个时枚举之前的所有可能的i-1状态，但是这i-1状态是不记录顺序的，只在i-1到i时考虑顺序，这样的话往往能够节约很多时间。【POJ331...
ZOJ 3963 Heap Partition 贪心
我是谁 --- 卡比兽
03-23 169
ZOJ 3963 题意 给你n个点 叫你用最少的二叉树构造法构造出 父亲节点比儿子节点权值小 然后要坐标大的点做儿子 这道题一开始写的是贪心 普通线段树找比他大的前两个点 然后不断递归下去 实际上是错的 因为链的情况更符合 于是应该写权值线段树找大于他的最小点 这样就可以用set写了 下面是赛中贪心错的代码 #include <iostream> #include <c...
[ACM_动态规划] 轮廓线动态规划——铺放骨牌(状态压缩1)
weixin_34111819的博客
01-18 148
Description Squares and rectangles fascinated the famous Dutch painter Piet Mondriaan. One night, after producing the drawings in his 'toilet series' (where he had to use his toilet paper to ...
【POJ2411】Mondriaan's Dream（轮廓线DP）
小蒟蒻yyb的博客
08-03 291
题面 Vjudge 题解 这题我会大力状压！！！ 时间复杂度大概是O(22nn2)O(22nn2)O(2^{2n}n^2)，设f[i][S]f[i][S]f[i][S]表示当前第iii行向下伸展出去的状态为SSS 那么每次枚举一下当前行的放法，进行转移就好了。 然后就长成了这个样子（不要在意我强行缩减代码长度） 尽管这不是我们本题的重点，然而我还是放份代码 #include&l...
HDU 2167 Pebbles（轮廓线，状态压缩DP）
know_heng的博客
07-24 463
Problem Description You're given an unlimited number of pebbles to distribute across an N x N game board (N drawn from [3, 15]), where each square on the board contains some positive point value be
POJ 2411(轮廓线dp or 状态压缩dp)
weixin_43958964的博客
03-24 481
题目：click 题意：给定一个长，宽为n,m的矩阵，放满1*2的矩形小块问有多少种方法填满。 范围很小，轮廓dp可以做。 k5 k4 k3 k2 k1 状态改写为如上k5k4k3k2k1去表示，1表示这个格子已经放了，0表示没放。（代码中的状态表示也是这个顺序） 什么是轮廓线？ 图中的红色线就是边缘线，由此来进行状态转移。我们对每一个格子进行状态的分析，我们从左上...
轮廓线dp
qq_43983809的博客
01-25 640
首先这道题可以使用常规的dp，枚举第i-1行的状态和第i行的状态进行状态转移。第一次接触轮廓线dp，在网上找一些资料感觉介绍很不清楚，看了这篇。然后是使用轮廓线dp，时间复杂度相比常规状压dp更优。后清晰了一些，主要是根据这道题目。使用图1的状态表示方式源程序。使用图3表示状态的方式源程序。
c语言dp状态转移方程,[总结-动态规划]经典DP状态设定和转移方程
weixin_42511201的博客
05-25 1044
马上区域赛，发现DP太弱，赶紧复习补上。#普通DPCodeForces-546D Soldier and Number Game 筛法+动态规划待补UVALive-8078 Bracket Sequence问以每个字符为左端点的最长括号序列是多长。(包括尖、花、中小括号)状态：设dp[i]为从i开始的括号序列最长长度。转移：以i+1为起点的最长串后边的字符若与左括号匹配，答案是加上这个字符后边的最...
动态规划之状态压缩专题
剑锋OI SharpSwordOI
11-27 1092
态压缩动态规划（简称状压dp）是另一类非常典型的动态规划，通常使用在NP问题的小规模求解中，虽然是指数级别的复杂度，但速度比搜索快，其思想非常值得借鉴。 为了更好的理解状压dp，首先介绍位运算相关的知识。 1.’&’符号，x&y，会将两个十进制数在二进制下进行与运算，然后返回其十进制下的值。例如3(11)&2(10)=2(10)。 2.’|’符号，x|y，会将两个十进制数在二进制下进行或运算，
zoj 1002_zoj1002_
09-29
【压缩包子文件的文件名称列表】仅有一个文件名 "zoj 1002"，这很可能是一个C++源代码文件，可能命名为`zoj1002.cpp`或`zoj1002.cc`，用于解决ZOJ 1002题目。通常，这种程序会包含输入输出处理、问题逻辑分析、算法...
HDU4949 Light (轮廓线dp)
JayYe的专栏
08-15 1806
题意：给你一个01矩阵，有两种操作： 第一种： 把a(i,j)的周围四个都异或一下 第二种： 把a(i, j)的周围四个和a(i,j)都异或一下 求把矩阵变成全0矩阵的最少操作次数 思路：如下图所示的轮廓线dp，逐格递推的，cur为当前决策的格子，红色线就是轮廓线，轮廓线以上的格子的操作状态都已经确定了，而对下面状态有影响的只有黄色格子，每个格子保存的是格子当前的数和它自己操作了多
HDU 1693 Eat the Trees（插头DP，多条回路）
weixin_34152820的博客
10-01 145
  Eat the Trees Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 1663    Accepted Submission(s): 783 Problem Description Most of us know that...
[POJ2411] [UImLocal2000] Mondriaan's Dream [状态压缩/轮廓线][dp]
最后一次修改
10-25 269
[Link\frak{Link}Link] 这道题属于棋盘覆盖问题的变形。 这一类问题在固定行数/列数的情况下是一类十分经典的构造dp； 不论行数/列数，这类问题一般都有明显的轮廓线概念。 同时在数据范围进一步扩大的情况下、可能还可以用矩阵来优化转移。 关于棋盘覆盖问题也有组合学公式，不过不会要求掌握。 这道题在（单纯的）状态压缩和轮廓线的意义下解法是不同的，并且本题的轮廓线解法优于一般...
POJ_P2411 Mondriaan's Dream(轮廓线DP)
BeiYu
01-30 524
Mondriaan’s Dream Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 13983 Accepted: 8086 DescriptionSquares and rectangles fascinated the famous Dutch painter Piet Mondriaan. On
轮廓线DP（插头DP 裸 经典骨牌）
可视化和机器学习
04-27 6849
引言：所谓轮廓线，不是某一行，或者某一列，而是指某一个特定轮廓的状态。 简单的例子： Tiling a Grid With Dominoes http://acm.hdu.edu.cn/showproblem.php?pid=1992 限制了棋盘宽度为4，数据不超过int，只有长度22以内满足答案。 手写状态找规律。 http://www.cnblogs.com/lzsz121
轮廓线dp：GYM103446C
zhangtingxiqwq的博客
11-06 115
考虑轮廓线dp，当我们枚举到蓝色格子的时候，我们记录红色格子的状态。发现枚举过程中还有再维护一个0/1状态d，表示此行有没有向左。当枚举到行末时，我们需要“换行”，把所有3变成1。他们会对所有3和同列的值造成影响。
轮廓线dp、插头dp
最新发布
Se_ren_di_pity的博客
09-08 742
轮廓线：对于一个n*m棋盘上的决策，当我们考虑到(i,j)这个格子时，第i行前j-1个格子的下边以及第j个格子的左边以及第j到第m-1-j个格子的上边所组成的一条线，我们称之为轮廓线轮廓线dp：当我们在棋盘第(i,j)格子上放置棋子时，我们的决策受到且仅受到(i,j)格子左边以及上边的状态的影响，我们这个时候就可以 采取轮廓线dp的方法来解决题目大意：给出n*m的棋盘，其中为1的格子必须铺线，为0的格子禁止铺线，棋盘上可形成多条回路，求所有合法方案数。
【BZOJ】4572: [Scoi2016]围棋-轮廓线DP
远行客
09-15 277
传送门：bzoj4572 题解 算可以匹配的比较麻烦，转化成求不能匹配的，取个补集即可。 c≤6,M≤12c≤6,M≤12c\leq 6,M\leq 12，数据范围明示轮廓线DP，设f[k][sta][i][j](0≤k≤n,0≤sta≤2M,0≤i,j&amp;amp;amp;lt;c)f[k][sta][i][j](0≤k≤n,0≤sta≤2M,0≤i,j&amp;amp;amp;lt;c)f[k][sta][i][j](0\...
div 2 199 D 状态压缩dp 轮廓线很好写
wintowanti
11-05 1671
div 2 199 D 状态压缩dp 轮廓线很好写
C++代码解决ZOJ问题1002, 1091, 1789
"C++代码实现解决ZOJ1002、1091、1789三道编程题" 这段代码是用C++编写的，主要目的是解决ZOJ（在线判题系统ZeroJudge）中的三道题目：1002、1091和1789。代码中使用了深度优先搜索（DFS）算法来寻找最优解。虽然...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值