1. 程式人生 > >POJ3761 Bubble Sort

POJ3761 Bubble Sort

art -a span turn %d def poj 當前 bubble

對1~n組成的序列進行冒泡排序,一共進行了k趟,問有幾個符合題意的序列。

註意:這裏指每一趟是指交換當前相鄰的全部逆序對,比如:2 1 4 3進行一趟交換就是1 2 3 4

假設我們細心觀察。就會發現。須要進行的趟數等於序列中對於某個最多逆序對數的數。

比如:在序列 3 2 1 4中。3的逆序對為0,2的逆序對為1,1的逆序對為2,4的逆序對為0,一共排序2次。

直接計算最大逆序對對數為k對的排列數不方便,這裏有一個很巧妙的方法,g[k]表示最大逆序對對數小於等於k對的排列總數,那麽終於的答案就是g[k] - g[k-1]。

對於前k個最大的數,他們的最大逆序對對數不會超過k個,因此對於前k大的數有k!種排列,則對於第k大之後的數,每一個數的位置有k+1種,那麽總共就有 k! * (k+1)^(n-k)種情況。

#include<cstdio>
#define tsy 20100713
using namespace std;
typedef long long LL;
LL jc[1000005];
void init()
{
	jc[0] = 1;
	jc[1] = 1;
	for(int i = 2; i <= 1000000; i++)
	{
		jc[i] = jc[i-1]*(i%tsy)%tsy;
	}
}
LL n;
LL ksm(LL x,LL k)
{
	LL res = 1;
	while(k)
	{
		if(k&1) res = (res*x)%tsy;
		k >>= 1;
		x = (x*x)%tsy;
	}
	return res;
}
LL getans(LL k)
{
	return (jc[k]*ksm(k+1,n-k))%tsy;
}
int main()
{
	init();
	int T;
	LL k;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%I64d%I64d",&n,&k);
		printf("%I64d\n",(getans(k) - getans(k-1) + tsy)%tsy);
	}
}


POJ3761 Bubble Sort