【模版】组合数学

#include<bits/stdc++.h> #define ll long long using namespace std; const int mx=2e6; const int mod=1e9+7; int n,m; ll res,fac[mx+5],inv[mx+5]; //相同的字符 x ,取序列中的最后一个 ll fpow(ll x,ll y) { ll mul=1; for(;y;y>>=1) { if(y&1) mul=mul*x%mod; x=x*x%mod; } return mul; } ll C(ll x,ll y) { return fac[x]*inv[x-y]%mod*inv[y]%mod; } ll A(ll x,ll y) { return fac[x]*inv[x-y]%mod; } ll G(ll x,ll y) { return C(x+y-1,y-1); } int main() { // freopen("data.in","r",stdin); fac[0]=1; for(int i=1;i<=mx;i++) fac[i]=fac[i-1]*i%mod; inv[mx]=fpow(fac[mx],mod-2); for(int i=mx;i>=1;i--) { inv[i-1]=inv[i]*i%mod; } }

__EOF__

本文作者仰望星空的蚂蚁
本文链接https://www.cnblogs.com/cqbzly/p/17530320.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   仰望星空的蚂蚁  阅读(4)  评论(0编辑  收藏  举报  
相关博文:
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」
点击右上角即可分享
微信分享提示