【AGC030F】Permutation and Minimum(DP)

题目链接

题解

首先可以想到分组后,去掉两边都填了数的组。
然后就会剩下\((-1,-1)\)和\((-1,x)\)或\((x,-1)\)这两种情况
因为是最小值序列的情况数,我们可以考虑从大到小填数,一个组填完了最小值就是当前填的数
设\(f[i][j][k]\)表示已经填到数\(i\)剩余\(j+k\)个填了一个数的组,其中有\(j\)个组是第一种情况填出来的,\(k\)个是第二种情况来的
分情况转移
如果这个数属于第一种情况,那么可以转移到:
\(f[i][j+1][k]:\)填到一个\((-1,-1)\)中
\(f[i][j-1][k]:\)和之前填了一个数的\((-1,-1)\)完成一组
\(f[i][j][k-1]:\)和\((x,-1)\)完成一组
如果属于第二种情况,那么可以转移到:
\(f[i][j][k+1]:\)不是最小值
\(f[i][j-1][k]:\)是最小值
最后因为\((-1,-1)\)是无序的,所以答案要乘上一个阶乘

Code

int cnt1 = 0, cnt2 = 0;
    for (int i = 1; i < 2 * n; i += 2) {
        if (a[i] == a[i + 1])
            cnt2++;
        else {
            if (a[i] > a[i + 1]) swap(a[i], a[i + 1]);
            if (a[i] == -1) cnt1++, flg[a[i + 1]] = 1;
            else vis[a[i]] = vis[a[i + 1]] = 1;
        }
    }
    int m = 0;
    f[0][0][0] = 1;
    for (int i = 2 * n; i >= 1; i--)
        if (!vis[i]) {
            m++;
            if (!flg[i]) {
                for (int j = 0; j <= cnt1 + cnt2; j++)
                    for (int k = 0; k <= cnt1; k++) {
                        pls(f[m][j + 1][k], f[m - 1][j][k]);
                        if (j) pls(f[m][j - 1][k], f[m - 1][j][k]);
                        if (k) pls(f[m][j][k - 1], 1ll * k * f[m - 1][j][k] % Mod);
                    }
            }
            else {
                for (int j = 0; j <= cnt1 + cnt2; j++)
                    for (int k = 0; k <= cnt1; k++) {
                        pls(f[m][j][k + 1], f[m - 1][j][k]);
                        if (j) pls(f[m][j - 1][k], f[m - 1][j][k]);
                    }
            }
        }
    int ans = f[m][0][0];
    for (int i = 1; i <= cnt2; i++)
        ans = 1ll * ans * i % Mod;

原文地址:https://www.cnblogs.com/zzy2005/p/12005563.html

时间: 2024-11-26 01:42:44

【AGC030F】Permutation and Minimum(DP)的相关文章

【BZOJ】1068: [SCOI2007]压缩(dp)

http://www.lydsy.com/JudgeOnline/problem.php?id=1068 发现如果只设一维的话无法转移 那么我们开第二维,发现对于前i个来说,如果确定了M在哪里,第i个是用R还是不用就能确定了(如果用R那么在中间一定变成了缓冲串) 那么可以转移了 设d[i,j]表示前i个串,最近的一个M在i的前边一个格子,的最短长度,有 d[1,1]=1 d[i,i]=min{d[i-1,j]}+2 //即用一次M又补上i,所以+2 d[i,j]=d[pos,j]+1,其中pos

【bzoj1925】地精部落[SDOI2010](dp)

题目传送门:1925: [Sdoi2010]地精部落 这道题,,,首先可以一眼看出他是要我们求由1~n的排列组成,并且抖来抖去的序列的方案数.然后再看一眼数据范围,,,似乎是O(n^2)的dp?然后各种撕烤,,,然而还是不会... 对于这道题,我第一眼的想法是用f[i][j]表示长度为i,最后一个数是j的抖动序列的方案数,然而这是个1~n的排列,似乎没法解决数字重复的问题..QAQ 于是看了一波题解,,,原来这个dp是这样子搞的:用f[i][j]表示i个数的排列.第一个数<=j且开头下降的抖动序

【51nod1519】拆方块[Codeforces](dp)

题目传送门:1519 拆方块 首先,我们可以发现,如果第i堆方块被消除,只有三种情况: 1.第i-1堆方块全部被消除: 2.第i+1堆方块全部被消除:(因为两侧的方块能够保护这一堆方块在两侧不暴露) 3.第i堆方块过了h[i]次操作后,从上到下被消除. 于是我们设l[i]为第i堆方块从左边开始消除的最小操作次数,设r[i]从右边开始消除的最小操作次数. 然后从左向右dp出l[i],从右向左dp出r[i],然后就能算答案了. 代码: #include<cstdio> #include<cs

【vijos】1764 Dual Matrices(dp)

https://vijos.org/p/1764 自从心态好了很多后,做题的确很轻松. 这种题直接考虑我当前拿了一个,剩余空间最大能拿多少即可. 显然我们枚举每一个点拿出一个矩形(这个点作为右下角),然后剩余空间我们只需要考虑i+1~n和j+1~m这一块空间即可(至于为什么多想想QAQ) 所以我们维护i+1~n和j+1~m能取到的最大矩形即可. 显然二维前缀和,然后再维护一个二维 mx[i][j]=max{get(i, j), mx[i+1][j], mx[i][j+1]},然后每一次找一个点i

【uva-11584】Partitioning by Palindromes(dp)

粗略的复杂度是L^3,长度最大是1000,,没敢做,之后发现其实这个复杂度的系数也不大,可以过,而且很快. dp[j] = dp[i - 1] + 1 (if(str[i] ~ str[j]为回文) 14327451 11584 Partitioning by Palindromes Accepted C++ 0.052 2014-10-09 09:33:17 #include<set> #include<map> #include<stack> #include<

【vijos】1757 逆序对(dp)

https://vijos.org/p/1757 有时候自己sb真的是不好说... 我竟然想了半天都没想到这个转移. 我是有多傻.... 我们设f[i][j]表示1~i的排列且逆序对恰好是j的方案数. 显然我们只需要将i放进i-1排列中就行了. 而且发现i始终大于i-1 那么就好做了,我们只要将所有i放到i-1序列的位置的方案全都加起来即可,即: f[i][j]=sum{f[i-1][k], max{0, j-i+1}<=k<=j} 用前缀和搞搞就行了. #include <cstdio

【Luogu】P2331最大子矩阵(DP)

题目链接 这题的状态转移方程真是粗鄙. f[i][j][k]表示前i行用了j个矩阵状态为k的时候的最大值. k=0:两列都不选. k=1:取左弃右. k=2:选右弃左. k=3:左右都选,但分属于两个独立矩阵. k=4:左右都选,且同属于一个矩阵. 参考题解:孤寂的时代 代码 #include<cstdio> #include<cctype> #include<cstring> inline long long max(long long x,long long y){

【UVa】Partitioning by Palindromes(dp)

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=27&page=show_problem&problem=2549 设w[i,j]为i-j能分割成的最少回文串 f[i]为前i个字符能够分成的最少回文串 w[i,j]=1 当w[i+1,j-1]==1 && s[i]==s[j] 或 i==j-1 && s[i]==s[j] w[i,j]=

【Luogu】P2051中国象棋(DP)

题目链接 去看STDCALL的题解吧 #include<cstdio> #include<cctype> #define mod 9999973 inline long long read(){ long long num=0,f=1; char ch=getchar(); while(!isdigit(ch)){ if(ch=='-') f=-1; ch=getchar(); } while(isdigit(ch)){ num=num*10+ch-'0'; ch=getchar(