矩阵生成树

#include <cstdio>
#define LL long long 
 
  const LL mo=1000003;
 
  LL a[101][101],fac[200001],num[101];
 
  LL qpow(LL bas,int powe){
      LL ret=1;
      for (;powe;bas*=bas,bas%=mo){
        if (powe&1) ret*=bas,ret%=mo;
      powe=powe>>1;    
    }
    return(ret);
  }
 
  LL work(int n){
    for (int i=1;i<=n;i++){
      for (int j=i+1;j<=n;j++){
          LL bas=a[j][i]*qpow(a[i][i],mo-2)%mo;
          for (int k=i;k<=n;k++)
            a[j][k]-=a[i][k]*bas%mo,a[j][k]%=mo,a[j][k]+=mo,a[j][k]%=mo;
      }    
    }
    LL ret=1;
    for (int i=1;i<n;i++) ret*=a[i][i],ret%=mo;
    return(ret);
  }
  
  int main(){
      fac[0]=1;
      for (int i=1;i<=200000;i++) fac[i]=fac[i-1]*i%mo;
      int n;
      while (scanf("%d",&n)!=EOF,n!=0){
         for (int i=1;i<=n;i++)
       for (int j=1;j<=n;j++)
         a[i][j]=0;
      for (int i=1;i<=n;i++){
          scanf("%d",&num[i]);
          for (int j=1;j<=num[i];j++){
            int tt;
            scanf("%d",&tt);
          a[i][tt]--;a[i][i]++;
        }
      }
      for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) a[i][j]%=mo,a[i][j]+=mo,a[i][j]%=mo;
      LL ans=1;
      for (int i=1;i<=n;i++) ans*=fac[num[i]-1],ans%=mo;
      ans*=work(n);ans%=mo;
      ans*=num[1];ans%=mo;
      printf("%lld\n",ans);
    }
  }