题目描述

有n个木块排成一行,从左到右依次编号为1~n。你有k种颜色的油漆,其中第i种颜色的油漆足够涂ci个木块。所有油漆刚好足够涂满所有木块,即c1+c2+…+ck=n。相邻两个木块涂相同色显得很难看,所以你希望统计任意两个相邻木块颜色不同的着色方案。

分析

很多网上的题解都讲的不怎么清楚,我将个人思考到一发WAAC的过程清楚地讲一遍。
首先考虑如何定义状态,一眼就发现了\(c_i<=5\),说明什么,我们可以直接将各个颜色还剩下多少的个数都压到状态中,也就是一个五维的状态。而且题目很明显的发现如果定义各个颜色还有多少个,那么需要用15维,数组直接爆炸。
但是题目要求不能相邻,因为我们已经采取了用个数作为前面的5个状态,那么参考chess的思路,排列组合来排除是否是相邻的,这个待会再说。说的明白一点,就是在定义\(1\)维状态,表示上一个格子我们用的是个数是\(lst\)颜色。
总结一下:\(f[i1][i2][i3][i4][i5][lst]\),\(i_i\)表示的是\(i\)颜色还有\(i_i\)个,前一个涂得颜色是\(lst\)个颜色时方案总数。

我们把这个放棋子的过程,按照从左到右的顺序。


接下来我们思考决策。我们举3个例子,分别是我们放的颜色是\(1\)和\(3\)和\(5\)。
如果我们当前放的颜色是个数为\(1\)的,那么分成两种情况:

  • 如果lst=2,那么说明原来是一个是从\(2\)个颜色中取出了一个颜色,变成了个数为\(1\)的颜色,从而得出我们现在放的颜色就是现在有\(1\)个的颜色,那么不能重复,就说明剩下来的\(i1-1\)个一个的颜色都可以放到这个颜色中。
  • 如果不是2的话,就和上面的情况刚好相反,也就是之前放的颜色不会和我们现在放的颜色不会冲突,乘法原理可以得到方程。

总结一下转移方程
\[if (lst=2) f[i1][i2][i3][i4][i5][lst] += (i1 - 1) * f[i1-1][i2][i3][i4][i5][1]\]
否则
\[f[i1][i2][i3][i4][i5][lst] += i1*f[i1-1][i2][i3][i4][i5][1]\]


如果我们放的颜色个数为\(3\),这种情况和\(2\)以及\(4\)是相同的,以3为例:

  • 如果lst=4,那么参照前面的思路,现在这个颜色就是个数为\(3\)个的颜色,那么如果要将\(3\)放进去,那么会是剩下来个数为\(2\)的颜色个数会\(+1\),让个数为\(3\)的颜色个数\(-1\),
  • 反之也就是之前放的颜色不会和我们现在放的颜色不会冲突,乘法原理可以得到方程。

总结一下转移方程
\[if (lst=4) f[i1][i2][i3][i4][i5][lst]+=(i3-1)*f[i1][i2+1][i3-1][i4][i5][3]\]
否则
\[f[i1][i2][i3][i4][i5][lst]+=i3*f[i1][i2+1][i3-1][i4][i5][3]\]


如果我们放的个数是\(5\):
因为不可能有颜色可以放了之后还有\(5\)个,那么就直接乘法原理就可以的到房成了。
总结一下转移方程
\[f[i1][i2][i3][i4][i5][lst] += i5*f[i1-1][i2][i3][i4][i5][5]\]


状态转移方程完成了,那么接下来就是状态的边界条件。很清楚的是,一开始的时候无论取哪一个颜色方案数都只有\(1\),那么也就是\(f[0][0][0][0][0][i]=1\),其他都赋值为0。
因为考虑到这道题目线性推常数较大,就采用了记忆化搜索来实现DP。

ac代码

#include <bits/stdc++.h>
#define ll long long
#define ms(a, b) memset(a, b, sizeof(a))
#define inf 0x3f3f3f3f
#define mod 1000000007
using namespace std;
template <typename T>
inline void read(T &x) {
    x = 0; T fl = 1;
    char ch = 0;
    while (ch < '0' || ch > '9') {
        if (ch == '-') fl = -1;
        ch = getchar();
    }
    while (ch >= '0' && ch <= '9') {
        x = (x << 1) + (x << 3) + (ch ^ 48);
        ch = getchar();
    }
    x *= fl;
}
ll f[16][16][16][16][16][6];
int d[6];
int n;
ll dfs(int i1, int i2, int i3, int i4, int i5, int lst) {
    if (f[i1][i2][i3][i4][i5][lst]) return f[i1][i2][i3][i4][i5][lst];
    ll res = 0;
    if (i1 > 0) {
        if (lst == 2) res = (res + (i1 - 1) * dfs(i1 - 1, i2, i3, i4, i5, 1)) % mod;
        else res = (res + i1 * dfs(i1 - 1, i2, i3, i4, i5, 1)) % mod;
    }
    if (i2 > 0) {
        if (lst == 3) res = (res + (i2 - 1) * dfs(i1 + 1, i2 - 1, i3, i4, i5, 2)) % mod;
        else res = (res + i2 * dfs(i1 + 1, i2 - 1, i3, i4, i5, 2)) % mod;
    }
    if (i3 > 0) {
        if (lst == 4) res = (res + (i3 - 1) * dfs(i1, i2 + 1, i3 - 1, i4, i5, 3)) % mod;
        else res = (res + i3 * dfs(i1, i2 + 1, i3 - 1, i4, i5, 3)) % mod;
    }
    if (i4 > 0) {
        if (lst == 5) res = (res + (i4 - 1) * dfs(i1, i2, i3 + 1, i4 - 1, i5, 4)) % mod;
        else res = (res + i4 * dfs(i1, i2, i3 + 1, i4 - 1, i5, 4)) % mod;
    }
    if (i5 > 0) {
        res = (res + i5 * dfs(i1, i2, i3, i4 + 1, i5 - 1, 5)) % mod;
    }
    return f[i1][i2][i3][i4][i5][lst] = res % mod;
}
int main() {
    int k;
    read(k);
    for (int i = 1; i <= k; i ++) {
        int x;
        read(x);
        d[x] ++;
    }
    for (int i = 1; i <= 5; i ++) f[0][0][0][0][0][i] = 1;
    printf("%lld\n", dfs(d[1], d[2], d[3], d[4], d[5], 0) % mod);
    return 0;
}

