题目链接

BZOJ3243

题解

模数只有\(2\)或\(3\),可以大力讨论

如果模数为\(2\),乘积结果只有\(1\)或\(0\)

如果一个向量和前面所有向量乘积都为\(1\),那么其和前面向量前缀和的乘积就唯一确定

我们维护向量前缀和,第一个乘积情况不符的向量一定是答案,然后再枚举另一个向量即

\(O(nd)\)

如果模数为\(3\),乘积如果不为\(0\),还可以为\(1\)或\(2\),我们讨论的方法就不适用了

其实还是可以的

\[1^2 = 2^2 = 1 \pmod 3
\]

我们只要维护平方和即可

如何维护平方和?

\[(\sum\limits_{i = 1}^{d} a_ib_i)^2 = \sum\limits_{i = 1}^{d} \sum\limits_{j = 1}^{d} a_ia_jb_ib_j
\]

就相当于原来的\(d\)维向量变成了\(d^2\)维,\(O(nd^2)\)也是可以过的

#include<iostream>
#include<cstdio>
#include<cmath>
#include<ctime>
#include<cstring>
#include<algorithm>
#define LL long long int
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<' '; puts("");
using namespace std;
const int maxn = 100005,maxm = 105,INF = 1000000000;
inline int read(){
int out = 0,flag = 1; char c = getchar();
while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
return out * flag;
}
int n,d,P;
int a[maxn][maxm],id[maxn];
int sum[maxm],Sum[maxm][maxm];
int mult(int* a,int* b){
int re = 0;
for (int i = 1; i <= d; i++) re = (re + a[i] * b[i] % P) % P;
return re;
}
int Mult(int a[],int b[][maxm]){
int re = 0;
for (int i = 1; i <= d; i++)
for (int j = 1; j <= d; j++)
re = (re + a[i] * a[j] * b[i][j] % P) % P;
return re;
}
void solve1(){
for (int i = 1; i <= d; i++) sum[i] = a[id[1]][i];
for (int i = 2; i <= n; i++){
int t = mult(a[id[i]],sum);
if (t != ((i - 1) & 1)){
for (int k = 1; k < i; k++)
if (!mult(a[id[k]],a[id[i]])){
printf("%d %d\n",min(id[k],id[i]),max(id[i],id[k]));
break;
}
return;
}
for (int j = 1; j <= d; j++) sum[j] = (sum[j] + a[id[i]][j]) % P;
}
printf("-1 -1\n");
}
void solve2(){
for (int i = 1; i <= d; i++)
for (int j = 1; j <= d; j++)
Sum[i][j] = a[id[1]][i] * a[id[1]][j] % P;
for (int i = 2; i <= n; i++){
int t = Mult(a[id[i]],Sum);
if (t != (i - 1) % P){
for (int k = 1; k < i; k++)
if (!mult(a[id[k]],a[id[i]])){
printf("%d %d\n",min(id[k],id[i]),max(id[i],id[k]));
break;
}
return;
}
for (int j = 1; j <= d; j++)
for (int k = 1; k <= d; k++)
Sum[j][k] = (Sum[j][k] + a[id[i]][j] * a[id[i]][k] % P) % P;
}
printf("-1 -1\n");
}
int main(){
srand(time(NULL));
n = read(); d = read(); P = read();
for (int i = 1; i <= n; i++)
for (int j = 1; j <= d; j++)
a[i][j] = read() % P;
for (int i = 1; i <= n; i++) id[i] = i;
random_shuffle(id + 1,id + 1 + n);
if (P == 2) solve1();
else solve2();
return 0;
}

