东林oj 1541 数字三角形加强版

题目链接 https://acm.nefu.edu.cn/JudgeOnline/problemShow.php?problem_id=1541
本题用到：杨辉三角，利用逆元求组合数，快速幂。
刚开始，想到的是矩阵连乘，按照题目要求写了6行图形，发现其实是杨辉三角加快速幂。
三角形中的数设为xA^x1B^x2 ,那么x在三角形中就是杨辉三角。F[i][j]=A∗F[i−1][j]+B∗F[i−1][j−1] 式子中除去A和B，与杨辉三角中每个数等于它上方两数之和是相似。
并且第n行每个数A和B的指数和为n-1，第一个为A^n-1B⁰,然后A的指数减一，B的指数加一，直到A⁰B^n-1(找规律得到的，不会推）。
综上所述，第n行第m个数的值为C(n-1,m-1)A^n-mB^m-1.
先打表求出范围内的阶乘，再求阶乘的逆元，0的阶乘的逆元学长的模板里没有（学长应该用了特判），当时样例输出为0，感觉慌得一批，幸亏我机智，想到0的阶乘为1，逆元也为1。（太机智了）
然后输入a b n m，利用c(n-1,m-1)=(n-1)!/((n-m)!*(m-1)!)。求出c(n-1,m-1)，再利用快速幂求出)A^n-mB^m-1。相乘取模结束。

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define N 100009
#define mod 1000000009
using namespace std;
long long fac[N],inv[N];
long long quickmod(long long n,long long m)
{
    long long b=1;
    while(m>0)
    {
        if(m&1)
            b=b*n%mod;
        n=n*n%mod;
        m=m>>1;
    }
    return b;
}
void  ints()
{
      fac[0]=1;
    for(int i=1;i<=N;i++)
    {
         fac[i]=fac[i-1]*i%mod;///求组合式
        // printf("%lld\n",fac[i]);
    }
    inv[0]=1;///学长的模板里没有  要注意
    inv[N]=quickmod(fac[N],mod-2);
    for(int i=N-1;i>=1;i--)
        inv[i]=inv[i+1]*(i+1)%mod;///求组合数的逆元
}
int main()
{
      long long t,x,x1,x2,a,b,n,m,sum;
      int i;
      memset(fac,0,sizeof(fac));
      memset(inv,0,sizeof(inv));
      ints();
      //for(i=1;i<25;i++)
       //printf("%lld\n",inv[i]);
      scanf("%lld",&t);
      while(t--)
      {
          scanf("%lld%lld%lld%lld",&a,&b,&n,&m);
          x=(fac[n-1]*inv[m-1])%mod;
          //printf("%lld\n",x);
          x=x*inv[n-m]%mod;
           //printf("%lld\n",x);
          x1=quickmod(a,n-m);
          x2=quickmod(b,m-1);
          sum=x*x1%mod;
          sum=sum*x2%mod;
          printf("%lld\n",sum);
      }
    return 0;