#defineLLlonglong#definenmax101#definennum20090126LLLLnum[nmax][nmax],fac[nmax];voidinit(){i" />

黄色网页视频 I 影音先锋日日狠狠久久 I 秋霞午夜毛片 I 秋霞一二三区 I 国产成人片无码视频 I 国产 精品 自在自线 I av免费观看网站 I 日本精品久久久久中文字幕5 I 91看视频 I 看全色黄大色黄女片18 I 精品不卡一区 I 亚洲最新精品 I 欧美 激情 在线 I 人妻少妇精品久久 I 国产99视频精品免费专区 I 欧美影院 I 欧美精品在欧美一区二区少妇 I av大片网站 I 国产精品黄色片 I 888久久 I 狠狠干最新 I 看看黄色一级片 I 黄色精品久久 I 三级av在线 I 69色综合 I 国产日韩欧美91 I 亚洲精品偷拍 I 激情小说亚洲图片 I 久久国产视频精品 I 国产综合精品一区二区三区 I 色婷婷国产 I 最新成人av在线 I 国产私拍精品 I 日韩成人影音 I 日日夜夜天天综合

hdu 2643 Rank hdu 2512 一卡通大冒險 stirling

系統 2170 0

hdu 2643
/*
?第二類Stirling數是把包含n個元素的集合劃分為正好k個非空子集的方法的數目。
?遞推公式為:
?S(n,k) = 0(n<k||k=0),
?S(n,n) = S(n,1) = 1,
?S(n,k) = S(n-1,k-1) + kS(n-1,k).
?*/

      #include<stdio.h>

#define LL long long

#define nmax 101

#define nnum 20090126LL

LL num[nmax][nmax], fac[nmax];

void init() {

	int i, j;

	for (i = 1, fac[0] = 1; i < nmax; i++) {

		fac[i] = fac[i - 1] * i % nnum;

	}

	for (i = 1; i < nmax; i++) {

		num[i][1] = 1;

		num[i][0] = 0;

	}

	for (i = 2; i < nmax; i++) {

		for (j = 1; j < nmax; j++) {

			if (i == j) {

				num[i][i] = 1;

			} else {

				num[i][j] = (num[i - 1][j - 1] + num[i - 1][j] * j) % nnum;

			}

		}

	}

}

int main() {

#ifndef ONLINE_JUDGE

	freopen("in.data", "r", stdin);

#endif

	init();

	int T, N, i;

	LL res;

	while (scanf("%d", &T) != EOF) {

		while (T--) {

			scanf("%d", &N);

			for (i = 1, res = 0; i <= N; i++) {

				res += num[N][i] * fac[i];

				res %= nnum;

			}

			printf("%I64d\n", res);

		}

	}



	return 0;

}


    

?

hdu 2512

      #include<stdio.h>

#define nmax 2001

#define nnum 1000

int num[nmax][nmax];

void init() {

	int i, j;

	for (i = 1; i < nmax; i++) {

		num[i][0] = 0, num[i][1] = 1;

	}

	for (i = 2; i < nmax; i++) {

		for (j = 1; j < nmax; j++) {

			if (i == j) {

				num[i][i] = 1;

				continue;

			}

			num[i][j] = (num[i - 1][j - 1] + num[i - 1][j] * j) % nnum;

		}

	}

}

int main() {

#ifndef ONLINE_JUDGE

	freopen("in.data", "r", stdin);

#endif

	int n, x, i, res;

	init();

	while (scanf("%d", &n) != EOF) {

		while (n--) {

			scanf("%d", &x);

			for (i = 1, res = 0; i <= x; i++) {

				res += num[x][i];

				res %= nnum;

			}

			printf("%d\n", res);

		}

	}

	return 0;

}


    
小知識:

Bell數,又稱為貝爾數。
是以埃里克·坦普爾·貝爾(Eric Temple Bell)為名的。

B(n)是包含n個元素的集合的劃分方法的數目。

B(0) = 1, B(1) = 1, B(2) = 2, B(3) = 5,
B(4) = 15, B(5) = 52, B(6) = 203,...

遞推公式為,
B(0) = 1,
B(n+1) = Sum(0,n) C(n,k)B(k). n = 1,2,...

其中,Sum(0,n)表示對k從0到n求和,C(n,k) = n!/[k!(n-k)!]

-------------------------
Stirling數,又稱為斯特靈數。
在組合數學,Stirling數可指兩類數,都是由18世紀數學家James Stirling提出的。

第一類Stirling數是有正負的,其絕對值是包含n個元素的集合分作k個環排列的方法數目。

遞推公式為,
S(n,0) = 0, S(1,1) = 1.
S(n+1,k) = S(n,k-1) + nS(n,k)。

第二類Stirling數是把包含n個元素的集合劃分為正好k個非空子集的方法的數目。

遞推公式為,
S(n,n) = S(n,1) = 1,
S(n,k) = S(n-1,k-1) + kS(n-1,k).
將n個有區別的球的球放入k個無標號的盒子中( n>=k>=1,且盒子不允許為空)的方案數就是stirling數.(即含 n 個元素的集合劃分為 k 個集合的情況數)

  遞推公式:

  S(n,k) = 0 (k > n)

  S(n,1) = 1 (k = 1)

  s(n,k)=1 (n=k)

  S(n,k) = S(n-1,k-1)+k*S(n-1,k) (n >= k >= 2)

  分析:設有n個不同的球,分別用b1,b2,...,bn表示。從中取出一個球bn,bn的放法有以下兩種:

  1.bn獨占一個盒子,那么剩下的球只能放在k-1個盒子里,方案數為S(n-1,k-1);

  2.bn與別的球共占一個盒子,那么可以將b1,b2,...,bn-1這n-1個球放入k個盒子里,然后將bn放入其中一個盒子中,方案數為k*S(n-1,m).


-------------
bell數和stirling數的關系為,

每個貝爾數都是"第二類Stirling數"的和。

B(n) = Sum(1,n) S(n,k).

參考資料: http://planetmath.org/?op=getobj&from=objects&id=9059

hdu 2643 Rank hdu 2512 一卡通大冒險 stirling (斯特靈數)的應用


更多文章、技術交流、商務合作、聯系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦!!!

發表我的評論
最新評論 總共0條評論