标签:space ios size blank += cst 滚动 namespace algorithm
题目:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1250
看了半天...
把第一问想成逆序对的话似乎很容易想了,新加入一个数,可以往前挪动,增加的逆序对数就是它后面那些数的个数;
所以 f[i][j] = ∑(k = max( 0 , j - i + 1)) f[i-1][k],用前缀和即可;
第二问正好用第一类斯特林数;
第一类斯特林数 str[i][j] 表示把 i 个数分成 j 个环,环有顺序的方案数,str[i][j] = str[i-1][j-1] (自成一环) + ( i - 1 ) * str[i-1][j] (跟到某个后面)
对应到这道题,原来每个数自己是一个环,表示它就在自己的位置上;
一个有顺序的环,按一个顺序往过数,一个数到下一个数相互交换,就对应了一种交换;
str[i][j] 中,i - j 个数与别的环合并,也就是发生了一次交换,所以答案就是 ∑(n - k <= i <= n ) str[n][i]
注意不要 MLE ...
用滚动数组求斯特林数得注意一点,初值 str[0][0] 不要赋成1,然后后面 i == j 时也求出来即可,感性理解一下的话~
这题...好像也不难啊。
代码如下:
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; typedef long long ll; int const maxn=3005,mod=1000000007; int n,m; ll f[5][maxn],s[5][maxn],str[5][maxn],ans1,ans2; void init() { // for(int i=0;i<=1;i++)str[i][i]=1; str[1][1]=1;//滚动数组!让初始化符合意义 for(int i=2;i<=n;i++) for(int j=1;j<=i;j++)//因为滚动所以 str[i][i]也在这里求 { int p=(i&1),q=!p; str[p][j]=(str[q][j-1]+(i-1)*str[q][j]%mod)%mod; } } int main() { scanf("%d%d",&n,&m); init(); for(int i=0;i<=1;i++) { f[i][0]=1,s[i][0]=1; for(int j=1;j<=m;j++)s[i][j]+=s[i][j-1]; } for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) { int p=(i&1),q=!p; if(j-i+1<=0)f[p][j]=s[q][j]; else f[p][j]=(s[q][j]-s[q][j-i]+mod)%mod; if(j)s[p][j]=(s[p][j-1]+f[p][j])%mod; else s[p][j]=f[p][j]; } int tmp=m; while(tmp>=0)ans1=(ans1+f[n&1][tmp])%mod,tmp-=2; for(int i=n-m;i<=n;i++)ans2=(ans2+str[n&1][i])%mod; printf("%lld %lld\n",ans1,ans2); return 0; }
标签:space ios size blank += cst 滚动 namespace algorithm
原文地址:https://www.cnblogs.com/Zinn/p/9606551.html