bzoj2729: [HNOI2012]排队

Description

某中学有 n 名男同学，m 名女同学和两名老师要排队参加体检。他们排成一条直线，并且任意两名女同学不能相邻，两名老师也不能相邻，那么一共有多少种排法呢？（注意：任意两个人都是不同的）

Input

只有一行且为用空格隔开的两个非负整数 n 和 m，其含义如上所述。

对于 30%的数据 n≤100,m≤100

对于 100%的数据 n≤2000,m≤2000
Output

输出文件 output.txt 仅包含一个非负整数，表示不同的排法个数。注意答案可能很大。
Sample Input

1 1
Sample Output

12

这题其实并不怎么难。。
又是在站军姿的时候想出来的。。
话说在开会的时候也想了一会，主任说我们班有人不认真，不知道是不是我
反正就想了很多东西
最后想到一种做法，我们可以先不管老师这个条件，然后拿老师看做两个男生，这种情况应该大家都会做啦，大概就是想把男生排好，然后把女生插进去公式不写了，不会的自己看代码。
然后我们现在的情况是不合法的，因为有可能老师可能站在一起了。所以要减去。怎么减呢？我们规定，老师一定要站在一起就好了啊，这样的话就把老师看做一个男生，然后别的和上面一样就好了
高精度又写挫了QAQ
CODE:

#include<cstdio>
#include<cstring>
int n,m;
struct Int
{
    int a[20008];
    int len;//位数 
     void clear() {memset(a,0,sizeof(a));len=1;}
};
Int Mul1 (Int x,int y)
{
    for (int u=1;u<=x.len;u++)
        x.a[u]*=y;
    x.len+=50;
    for (int u=1;u<=x.len;u++)
        if (x.a[u]>=10)
        {
            x.a[u+1]=x.a[u+1]+x.a[u]/10;
            x.a[u]%=10;
        }
    while (x.len>0&&x.a[x.len]==0) x.len--;
    return x;

}
Int get (int x)
{
    Int ooo;ooo.clear();ooo.a[1]=1;
    for (int u=1;u<=x;u++) ooo=Mul1(ooo,u);
    return ooo;
}
Int Mul (Int a,Int b)
{
    Int c;c.clear();
    for (int u=1;u<=a.len;u++)
        for (int i=1;i<=b.len;i++)
            c.a[u+i-1]=c.a[u+i-1]+a.a[u]*b.a[i];
    c.len=0;
    while (c.a[c.len+1]>0||c.len<=a.len+b.len+5)
    {
        c.len++;
        c.a[c.len+1]=c.a[c.len+1]+(c.a[c.len]/10);
        c.a[c.len]%=10;
    }
    while (c.len>0&&c.a[c.len]==0) c.len--;
    return c;
}
Int P(int x,int y)
{
    Int tt;tt.clear();tt.a[1]=1;
    for (int u=1;u<=y;u++) tt=Mul1(tt,x-u+1);
    return tt;
}
Int minus (Int a,Int b)
{
    for (int u=1;u<=a.len;u++)
    {
        a.a[u]-=b.a[u];
        if (a.a[u]<0)
        {
            a.a[u]+=10;
            a.a[u+1]--;
        }
    }
    while (a.len>0&&a.a[a.len]==0) a.len--;
    return a;
}
int main()
{
    scanf("%d%d",&n,&m);
    Int ans;
    ans=Mul(get(n+2),P(n+3,m));//先把老师的条件忽略
/*  for (int u=ans.len;u>=1;u--) printf("%d",ans.a[u]);
    printf("\n");*/
    ans=minus(ans,Mul1(Mul(get(n+1),P(n+2,m)),2));
    if (ans.len<=1&&ans.a[1]==0) printf("0\n");
    else for (int u=ans.len;u>=1;u--) printf("%d",ans.a[u]);
    return 0;
}