[luogu2476][bzoj1079][SCOI2008]着色方案【动态规划】的更多相关文章

  1. BZOJ1079 [SCOI2008]着色方案 动态规划

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ1079 题目概括 有n个木块排成一行,从左到右依次编号为1~n.你有k种颜色的油漆,其中第i种颜色的 ...

  2. BZOJ1079:[SCOI2008]着色方案(DP)

    Description 有n个木块排成一行,从左到右依次编号为1~n.你有k种颜色的油漆,其中第i种颜色的油漆足够涂ci个木块. 所有油漆刚好足够涂满所有木块,即c1+c2+...+ck=n.相邻两个 ...

  3. BZOJ1079 [SCOI2008]着色方案 【dp记忆化搜索】

    题目 有n个木块排成一行,从左到右依次编号为1~n.你有k种颜色的油漆,其中第i种颜色的油漆足够涂ci个木块. 所有油漆刚好足够涂满所有木块,即c1+c2+-+ck=n.相邻两个木块涂相同色显得很难看 ...

  4. BZOJ1079: [SCOI2008]着色方案 (记忆化搜索)

    题意:有n个木块排成一行,从左到右依次编号为1~n.你有k种颜色的油漆,其中第i种颜色的油漆足够涂ci个木块. 所有油漆刚好足够涂满所有木块,即c1+c2+...+ck=n.相邻两个木块涂相同色显得很 ...

  5. 2018.10.20 bzoj1079: [SCOI2008]着色方案(多维dp)

    传送门 dp妙题. f[a][b][c][d][e][last]f[a][b][c][d][e][last]f[a][b][c][d][e][last]表示还剩下aaa个可以用一次的,还剩下bbb个可 ...

  6. BZOJ1079 [SCOI2008]着色方案[组合计数DP]

    $有a_{1}个1,a_{2}个2,...,a_{n}个n(n<=15,a_{n}<=5),求排成一列相邻位不相同的方案数.$ 关于这题的教训记录: 学会对于复杂的影响分开计,善于发现整体 ...

  7. bzoj1079: [SCOI2008]着色方案

    dp.以上次染色时用的颜色的数量和每种数量所含有的颜色作状态. #include<cstdio> #include<algorithm> #include<cstring ...

  8. 【记忆化搜索】bzoj1079 [SCOI2008]着色方案

    #include<cstring> #include<cstdio> using namespace std; #define MOD 1000000007 typedef l ...

  9. bzoj1079: [SCOI2008]着色方案

    ci<=5直接想到的就是5维dp了...dp方程YY起来很好玩...写成记忆化搜索比较容易 #include<cstdio> #include<cstring> #inc ...

随机推荐

  1. 2019省赛训练组队赛3.26周二---FJUT 2016

    A.Minimum’s Revenge There is a graph of n vertices which are indexed from 1 to n. For any pair of di ...

  2. PHP导出CSV文件出现乱码的解决方法

    在做项目时碰到使用外语的情况下,我们就会使用UTF-8编码.但是,在用PHP导出CSV文件时,如果写入的数据是使用UTF-8编码的日语.韩语之类的外文,就会出现乱码. 要解决PHP生成CSV文件的乱码 ...

  3. 6 Prefer and Would rather

    1 prefer 使用 "prefer" 用来表明通常喜欢某件事甚于另一件事.说话者喜欢打高尔夫球更甚于喜欢打网球."prefer" 的后面可以接名词(&quo ...

  4. Windows 下面 redis 发布为服务的官方方法

    除了 NSSM 之外 另外一种方式 感觉还是很好用的 redis-server --service-install redis.windows.conf --loglevel verbose 感觉也可 ...

  5. python3 写的一个压测脚本(有待开发)

    import requests import queue import threading import time status_code_list = [] exec_time = 0 class ...

  6. Guava Cache源码详解

    目录 一.引子 二.使用方法 2.1 CacheBuilder有3种失效重载模式 2.2 测试验证 三.源码剖析 3.1 简介 3.2 源码剖析 四.总结 优点: 缺点: 正文 回到顶部 一.引子 缓 ...

  7. 老男孩python学习自修第十九天【异常处理】

    1.常见的错误 TypeError 类型错误 NameError 没有该变量 ValueError 不期望的值 AttributeError 没有该属性 UnboundLocalError 没有该局部 ...

  8. tomcat 与 nginx,apache的区别

    tomcat 与 nginx,apache的有什么区别 回答一: 题主说的Apache,指的应该是Apache软件基金会下的一个项目——Apache HTTP Server Project:Nginx ...

  9. delphi怎样在关闭程序时弹出窗口?

    我想在关闭delphi编译的程序时,弹出“您是否确实要退出的窗口”点击否不退出,是退出 在主窗体的CloseQuery事件里,使用messagebox进行提示,根据选择的按钮对Canclose进行设置 ...

  10. java 中的Collection

    /* *一. Collection?-------->容器! * * 1.来源于java.util包 非常实用的数据结构! * *二. 方法? * * void clear()删除集合中所有元素 ...