排列组合+组合数取模 HDU 5894
1 // 排列组合+组合数取模 HDU 5894 2 // 题意：n个座位不同，m个人去坐(人是一样的)，每个人之间至少相隔k个座位问方案数 3 // 思路： 4 // 定好m个人 相邻人之间k个座位 剩下就剩n-(m+1)*k个座位 5 // 剩下座位去插m个不同的盒子==就等价n个相同的球放m个不同的盒子 6 // 然后组合数出来了 7 // 乘n的话是枚举座位，除m是去掉枚举第一个座位的时候，剩下人相邻的座位相对不变的情况 8 9 #include <iostream> 10 #include <algorithm> 11 #include <cstring> 12 #include <cstdio> 13 #include <vector> 14 #include <cmath> 15 #include <map> 16 #include <queue> 17 using namespace std; 18 #define LL long long 19 typedef pair<int,int> pii; 20 const int inf = 0x3f3f3f3f; 21 const int mod = 1e9+7; 22 const int N = 1e6+10; 23 const int maxx = 200010; 24 #define clc(a,b) memset(a,b,sizeof(a)) 25 const double eps = 1e-8; 26 void fre() {freopen("in.txt","r",stdin);} 27 void freout() {freopen("out.txt","w",stdout);} 28 inline int read() {int x=0,f=1;char ch=getchar();while(ch>'9'||ch<'0') {if(ch=='-') f=-1; ch=getchar();}while(ch>='0'&&ch<='9') {x=x*10+ch-'0';ch=getchar();}return x*f;} 29 30 int f[N]; 31 int inv(int x){ 32 int ret=1,y=mod-2; 33 while(y){ 34 if(y&1)ret=1ll*ret*x%mod; 35 y>>=1;x=1ll*x*x%mod; 36 } 37 return ret; 38 } 39 int C(int n,int m){ 40 if(n<m)return 0; 41 int ret=1ll*f[n]*inv(f[m])%mod; 42 ret=1ll*ret*inv(f[n-m])%mod; 43 return ret; 44 } 45 int lucas(int n,int m){ 46 if(m == 0) return 1; 47 return 1ll*C(n % mod, m % mod) * lucas(n / mod, m / mod) % mod; 48 } 49 50 void init(){ 51 f[0]=1; 52 for(int i=1;i<=N-5;++i)f[i]=1ll*i*f[i-1]%mod; 53 } 54 55 int main(){ 56 init(); 57 int T; 58 scanf("%d",&T); 59 while(T--){ 60 int n,m,k; 61 scanf("%d%d%d",&n,&m,&k); 62 if(m==1){ 63 printf("%d\n",n); 64 continue; 65 } 66 LL res=n-(m+1)*k; 67 if(res<0){ 68 printf("0\n"); 69 continue; 70 } 71 LL tem=lucas(n-m*k-1,m-1); 72 LL ans=((tem)%mod)*n%mod; 73 printf("%d\n",(ans*inv(m))%mod); 74 } 75 return 0; 76 }
修正一下，加一组样列：
n=5,m=1,k=1000
这时候用res的话 必须判m==1
其实可以一起去掉，直接lucas也a
数据太重要了
