基准时间限制:2 秒 空间限制:131072 KB 分值: 80 难度:5级算法题
 收藏
 关注
在一个排列中,如果一对数的前后位置与大小顺序相反,即前面的数大于后面的数,那么它们就称为一个逆序。一个排列中逆序的总数就称为这个排列的逆序数。
如2 4 3 1中,2 1,4 3,4 1,3 1是逆序,逆序数是4。
 
1-n的全排列中,逆序数最小为0(正序),最大为n*(n-1) / 2(倒序)
给出2个数n和k,求1-n的全排列中,逆序数为k的排列有多少种?
例如:n = 4 k = 3。
 
1 2 3 4的排列中逆序为3的共有6个,分别是:
1 4 3 2
2 3 4 1
2 4 1 3
3 1 4 2
3 2 1 4
4 1 2 3
 
由于逆序排列的数量非常大,因此只需计算并输出该数 Mod 10^9 + 7的结果就可以了。
 
Input
第1行:一个数T,表示后面用作输入测试的数的数量。(1 <= T <= 10000)
第2 - T + 1行:每行2个数n,k。中间用空格分隔。(2 <= n <= 1000, 0 <= k <= 20000)
Output
共T行,对应逆序排列的数量 Mod (10^9 + 7)
Input示例
1
4 3
Output示例
6

设f(n,k)表示n个数的排列中逆序数个数为k的排列数。

最大的数n可能会排在第n-i位,从而产生i个与n有关的逆序对,去掉n之后,剩下的n-1个数的排列有k-i个逆序对。所以,f(n,k)=求和(f(n-1,k-i))(0<=i<n)。
同理有f(n,k-1)=求和(f(n-1,k-1-i))(0<=i<n)。
两式相减,可得f(n,k)-f(n,k-1)=f(n-1,k)-f(n-1,k-n)。
递推公式为f(n,k)=f(n,k-1)+f(n-1,k)-f(n-1,k-n)。
然后动态规划可得。
#include<cstdio>
using namespace std; typedef long long ll;
const int MAXK = 2e4+;
const int MAXN = 1e3+;
const int mod = 1e9+;
#define min(a,b) (a<b)?a:b int n,k,dp[MAXN][MAXK]; // dp[n,k] = dp[n,k-1] + dp[n-1,k] - dp[n-1,k-n]; int getMod(ll t) {
if(t >= mod) return t-mod;
if(t<) return t+mod;
return t;
} void init() {
int i,j;
for(i=;i<=;i++) {
dp[i][]=;
for(j=;j<=i*(i-)/&&j<=;j++) {
ll tmp=;
ll tmp1=dp[i][j-];
ll tmp2=dp[i-][j];
ll tmp3=(j>=i)?dp[i-][j-i]:;
tmp = tmp1+tmp2-tmp3;
dp[i][j] = getMod(tmp);
}
}
}
int main () {
init();
int T; scanf("%d",&T);
while (T--){
scanf("%d %d",&n,&k);
printf("%d\n",dp[n][k]);
}
return ;
}

最新文章

  1. 如何在网页标题栏加入logo图标?
  2. 【python】如何在某.py文件中调用其他.py内的函数
  3. Win10怎么设置始终以管理员身份运行应用程序
  4. Jenkins+svn+maven
  5. 类中用const限定的成员函数
  6. 开发环境安装 Java Mysql MyEclipse Android Adt
  7. 排序算法_HeapSort
  8. NDK开发之访问域
  9. s​s​h​配​置​公​钥​和​私​钥​登​陆​S​e​c​u​r​e​C​R​T
  10. Java阶段性测试--知识点:数组,面向对象,集合、线程,IO流
  11. 一起学习c++11——c++11中的新语法
  12. MySQL保留字不能作为字段名使用
  13. 三十二、Linux 进程与信号——不可靠信号
  14. python 多线程中子线程和主线程相互通信
  15. linux matlab2016 安装
  16. java-信息安全(十五)-单向认证
  17. FFmpeg在JAVA中的使用-音频提取&amp;字幕压缩
  18. CentOS 7.4 系统安装 git
  19. SQLServer2005重建索引前后对比【转】
  20. oracle 一对多数据分页查询筛选

热门文章

  1. Elasticsearch教程-从入门到精通(转)
  2. python2.X编码
  3. 高并发秒杀系统方案(分布式session)
  4. HTML状态消息
  5. hdu2159FATE(二维背包)
  6. Vue.Js添加自定义插件
  7. 012-centos6.5配置静态ip
  8. CTC(Connectionist Temporal Classification)介绍
  9. byte、二进制、十进制数值之间的转换
  10. 由于防火墙限制无法访问linux服务器上的tomcat应用