組合數(shù)的計(jì)算方法

組合數(shù):從n個(gè)不同元素中取出m(m≤n)個(gè)元素的所有組合的個(gè)數(shù),叫做從n個(gè)不同元素中取出m個(gè)元素的組合數(shù)。計(jì)算公式為:
C(n,m)=n!/((n-m)!\times m!), \text{where}\ m\leq n

  • 性質(zhì)1:C(n,m)= C(n,n-m)
  • 性質(zhì)2:C(n,m)=C(n-1,m-1)+C(n-1,m)

第一種方法:打表

根據(jù)性質(zhì)2直接構(gòu)建一個(gè)n\times n的矩陣進(jìn)行計(jì)算:

public class Template {
    static int mod = (int) 1e9 + 7;
    static int max = 110;
    static long[][] com = new long[max][max];

    public static void main(String[] args) {
        int n = 100, m = 30;
        for (int i = 0; i < max; i++) {
            com[i][0] = com[i][i] = 1;
            for (int j = 1; j < i; j++) {
                com[i][j] = (com[i - 1][j - 1] + com[i - 1][j]) % mod;
            }
        }
        System.out.println(com[n][m]);
    }
}

空間復(fù)雜度:O(n^2)

預(yù)處理時(shí)間復(fù)雜度:O(n^2),查詢時(shí)間復(fù)雜度:O(1)

第二種方法:階乘無模

根據(jù)組合的組合數(shù)的計(jì)算公式C(n,m)=n!/((n-m)!\times m!)進(jìn)行:

public class Template {
    static int mod = (int) 1e9 + 7;
    static int max = 110;
    static long[] fac = new long[max];

    public static void main(String[] args) {
        int n = 20, m = 10;
        fac[0] = 1;
        for (int i = 1; i < max; i++) {
            fac[i] = (fac[i - 1] * i);
        }
        System.out.println(fac[n] / fac[m] / fac[n - m]);
    }
}

空間復(fù)雜度:O(n)

預(yù)處理時(shí)間復(fù)雜度:O(1)濒旦,查詢時(shí)間復(fù)雜度:O(n)

由于涉及除法误墓,無法直接取模,所以引入乘法逆元出吹。

第三種方法:乘法逆元

逆元:對(duì)于apap互素),若a*b\%p\equiv1辙喂,則稱b的最小正整數(shù)解為a%p的逆元捶牢。

當(dāng)求解(a/b)\%p,如果知道b\%p的逆元為c巍耗,那么可以轉(zhuǎn)化為(a/b)\%p=a*c\%p=(a\%p)(c\%p)\%p秋麸。暴力做法:

public class Template {
    static int mod = (int) 1e9 + 7;
    static int max = 110;
    static long[] fac = new long[max];
    static long[] inv = new long[max];

    public static void main(String[] args) {
        int n = 100, m = 30;
        inv[0] = fac[0] = 1;
        for (int i = 1; i < max; i++) {
            fac[i] = (fac[i - 1] * i) % mod;
            inv[i] = inv(fac[i]);
        }
        System.out.println(((fac[n] * inv[m]) % mod * inv[n - m]) % mod);
    }

    public static long inv(long a) {
        for (int x = 1; x < mod; x++) {
            if (a * x % mod == 1) return x;
        }
        return 0;
    }
}

空間復(fù)雜度:O(n)

預(yù)處理時(shí)間復(fù)雜度:O(p),其中p=mod炬太,查詢時(shí)間復(fù)雜度:O(1)

第四種方法:乘法逆元+快速冪+階乘

費(fèi)馬小定理:對(duì)于a和素?cái)?shù)p灸蟆,滿足a^{p-1}\%p\equiv 1

因?yàn)?img class="math-inline" src="https://math.jianshu.com/math?formula=a%5E%7Bp-1%7D%3Da%5E%7Bp-2%7D*a" alt="a^{p-1}=a^{p-2}*a" mathimg="1">亲族,所以有a^{p-2}*a\%p\equiv 1炒考。根據(jù)逆元的定義可知,a^{p-2}a的逆元霎迫。因此可以將求解逆元的問題轉(zhuǎn)換為a^{p-2}的快速冪問題票腰。

public class Template {
    static int mod = (int) 1e9 + 7;
    static int max = 110;
    static long[] fac = new long[max];
    static long[] inv = new long[max];

