2227: [Zjoi2011]看电影(movie)

来源：互联网发布：兰芝睡眠面膜知乎编辑：程序博客网时间：2024/06/06 18:30

经典思路：线转化为圈，利用对称性
假设有k+1个椅子围成一圈，每张票上写着1~k+1（而不是1~k)，因为概率均等，我们统计有多少种方案导致椅子k+1上坐着人
一共有(k+1)n种可能的情况。为了方便叙述，我们假定每种情况模拟一遍，每次在没人坐的椅子上放一枚硬币（每次恰好有k+1-n把椅子没人坐），则所有情况都模拟完之后每把椅子上都有(k+1-n)(k+1)n-1枚硬币，也就是说，有(k+1-n)(k+1)n-1种情况下，椅子k+1上没人坐。
改回每张票写1~k之后，这个结论不变。——————-from liurujia

#include<bits/stdc++.h>#define rep(i,k,n) for(int i=k;i<=n;i++)#define rep2(i,k,n) for(int i=k;i>=n;i--)using namespace std;const int N=1000005;int gcd(int a,int b){return !b ? a : gcd(b,a%b);}struct gao{    int num[N],len;    void init(){len=1;memset(num,0,sizeof(num));}    void put(){        printf("%d",num[len]);        rep2(i,len-1,1)printf("%04d",num[i]);    }    void operator*=(int x){        rep(i,1,len)num[i]*=x;        for(int i=1;i<=len || num[i];i++)num[i+1]+=num[i]/10000,num[i]%=10000,len=max(len,i);    }}A,B;int T,k,n,a[N],b[N];int main(){    scanf("%d",&T);    while(T--){        scanf("%d%d",&n,&k);        if(n>k){printf("0 1\n");continue;}        rep(i,1,n)a[i]=k+1,b[i]=k;        a[n]=(k+1-n);        rep(i,1,n){int d=gcd(a[n],b[i]);a[n]/=d,b[i]/=d;}        A.init(),B.init();        A.num[1]=1,B.num[1]=1;        rep(i,1,n){            A*=a[i];            B*=b[i];        }        A.put();printf(" ");B.put();puts("");    }}

0 0