“蔚来杯“2022牛客暑期多校训练营7 J： Melborp Elcissalc

sliver-awa

已于 2022-08-14 13:37:50 修改

阅读量238

点赞数 1

文章标签：算法 c++ 动态规划

于 2022-08-09 16:46:30 首次发布

本文链接：https://blog.csdn.net/Penny_rd/article/details/126248805

版权

博客详细解析了“蔚来杯”2022年牛客暑期多校训练营第七场的J题——Melborp Elcissalc。由于纯暴力会导致超时，作者采用了前缀和与动态规划的方法来解决。通过计算sum数组并求模k，确定合法区间，然后用四维动态规划求解方案数，最终实现勉强通过的时间复杂度。

摘要由CSDN通过智能技术生成

"蔚来杯"2022牛客暑期多校训练营7 J： Melborp Elcissalc

原题：https://ac.nowcoder.com/acm/contest/33192/J

题目大意

在这里插入图片描述

题解

首先纯纯暴力肯定不行，必会超时
先求前缀和，以sum表示，对于每个sum取模k，不影响结果。即
$sum_i=(sum_i-_1+a_i)\%k$
若区间 $[i, j]$ 中 $sum_i=sum_j$ ，则该区间的和为 $k$ 的倍数，不难证明
接下来，设 $num_i$ 表示 $sum_x=i$ 的数量，则方案数为在这里插入图片描述

然后就可以用dp求方案数了
设 $dp_i,_j,_k$ ，表示已用了 $[0, i]$ 区间内的数，填了 $i$ 个位置，合法数为 $k$ 的方案数，则转移式为
在这里插入图片描述

其中 $s=num_i$
所以最终结果为 $ans=dp_k-_1,_n,_t$
进行四维暴力枚举，时间复杂度 $\large O(64^5)$ 大约等于 $10^9$ ，勉强能过

参考代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=110;
const int M=(1<<12)+10;
const int mod=998244353;
ll c[N][N];
ll dp[N][N][M];
ll n,k,t;
void init(){
    for(int i=0;i<=100;i++) {
        for (int j=0;j<=i;j++) {
            if(!j)c[i][j]=1;
            else c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
        }
    }
    for(int j=0;j<=n;j++)
        dp[0][j][c[j+1][2]]=1;
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    cin>>n>>k>>t;
    init();
    for (int i=1;i<k;i++)
    for (int j=0;j<=n;j++)
    for (int u=0;u<=t;u++)
    for (int l=0;l+j<=n;l++)
    (dp[i][j+l][u+c[l][2]]+=dp[i-1][j][u]*c[j+l][l])%=mod;
    cout<<dp[k-1][n][t];
    return 0;
}