神TM有是結論題，我討厭結論題mmp。

楊氏矩陣了解一下(建議去維基百科)。

反正就是推柿子，使勁推，最後寫起來有一點小麻煩，但是在草稿紙(然鵝我木有啊)上思路清晰的話還是沒問題的。

#include<cstdio>
#include<cctype>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=2000000,ha=1e9+7;

inline int read(){
	int x=0; char ch=getchar();
	for(;!isdigit(ch);ch=getchar());
	for(;isdigit(ch);ch=getchar()) x=x*10+ch-‘0‘;
	return x;
}

inline int add(int x,int y){ x+=y; return x>=ha?x-ha:x;}
inline void ADD(int &x,int y){ x+=y; if(x>=ha) x-=ha;}

inline int ksm(int x,int y){
	int an=1;
	for(;y;y>>=1,x=x*(ll)x%ha) if(y&1) an=an*(ll)x%ha;
	return an;
}

int jc[maxn+5],qz[maxn+5],hz[maxn+5];
int A,B,C,D,T,ans,MX,MY;

inline void init(){
	jc[0]=qz[0]=hz[0]=1;
	for(int i=1;i<=maxn;i++){
		jc[i]=jc[i-1]*(ll)i%ha;
		qz[i]=qz[i-1]*(ll)ksm(i,i)%ha;
	}
	
	for(int i=1;i<=maxn;i++) hz[i]=hz[i-1]*(ll)jc[i]%ha;
}

inline int gethz(int x,int len){
	if(!len) return 1;
    return hz[x+len-1]*(ll)ksm(hz[x-1]*(ll)ksm(jc[x-1],len)%ha,ha-2)%ha;
}

inline int Get(int x,int y,int derta){
	if(x>y) swap(x,y);

	if(!x) return 1;
	
	int now=qz[x+derta]*(ll)ksm(qz[derta]*(ll)ksm(jc[x+derta]*(ll)ksm(jc[derta],ha-2)%ha,derta)%ha,ha-2)%ha;
	now=now*(ll)ksm(jc[y+derta]*(ll)ksm(jc[x+derta],ha-2)%ha,x)%ha*(ll)gethz(y+derta+1,x-1)%ha;
	return now;
}

int main(){
	init();
	
	T=read();
	while(T--){
		A=read(),B=read(),C=read(),D=read();
		A-=B-1,C-=D-1,MX=min(A,C),MY=min(B,D);
		
		if((A>C&&B<D)||(A<C&&B>D))
		    ans=Get(A+C-MX*2,MY,0)*(ll)Get(MX,B+D-MY*2,0)%ha*(ll)Get(MX,MY,A+C+B+D-MX*2-MY*2)%ha;
        else ans=Get(A+C-MX,B+D-MY,0);
        
		ans=jc[A*(ll)B+C*(ll)D-MX*(ll)MY]*(ll)ksm(ans,ha-2)%ha;
		
		printf("%d\n",ans);
	}
	
	return 0;
}

SD 一輪集訓 day1 lose