    public static void main(String[] args) {
        int n = 100, m = 30;
        inv[0] = fac[0] = 1;
        for (int i = 1; i < max; i++) {
            fac[i] = (fac[i - 1] * i) % mod;
            inv[i] = inv(fac[i]);
        }
        System.out.println(((fac[n] * inv[m]) % mod * inv[n - m]) % mod);
    }

    public static long pow(long a, long b) {
        long ans = 1;
        while (b > 0) {
            if ((b & 1) == 1) ans = (ans * a) % mod;
            a = a * a % mod;
            b = b >> 1;
        }
        return ans;
    }

    public static long inv(long a) {
        return pow(a, mod - 2);
    }
}

空間復(fù)雜度:O(n)

預(yù)處理時(shí)間復(fù)雜度:O(\log p),其中p=mod女气,查詢時(shí)間復(fù)雜度為O(1)

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末杏慰,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子炼鞠,更是在濱河造成了極大的恐慌缘滥,老刑警劉巖,帶你破解...
    沈念sama閱讀 218,386評(píng)論 6 506
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件谒主,死亡現(xiàn)場(chǎng)離奇詭異朝扼,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)霎肯,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,142評(píng)論 3 394
  • 文/潘曉璐 我一進(jìn)店門擎颖,熙熙樓的掌柜王于貴愁眉苦臉地迎上來榛斯,“玉大人,你說我怎么就攤上這事搂捧⊥运祝” “怎么了?”我有些...
    開封第一講書人閱讀 164,704評(píng)論 0 353
  • 文/不壞的土叔 我叫張陵允跑,是天一觀的道長(zhǎng)王凑。 經(jīng)常有香客問我,道長(zhǎng)聋丝,這世上最難降的妖魔是什么索烹? 我笑而不...
    開封第一講書人閱讀 58,702評(píng)論 1 294
  • 正文 為了忘掉前任,我火速辦了婚禮弱睦,結(jié)果婚禮上百姓,老公的妹妹穿的比我還像新娘。我一直安慰自己况木,他們只是感情好瓣戚,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,716評(píng)論 6 392
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著焦读,像睡著了一般。 火紅的嫁衣襯著肌膚如雪舱权。 梳的紋絲不亂的頭發(fā)上矗晃,一...
    開封第一講書人閱讀 51,573評(píng)論 1 305
  • 那天,我揣著相機(jī)與錄音宴倍,去河邊找鬼张症。 笑死,一個(gè)胖子當(dāng)著我的面吹牛鸵贬,可吹牛的內(nèi)容都是我干的俗他。 我是一名探鬼主播,決...
    沈念sama閱讀 40,314評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼阔逼,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼兆衅!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起嗜浮,我...
    開封第一講書人閱讀 39,230評(píng)論 0 276
  • 序言:老撾萬榮一對(duì)情侶失蹤羡亩,失蹤者是張志新(化名)和其女友劉穎,沒想到半個(gè)月后危融,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體畏铆,經(jīng)...
    沈念sama閱讀 45,680評(píng)論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,873評(píng)論 3 336
  • 正文 我和宋清朗相戀三年吉殃,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了辞居。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片楷怒。...
    茶點(diǎn)故事閱讀 39,991評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖瓦灶,靈堂內(nèi)的尸體忽然破棺而出鸠删,到底是詐尸還是另有隱情,我是刑警寧澤倚搬,帶...
    沈念sama閱讀 35,706評(píng)論 5 346
  • 正文 年R本政府宣布冶共,位于F島的核電站,受9級(jí)特大地震影響每界,放射性物質(zhì)發(fā)生泄漏捅僵。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,329評(píng)論 3 330
  • 文/蒙蒙 一眨层、第九天 我趴在偏房一處隱蔽的房頂上張望庙楚。 院中可真熱鬧,春花似錦趴樱、人聲如沸馒闷。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,910評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽纳账。三九已至,卻和暖如春捺疼,著一層夾襖步出監(jiān)牢的瞬間疏虫,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,038評(píng)論 1 270
  • 我被黑心中介騙來泰國打工啤呼, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留卧秘,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 48,158評(píng)論 3 370
  • 正文 我出身青樓官扣,卻偏偏與公主長(zhǎng)得像翅敌,于是被迫代替她去往敵國和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子惕蹄,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,941評(píng)論 2 355

推薦閱讀更多精彩內(nèi)容