BZOJ 4402 Claris的剑

Description

Claris想要铸一把剑，这把剑必须符合他的审美，具体来说，我们可以把这把剑的不同地方的宽度看成一个序列，这个序列要满足以下条件：
1.每个元素都是正整数（你的宽度不可能是负数吧）
2.每个元素不能超过M，太宽了如果比Claris身高还高怎么办（你可以认为Claris的身高就是M）
3.相邻两个元素的差的绝对值必须是1（如果是0，则这个地方不是锯齿，杀伤力不够，如果太大，又太丑了）
4.第一个元素的值必须是1（剑尖必须是最窄的地方）
他想知道有多少把长度不超过N（即宽度的序列长度不超过N）的合法的本质不同的剑。
我们认为两把剑本质不同，当且仅当存在至少一个宽度，在两把剑的宽度序列里面出现次数不一样。
比如{1,2,3}和{1,3,2}是本质相同的
{1,2,3}和{1,2,1}则是本质不同的

Input

只有两个整数，表示N,M (数据保证$N,M \leq 2000000$)

Output

只有一个整数，表示答案对$10^9+7$取模的结果

Sample Input

5 3

Sample Output

9

HINT

样例解释



所有本质不同的合法的剑有如下：



{1}

{1,2}

{1,2,3}

{1,2,1}

{1,2,3,2}

{1,2,1,2}

{1,2,3,2,3}

{1,2,3,2,1}

{1,2,1,2,1}

Source

By lyzy

~~~~~~~~~~~~~~~~~~~~~~~~~~

数论+思路~

因为本质不同的情况是至少有一种长度出现次数不同，所以前后顺序是一样的。我们只考虑字典序最小的情况。

对于我们给出的任意m，情况都是 1,2,1,2,3,2,...,m-1,m,m-1,m 或者 1,2,1,2,3,2,...,m-1,m,m-1，后一种可以看成前一种中m=m-1时再添加一项m-1，所以只需要考虑前一种。

前一种可以看成1,2,3,...,m中间插入了(n-m)/2个(i,i+1)组成，而(i,i+1)又有m-1种，所以这里相当于是(m-1)个球放入(n-m)/2个筐中，球可以不放，筐可以为空，所以这部分答案是c((n-m)/2+m-1,m-1)。

取模的过程要强制转化成long long型！

#include<cstdio>
#include<algorithm>
using namespace std;
#define ll long long
#define modd 1000000007

int n,m,sheng[2000001],jiang[2000001],maxx,ans;

int mi(int u,int v)
{
	int now=1;u%=modd;
	while(v)
	{
		if(v&1) now=(ll)now*u%modd;
		u=(ll)u*u%modd;v>>=1;
	}
	return now;
}

int c(int n,int m)
{
	if(m>n) return 0;
	return ((ll)sheng[n]*jiang[n-m]%modd)*jiang[m]%modd;
}

int he(int u,int v)
{
	return (ll)(u+v)%modd;
}

int main()
{
	scanf("%d%d",&n,&m);
	if(!((ll)n*m))
	{
		puts("0");return 0;
	}
	m=min(n,m);maxx=(n-m)/2+m-1;
	sheng[0]=jiang[0]=ans=1;
	for(int i=1;i<=maxx;i++) sheng[i]=(ll)sheng[i-1]*i%modd;
	jiang[maxx]=mi(sheng[maxx],modd-2);
	for(int i=maxx-1;i;i--) jiang[i]=(ll)jiang[i+1]*(i+1)%modd;
	for(int i=2;i<=m;i++) ans=he(ans,he(n-i-1>=0 ? c((n-i-1)/2+i-1,i-1):0,c((n-i)/2+i-1,i-1)));
	printf("%d\n",ans);
	return 0;
}