bzoj2729: [HNOI2012]排队(高精度+组合数学)

题目传送门

解法:
本来我想的是这样的。。。
n个男生插老师。
然后插女生。
然后怎么算都出不来12。。。
%题解。。
原来还有一种叫做女生把老师隔开的东西。

分情况讨论:
老师之间有男生。
那么就是n个男生有n+1个间隔来插老师。
然后一共有n+2个人了。那么又有n+3个间隔来插女生。
排列组合即可。

老师之间绑一个女生。
那么这个小团体可以看作一个男生因为他两边随便放什么人都行。
这个团体有2*m种可能。
然后一共n+1个男生。
还剩m-1个女生。
此时剩下n+2个间隔。
又是排列组合。

代码实现:

#include
#include
#include
#include
#include
#include
#include
using namespace std;
typedef long long ll;
struct node {
    int len;ll a[2100];
    node() {len=0;memset(a,0,sizeof(a));}
}per;const ll mod=1000000000;
node chengfa(node n1,ll x) {
    node no;no.len=n1.len;
    for(int i=1;i<=no.len;i++)no.a[i]=n1.a[i]*x;
    for(int i=1;i<=no.len;i++) {no.a[i+1]+=no.a[i]/mod;no.a[i]%=mod;}
    int i=no.len;
    while(no.a[i+1]>0) {
        i++;no.a[i+1]+=no.a[i]/mod;no.a[i]%=mod;
    }no.len=i;return no;
}
node jiafa(node n1,node n2) {
    node no;no.len=max(n1.len,n2.len);
    for(int i=1;i<=no.len;i++)no.a[i]=n1.a[i]+n2.a[i];
    for(int i=1;i<=no.len;i++) {no.a[i+1]+=no.a[i]/mod;no.a[i]%=mod;}
    int i=no.len;
    while(no.a[i+1]>0) {
        i++;no.a[i+1]+=no.a[i]/mod;no.a[i]%=mod;
    }no.len=i;return no;
}
int main() {
    int n,m;scanf("%d%d",&n,&m);per.len=1;per.a[1]=1ll;node ans1,ans2;
    if(n+3<m){printf("0\n");return 0;}
    ans1=per;for(int i=1;i<=2;i++)ans1=chengfa(ans1,ll(n-i+2));
    for(int i=1;i<=m;i++)ans1=chengfa(ans1,ll(n+4-i));
    for(int i=1;i<=n;i++)ans1=chengfa(ans1,ll(i));
    ans2=per;ans2=chengfa(ans2,ll(2));ans2=chengfa(ans2,ll(m));
    for(int i=1;i<m;i++)ans2=chengfa(ans2,ll(n+3-i));
    for(int i=1;i<=n+1;i++)ans2=chengfa(ans2,ll(i));
    node ans=jiafa(ans1,ans2);
    for(int i=ans.len;i>=1;i--)printf("%lld",ans.a[i]);printf("\n");
    return 0;
}

你可能感兴趣的:(BZOJ,杂,数论)