//泥の分際で私だけの大切を奪おうだなん #include<bits/stdc++.h> using namespace std; inline int read(){ int s=0,w=1; char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();} while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar(); return s*w; } const int p=1e9+7; int qp(int x,int y) { int res=1; for(int t=x; y; y>>=1,t=1ll*t*t%p) if(y&1) res=1ll*res*t%p; return res; } int a[3003][3003]; bool f[3003][3003]; signed main() { int n=read(),m=read(); for(int i=1; i<=m; ++i) for(int j=1; j<=n; ++j) a[i][j]=read(); queue<pair<int,int>> q; int ans=0; for(int i=1; i<=n; ++i) for(int j=1; j<=n; ++j) if(i!=j) if(!f[i][j]) { f[i][j]=1;q.push({i,j}); int A=0,B=0; while(!q.empty()) { auto [x,y]=q.front();q.pop(); if(x>y) ++A; ++B; for(int k=1; k<=m; ++k) if(!f[a[k][x]][a[k][y]]) f[a[k][x]][a[k][y]]=1, q.push({a[k][x],a[k][y]}); } ans=(ans+1ll*A*(B-A)%p*qp(B,p-2))%p; } printf("%d\n",ans); return 0; }
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 | //泥の分際で私だけの大切を奪おうだなん #include<bits/stdc++.h> using namespace std; inline int read(){ int s=0,w=1; char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();} while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar(); return s*w; } const int p=1e9+7; int qp(int x,int y) { int res=1; for(int t=x; y; y>>=1,t=1ll*t*t%p) if(y&1) res=1ll*res*t%p; return res; } int a[3003][3003]; bool f[3003][3003]; signed main() { int n=read(),m=read(); for(int i=1; i<=m; ++i) for(int j=1; j<=n; ++j) a[i][j]=read(); queue<pair<int,int>> q; int ans=0; for(int i=1; i<=n; ++i) for(int j=1; j<=n; ++j) if(i!=j) if(!f[i][j]) { f[i][j]=1;q.push({i,j}); int A=0,B=0; while(!q.empty()) { auto [x,y]=q.front();q.pop(); if(x>y) ++A; ++B; for(int k=1; k<=m; ++k) if(!f[a[k][x]][a[k][y]]) f[a[k][x]][a[k][y]]=1, q.push({a[k][x],a[k][y]}); } ans=(ans+1ll*A*(B-A)%p*qp(B,p-2))%p; } printf("%d\n",ans); return 0; } |