AtCoder Beginner Contest 207 E - Mod i(dp)

题意:

给出一个长度为 n n n的数组 a a a,假如将这个数组分成 k k k个连续段,第 i i i段的和为 b i b_i bi,如果对于每一个 i i i b i % i = = 0 b_i\%i==0 bi%i==0,那么该分割方法称为合法的,求有多少个合法的分割方法。

题解:

假如不考虑复杂度,那么可以想出一个 O ( n 3 ) O(n^3) O(n3)的做法:设 d p [ i ] [ j ] dp[i][j] dp[i][j]表示前 i i i个数分为 j j j段的合法方案数,那么转移方程就是

d p [ i ] [ j ] + = d p [ k ] [ j − 1 ] dp[i][j]+=dp[k][j-1] dp[i][j]+=dp[k][j1] ( s u m ( k + 1 , i ) % j = = 0 ) (sum(k+1,i)\%j==0) (sum(k+1,i)%j==0)

考虑怎么优化,因为 s u m ( k + 1 , i ) % j = = 0 sum(k+1,i)\%j==0 sum(k+1,i)%j==0 所以 ( s u m ( 1 , i ) − s u m ( 1 , k ) ) % j = = 0 (sum(1,i)-sum(1,k))\%j==0 (sum(1,i)sum(1,k))%j==0

又可以推出 s u m ( 1 , i ) % j = = s u m ( 1 , k ) % j sum(1,i)\%j==sum(1,k)\%j sum(1,i)%j==sum(1,k)%j

那么我们可以先枚举 j j j,再枚举 i i i,每次将 d p [ i ] [ j − 1 ] dp[i][j-1] dp[i][j1]放入数组 n u m [ s u m ( 1 , i ) % j ] num[sum(1,i)\%j] num[sum(1,i)%j], 那么 d p [ i ] [ j ] = n u m [ s u m ( 1 , i ) % j ] dp[i][j]=num[sum(1,i)\%j] dp[i][j]=num[sum(1,i)%j]

为什么呢?上面 n 3 n^3 n3的转移式中,将合适位置的 k k k d p [ k ] [ j − 1 ] dp[k][j-1] dp[k][j1]算入贡献,而 n u m num num数组里面存的正是所有满足要求的 d p [ k ] [ j − 1 ] dp[k][j-1] dp[k][j1]的和。

代码:

#pragma GCC diagnostic error "-std=c++11"
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<stack>
#include<set>
#include<ctime>
#define iss ios::sync_with_stdio(false)
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
typedef pair<int,int> pii;
const int mod=1e9+7;
const int MAXN=3e3+5;
const int inf=0x3f3f3f3f;
int n;
ll a[MAXN];
ll sum[MAXN];
ll dp[MAXN][MAXN];
ll f[MAXN][MAXN];
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        sum[i]=sum[i-1]+a[i];
    }
    f[0][0]=1;
    for(int j=1;j<=n;j++)
    {
        for(int i=1;i<=n;i++)
        {
            dp[i][j]=f[sum[i]%j][j-1];
            f[sum[i]%j][j-1]=(f[sum[i]%j][j-1]+dp[i][j-1])%mod;
        }
    }
    ll ans=0;
    for(int i=1;i<=n;i++)
    {
        ans=(ans+dp[n][i])%mod;
    }
    cout<<ans<<endl;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值