LOJ6077「2017 山东一轮集训 Day7」逆序对 生成函数+组合数+DP

Posted CQzhangyu

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了LOJ6077「2017 山东一轮集训 Day7」逆序对 生成函数+组合数+DP相关的知识,希望对你有一定的参考价值。

【LOJ6077】「2017 山东一轮集训 Day7」逆序对

题目描述

给定 n,k ,请求出长度为 n的逆序对数恰好为 k 的排列的个数。答案对 109+7 取模。

对于一个长度为 n 的排列 p ,其逆序对数即满足 i<j 且 pi>pj 的二元组 (i,j)的数量。

输入格式

一行两个整数 n,k

输出格式

一行,表示答案。

样例输入

7 12

样例输出

531

数据范围与提示

对于 20% 的数据,n,k≤20
对于 40% 的数据,n,k≤100
对于 60% 的数据,n,k≤5000
对于 100% 的数据,$1 \leq n, k \leq 100000, 1 \leq k \leq \binom{n}{2}$

题解:本人第一思路是生成函数,但是想了想模数1e9+7没法搞,后来发现这个思路还真的是对的。(还真的有人拿生成函数A了,太神了)

首先从小到大插入第i个数时,逆序对数可能增加0,1,2,...i-1,所以最终得到的生成函数就是

$f(n)=1\times(1+x)\times(1+x+x^2)\times(1+x+x^2+x^3)...$

$f(n)={\prod\limits_{i=1}^n(1-x^i)\over(1-x)^n}$

下面那个东西很好求,${1\over (1-x)^n}=(1+x+x^2+...)^n=\sum C_{i+n-1}^{n-1}x^i$,然后我们考虑上面那个东西有什么意义。

你可以理解为第i项的系数是:有n个数,1,2,3...n,从中选出j个数使得总和为i的方案数$\times(-1)^j$。

这就大大简化了我们的问题,我们令f[j][i]表示选出j个数总和为i的方案数,显然j是$\sqrt{i}$级别的。

但是我们选出来的j个数并不能重复,所以这个问题还是比较难处理的,我们可以再转化一下,求长度为j,每个数在[1,n]之间,总和为i的上升序列的方案数。

如何构造出所有的上升序列呢?我们考虑将这个序列逆向差分$(b_i=a_i-a_{i+1})$,于是这个序列的总和就变成了$\sum\limits_{k=1}^jb_k\times k$。我们只需要满足$b_k>0$即可。

这时就容易DP了,f[i][j]可以由这几种状态转移而来:

如果i>=j,我们可以将bj++,那么f[j][i]+=f[j][i-j];我们还可以在bj后面增加一个1,那么f[j][i]+=f[j-1][i-j]。
如果i>n,此时可能出现a序列的最后一项>n的情况,即b序列的总和>n的情况,那么f[j][i]-=f[j-1][i-n-1]即可。

最后统计一下答案即可,时间复杂度$O(k\sqrt{k})$

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
const int N=100010;
const int M=450;
const int P=1000000007;
typedef long long ll;
ll jc[N<<1],jcc[N<<1],ine[N<<1];
ll ans;
int f[M][N];
int n,k;
inline ll c(int a,int b)
{
	if(a<b)	return 0;
	return jc[a]*jcc[b]%P*jcc[a-b]%P;
}
inline void upd(int &x,int y)
{
	x+=y;
	if(x>=P)	x-=P;
}
int main()
{
	scanf("%d%d",&n,&k);
	int i,j;
	ine[0]=ine[1]=jc[0]=jc[1]=jcc[0]=jcc[1]=1;
	for(i=2;i<=n+k;i++)	jc[i]=jc[i-1]*i%P,ine[i]=P-(P/i)*ine[P%i]%P,jcc[i]=jcc[i-1]*ine[i]%P;
	f[0][0]=1;
	for(i=1;i<M;i++)
	{
		for(j=i;j<=k;j++)
		{
			if(j>=i)	upd(f[i][j],f[i][j-i]),upd(f[i][j],f[i-1][j-i]);
			if(j>n)	upd(f[i][j],P-f[i-1][j-n-1]);
		}
	}
	for(i=0;i<=k;i++)
	{
		ll tmp=0;
		for(j=0;j<M;j++)	tmp+=((j&1)?-1:1)*f[j][i];
		tmp=(tmp%P+P)%P;
		ans=(ans+tmp*c(k-i+n-1,n-1))%P;
	}
	printf("%lld",ans);
	return 0;
}

 

以上是关于LOJ6077「2017 山东一轮集训 Day7」逆序对 生成函数+组合数+DP的主要内容,如果未能解决你的问题,请参考以下文章

loj6102 「2017 山东二轮集训 Day1」第三题

Loj #6142. 「2017 山东三轮集训 Day6」A

「2017 山东一轮集训 Day6」子序列

「2018山东一轮集训」 Tree

[LOJ#2328]「清华集训 2017」避难所

「6月雅礼集训 2017 Day7」电报