百度之星初赛B——Chess

Chess Accepts: 1805
Submissions: 5738
Time Limit: 2000/1000 MS (Java/Others)
Memory Limit: 32768/32768 K (Java/Others)
Problem Description

車是中国象棋中的一种棋子，它能攻击同一行或同一列中没有其他棋子阻隔的棋子。一天，小度在棋盘上摆起了许多車……他想知道，在一共N×M个点的矩形棋盘中摆最多个数的車使其互不攻击的方案数。他经过思考，得出了答案。但他仍不满足，想增加一个条件：对于任何一个車A，如果有其他一个車B在它的上方（車B行号小于車A），那么車A必须在車B的右边（車A列号大于車B）。

现在要问问你，满足要求的方案数是多少。
Input

第一行一个正整数T，表示数据组数。

对于每组数据：一行，两个正整数N和M（N<=1000，M<=1000）。
Output

对于每组数据输出一行，代表方案数模1000000007（1e9+7）。
Sample Input

1
1 1

Sample Output

1

【题解】：

大佬A题超快，跟我们A那些A+B差不多快；

【大佬的代码】：

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const LL mod = 1e9+7;
LL C[1005][1005];
void init()
{
    for(int i = 0; i <= 1000; ++i)
        C[i][0] = 1;
    for(int i = 1; i <= 1000; ++i)
        for(int j = 1; j <= i; ++j)
            C[i][j] = (C[i-1][j]+C[i-1][j-1])%mod;
}
int n, m;
int main()
{
    int re; init();
    scanf("%d", &re);
    while(re--)
    {
        scanf("%d %d", &n, &m);
        int minn = min(n, m);
        int maxx = max(n, m);
        printf("%I64d\n", C[maxx][minn]);
    }
    return 0;
}

【我们队的代码】（其实我都看不懂）

#include<cstdio>
#include<cstring>
#define N 1000005
int n,m,t,pri[N],que[N],cnt[N],ans[500],tot,mor;
long long b=1e9+7;
void calc(int lim,int flg){
    for(int i=2;i<=lim;++i)
        for(int num=i;num!=1;)
            cnt[pri[num]]+=flg,num/=pri[num];
}
int quyu(){
	long long c=0;
    for(int i=tot;i>0;i--)
    {c*=10;c+=ans[i];c%=b;}
    printf("%I64d\n",c);
}
main(){
	int T;
	scanf("%d",&T);
	while(T--){
		memset(ans,0,sizeof(ans));
		memset(pri,0,sizeof(pri));
		memset(que,0,sizeof(que));
		memset(cnt,0,sizeof(cnt));
		tot=1;t=0;
		scanf("%d%d",&n,&m);
	    if(n<m)n^=m^=n^=m;
	    for(int i=2;i<=n;++i){
	        if(!pri[i])que[++t]=pri[i]=i;
	        for(int j=1;j<=t&&i*que[j]<=n;++j){
	            pri[que[j]*i]=que[j];
	            if(!(i%que[j]))break;
	        }
	    }
	    ans[1]=1;
	    calc(n,+1);
	    calc(m,-1);
	    calc(n-m,-1);
	    for(int i=1;i<=t;++i)
	        for(int j=1;j<=cnt[que[i]];++j){tot+=7;
	            for(int k=1;k<=tot;++k)ans[k]*=que[i];
	            for(int k=1;k<=tot;++k)ans[k+1]+=ans[k]/10,ans[k]%=10;
	            while(tot&&!ans[tot])--tot;
	        }
	    while(tot&&!ans[tot])--tot;
	    quyu();
	}
}