BZOJ3243 [Noi2013]向量内积 【乱搞】的更多相关文章

  1. BZOJ3243 NOI2013向量内积(随机化)

    考虑奇技淫巧. 首先是k=2.对向量维护一个前缀和,每次将当前向量与前缀和点乘.如果点乘结果不等于i-1&1,说明当前向量至少和之前的某个向量的数量积是2的倍数,暴力找就可以了.当然等于i-1 ...

  2. 【BZOJ-3243】向量内积 随机化 + 矩阵

    3243: [Noi2013]向量内积 Time Limit: 10 Sec  Memory Limit: 256 MBSec  Special JudgeSubmit: 1249  Solved:  ...

  3. 【fake题解】[NOI2013]向量内积

    [fake题解][NOI2013]向量内积 做法1 大暴力.哪里不会T哪里. 做法2 所有数都%=k不影响结果.(废话 k的取值只有2和3,所以肯定是要分类讨论的.k=2肯定简单些啦. k=2 出现的 ...

  4. [Noi2013]向量内积

    来自FallDream的博客,未经允许,请勿转载,谢谢. 两个d 维向量A=[a1,a2,...,ad]与B=[b1,b2,...,bd]的内积为其相对应维度的权值的乘积和,即: $\sum_{i=1 ...

  5. P1224 [NOI2013]向量内积

    传送门 发现这个内积和矩乘有点像,考虑构造一个 $n$ 行 $m$ 列的矩阵 $A$,每一行都是一个题目给定的 $m$ 维向量 设 $B=AA^T$ ,其中 $A^T$ 为 $A$ 的转置矩阵,那么对 ...

  6. luogu P1224 [NOI2013]向量内积

    传送门 挺有意思的一道题 暴力60就是枚举每个向量暴力check,随机选向量就能多骗一些分 然后两个向量内积要模\(k\)为\(0\),那么如果全部不为\(0\)就不合法.先考虑\(k=2\),对于向 ...

  7. BZOJ3243/UOJ121 [Noi2013]向量内积

    本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ ...

  8. 3243: [Noi2013]向量内积 - BZOJ

    Description 两个d 维向量A=[a1,a2,...,ad]与B=[b1,b2,...,bd]的内积为其相对应维度的权值的乘积和,即: 现有 n 个d 维向量x1,...,xn ,小喵喵想知 ...

  9. 【uoj121】 NOI2013—向量内积

    http://uoj.ac/problem/121 (题目链接) 题意 给出${n}$个${d}$维向量,问是否有两个不同的向量的内积是${k}$的倍数. Solution 又卡了一上午常数,我弃了T ...

随机推荐

  1. Hadoop(2)-CentOS下的jdk和hadoop的安装与配置

    准备工作 下载jdk8和hadoop2.7.2 使用sftp的方式传到hadoop100上的/opt/software目录中 配置环境 如果安装虚拟机时选择了open java,请先卸载 rpm -q ...

  2. python_day4_shopping

    购物车例子,实现显示商品信息,输入商品编号并且可以减去自己的存入余额,当商品价格大于自己的余额的时候,直接退出:当不再选择商品的时候,退出显示余额和已经添加的商品. #购物车程序 product_li ...

  3. python中string,time,datetime三者之间的转化

    这里time特指import time中的对象,datetime 特指from datetime import datetime中的对象,string指python自带的字符数据类型. 从使用的情况来 ...

  4. ESP32 学习笔记 - 环境搭建

    打开终端 输入命令 sudo apt-get install gcc git wget make libncurses-dev flex bison gperf python python-seria ...

  5. centos 安装java1.8

    https://www.cnblogs.com/xuliangxing/p/7066913.html

  6. B -- POJ 1208 The Blocks Problem

    参考:https://blog.csdn.net/yxz8102/article/details/53098575 https://www.cnblogs.com/tanjuntao/p/867892 ...

  7. HDU 5530:Pipes Selection

    题意: 给定长度为\(L\),元素总和为\(S\)的非负整数序列\(A\),对于每一个\(1 \leq i \leq S\),求出:所有满足\(\sum_{j=l}^rA_j=i\)的二元组\((l, ...

  8. js倒计时页面跳转

    HTML: <p><span id="timer">60</span>s 后跳转到百度首页</p> JS: //倒计时方法 func ...

  9. Java基础——内部类

    一.什么是内部类 将一个类定义在另一个类里面或者一个方法里面,这样的类称为内部类 内部类所在的类在编译成功后,会出现这样两个class文件:OuterClass.class和OuterClass$In ...

  10. 分支push不上去的问题

    还原一下现场,我在自己的项目里面,从master里面checkout的一个分支,当我在我这个分支里面进行 push代码的操作,我突然发现我的代码不能执行push的操作,如图 这个原因是由于远端的仓库没 ...