BZOJ 4807(車-高精度)

Description
众所周知，車是中国象棋中最厉害的一子之一，它能吃到同一行或同一列中的其他棋子。車跟車显然不能在一起打
起来，于是rly一天又借来了许多许多的車在棋盘上摆了起来……他想知道，在N×M的矩形方格中摆最多个数的車
使其互不吃到的情况下方案数有几种。但是，由于上次摆炮摆得实在太累，他为了偷懒，打算增加一个条件：对于
任何一个車A，如果有其他一个車B在它的上面（車B行号小于車A），那么車A必须在車B的右边（車A列号大于車B）
。
棋子都是相同的。

Input
一行，两个正整数N和M。
N<=1000000，M<=1000000

Output
一行，输出方案数的末尾50位（不足则直接输出）。

Sample Input
2 2
Sample Output
1
HINT

Source

By FancyCoder

答案是$C(max(n,m),min(n,m))$

#include<bits/stdc++.h> 
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,0x3f,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define MEMx(a,b) memset(a,b,sizeof(a));
#define INF (0x3f3f3f3f)
#define F (1000000007)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %lld\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
#define ALL(x) (x).begin(),(x).end()
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define MAXN (1000000+10)
int p[MAXN],tot,mn[MAXN];
bool b[MAXN]={0};
void make_prime(int n)
{
    tot=0;
    Fork(i,2,n)
    {
        if (!b[i]) p[++tot]=i,mn[i]=i;
        For(j,tot)
        {
            if (i*p[j]>n) break;
            b[i*p[j]]=1;mn[i*p[j]]=p[j];
            if (i%p[j]==0) break;  
        }
    }
}
int cnt[MAXN]={0};
void modify(int x,int c) {
    while(x!=1) {
        cnt[mn[x]]+=c;
        x/=mn[x];
    }
}

#define MAXN (1000)
struct BigInteger {
    int n;
    int a[MAXN];
    enum {MOD=10};
    BigInteger() {
        MEM(a) a[0]=1;
    }
    int& operator [] (int p) {return a[p];} 
    const int& operator [] (int p) const {return a[p];}  
    BigInteger(int x) {
        a[0]=1;
        a[1]=x%MOD; 
        x/=MOD;
        while(x) a[++a[0]]=x%MOD,x/=MOD;
    }
    friend BigInteger operator*(BigInteger a,BigInteger b) {
        BigInteger c;
        c[0]=a[0]+b[0];
        For(i,a[0])
            For(j,b[0]){  
                c[i+j-1]+=a[i]*b[j], c[i+j]+=c[i+j-1]/MOD, c[i+j-1]%=MOD;  
            }  
        if (c[c[0]] == 0 ) c[0]--;  
        c[0]=min(c[0],50);
        return c;
    }
    void print() {
        while (a[0]&&a[a[0]] == 0 ) a[0]--;
        ForD(i,a[0]) {
            if(i==a[0]) printf("%d",a[i]);
            else printf("%d",a[i]);
        }
    }
};

int main()
{
//  freopen("bzoj4807.in","r",stdin);
//  freopen(".out","w",stdout);
    int n=read(),m=read();
    if (n<m) swap(n,m);
    make_prime(1000000);
    Fork(i,m+1,n) modify(i,1);
    For(i,n-m) modify(i,-1);
    BigInteger a(1);
    For(i,n)
        For(j,cnt[i])
            a=a*i;
    a.print();
    return 0;
}