bzoj2431:[HAOI2009]逆序对数列[DP]

最新推荐文章于 2022-10-21 14:46:57 发布

AC率最低的傻X

最新推荐文章于 2022-10-21 14:46:57 发布

阅读量1.5k

点赞数

分类专栏：动态规划

本文链接：https://blog.csdn.net/baidu_20126217/article/details/39505947

版权

动态规划专栏收录该内容

9 篇文章 0 订阅

订阅专栏

首先为了使状态更方便转移，我们显然要从大到小或从小到大往当前序列插入进行转移。

我们尝试在一个已经有j个逆序对和前i-1个数字组成的数列中插入，我们会发现将数字i插入可以很轻松的算出新的逆序对贡献。

作为这个数列中最大的数字，显然呆在它后面的都是逆序对咯。

于是得到动归方程:

f[i][j]=Sigma(f[i-1][j-i+k]) k∈[1,i]

首先发现一个性质，用的总是i-1状态，那么我们滚动一下就轻松剩下了一大堆空间。（直接去掉了一维）

然后重点来了！原来这个DP是O(n^3)的，但是这个递推式居然是用连续区间和，那么加个前缀和维护不就秒杀了。

------------------------------------------------------------------------------------------------------------------

Code:

#include<iostream>
#include<cstdio>
#include<cstring>
#define For(i,l,r) for(int (i)=(l);(i)<=(r);(i)++)
using namespace std;
const int N=1010,MOD=10000;
int n,k,f[N],sum[N];
int main(){
    scanf("%d%d",&n,&k);
    For(i,0,k) sum[i]=1;
    For(i,2,n){
        f[0]=1;
        For(j,1,k){
            if(j-i>=0){
                f[j]=(sum[j]-sum[j-i])%MOD;
                while(f[j]<0) f[j]+=MOD; 
            }else f[j]=sum[j]%MOD;
        }
        For(j,1,k){
            sum[j]=sum[j-1]+f[j];
            sum[j]%=MOD;
        }
    }printf("%d\n",f[k]);
    return 0;
}

要注意一下算术溢出，我也是对拍才发现的QAQ