今日打卡：洛谷：P1224 [NOI2013] 向量内积

哈哈嗨！我又来了！！！今天是大吉好吧！

#include<iostream>
#include<stdio.h>
#include<ctype.h>
#include<bitset>
#include<random>
#include<stdlib.h>
using namespace std;
inline int read(){
	int x=0,f=0; char ch=getchar();
	while(!isdigit(ch)) f|=(ch==45),ch=getchar();
	while(isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
	return f?-x:x;
}
int n,m,k;
struct matrix{
	vector<vector<int> > a;
	int n,m;
	inline void build(){
		a.resize(n+1);
		for(int i=1;i<=n;++i) a[i].resize(m+1);
	}
	matrix operator * (matrix B){
		matrix A=*this,C;
		C.n=A.n,C.m=B.m;
		C.a.resize(C.n+1);
		for(int i=1;i<=C.n;++i) C.a[i].resize(C.m+1);
		for(int i=1;i<=C.n;++i)
			for(int j=1;j<=C.m;++j){
				C.a[i][j]=0;
				for(int k=1;k<=A.m;++k){
					C.a[i][j]+=A.a[i][k]*B.a[k][j];
				}
				C.a[i][j]%=k;
			}
		return C;
	}
};
namespace sub2{
	inline void main(){
		mt19937 rnd(114514);
		matrix A,B;
		A.n=n,A.m=m,A.build();
		B.n=m,B.m=n,B.build();
		for(int i=1;i<=n;++i){
			for(int j=1;j<=m;++j){
				B.a[j][i]=A.a[i][j]=read()%k;
			}
		}
		for(int fick=1;fick<=10;++fick){
			matrix R;
			R.n=n,R.m=1,R.build();
			int sum=0;
			for(int i=1;i<=n;++i) R.a[i][1]=(rnd()&1),(sum+=R.a[i][1])%=2;
			R=B*R;
			R=A*R;
			for(int i=1;i<=n;++i){
				if(R.a[i][1]!=sum){
					for(int j=1;j<=n;++j){
						if(i==j) continue;
						int s=0;
						for(int k=1;k<=m;++k){
							(s+=A.a[i][k]*A.a[j][k])%=2;
						}
						if(s==0){
							printf("%d %d\n",min(i,j),max(i,j));
							return;
						}
					}
				}
			}
		}
		puts("-1 -1");
	}
}
namespace sub3{
	inline void main(){
		mt19937 rnd(114514);
		matrix A,B,C,D;
		A.n=n,A.m=m,A.build();
		B.n=m,B.m=n,B.build();
		for(int i=1;i<=n;++i){
			for(int j=1;j<=m;++j){
				B.a[j][i]=A.a[i][j]=read()%k;
			}
		}
		C=A,D=B;
		for(int fick=1;fick<=10;++fick){
			int sum=0;
			matrix R=C;
			for(int i=1;i<=n;++i){
				int x=rnd()%3;
				(sum+=x)%=3;
				for(int j=1;j<=m;++j){
					(R.a[i][j]*=x)%=3;
				}
			}
			R=B*R;
			for(int i=1;i<=n;++i){
				int s=0;
				for(int j=1;j<=m;++j){
					for(int k=1;k<=m;++k){
						s+=A.a[i][j]*R.a[j][k]*D.a[k][i];
					}
				}
				if(s%3!=sum){
					for(int j=1;j<=n;++j){
						if(i==j) continue;
						int tmp=0;
						for(int k=1;k<=m;++k){
							(tmp+=A.a[i][k]*A.a[j][k])%=3;
						}
						if(tmp==0){
							printf("%d %d\n",min(i,j),max(i,j));
							return;
						}
					}
				}
			}
		}
		puts("-1 -1");
	}
}
int main(){
	n=read(),m=read(),k=read();
	if(k==2) sub2::main();
	else sub3::main();
	return 0;
}

记录：