洛谷P1224 向量内积
什么毒瘤......
题意:给定n个d维向量,定义向量a和b的内积为
求是否存在两个向量使得它们的内积为k的倍数,并给出任一种方案。k <= 3。
解:很容易想到一个暴力是n2d的。显然我们不能n2枚举,所以要一次性把一个向量跟多个向量判断。
先思考k = 2的情况,显然每个位置和内积非0即1,这启发我们使用二进制。
假如把一个内积看成一个B进制数或者一个多项式,变量是B,我们就能发现,如果两个向量的内积为x,那么这个多项式的值也是x。
这种情况只要B取一个奇数就行了。理由是内积每一项非0即1,而进制为奇数的话,每一项的xi % 2 = 1,奇偶性不变。所以最后加起来和直接加起来的奇偶性相同。
k = 3的时候只要进制为3a + 1就行了。所以最终我们选择7进制。
然后有个很严峻的问题:我们要找一个运算使之与按位乘相对应。先想到了转成指标加法,经过一番推倒之后发现不可行。然后陷入江局......
正解:再观察一波内积式子,您就会发现这个其实是矩阵乘法中的一个位置的计算式......反正我是没发现。
那么令A = n × d的矩阵,B = A * AT,则Bi,j就是i和j的内积。
然后我们只需检验B和全1矩阵(对角线不一定是1)是否相同即可。这有一个经典算法:随机向量法。
随机出来的向量哪一位不同,就表明在全1矩阵的哪一列中存在差异。枚举跟这个向量匹配的向量即可。
k = 3的时候,我们把B中每一个元素都取平方,这样1和2都会变成1。
那么怎么把B中的每个元素取平方呢?
把B中某个元素的式子化开,会有:
然后就做完了......
#include <cstdio>
#include <algorithm>
#include <ctime>
#include <iostream> const int N = ; int a[N][], now[N], C[N], D[N], E[N], MO, F[N];
int n, d; inline bool check(int i, int j) {
int ans = ;
for(int k = ; k <= d; k++) {
(ans += a[i][k] * a[j][k]) %= MO;
}
return ans;
} inline void solve1() {
int T = , f = -;
while((T--) && (f == -)) {
int Sum = ;
for(int i = ; i <= n; i++) {
C[i] = rand() & ;
Sum += C[i];
}
//mul(1, n, d, C, a, D);
for(int i = ; i <= d; i++) {
D[i] = ;
for(int j = ; j <= n; j++) {
D[i] += C[j] * a[j][i];
D[i] &= ;
}
}
//mul(1, d, n, D, aT, E);
for(int i = ; i <= n; i++) {
E[i] = ;
for(int j = ; j <= d; j++) {
E[i] += D[j] * a[i][j];
E[i] &= ;
}
}
//mul_one(1, n, n, C, F);
for(int i = ; i <= n; i++) {
F[i] = ((Sum - C[i]) + C[i] * now[i]) & ;
if(E[i] != F[i]) {
f = i;
break;
}
}
}
if(f == -) {
printf("%d %d\n", f, f);
return;
}
for(int i = ; i <= n; i++) {
if(i == f) {
continue;
}
if(!check(i, f)) {
printf("%d %d\n", std::min(i, f), std::max(i, f));
return;
}
}
return;
} inline void solve2() {
int T = , f = -;
while((T--) && (f == -)) {
int Sum = ;
for(int i = ; i <= n; i++) {
C[i] = rand() % MO;
Sum += C[i];
}
//mul(1, n, d, C, a, D);
for(int i = ; i <= d; i++) {
for(int ii = ; ii <= d; ii++) {
int pos = (i - ) * d + ii;
D[pos] = ;
for(int j = ; j <= n; j++) {
D[pos] += C[j] * a[j][i] * a[j][ii];
D[pos] %= MO;
}
}
}
//mul(1, d, n, D, aT, E);
for(int i = ; i <= n; i++) {
E[i] = ;
for(int j = ; j <= d; j++) {
for(int jj = ; jj <= d; jj++) {
int pos = (j - ) * d + jj;
E[i] += D[pos] * a[i][j] * a[i][jj];
E[i] %= MO;
}
}
}
//mul_one(1, n, n, C, F);
for(int i = ; i <= n; i++) {
F[i] = ((Sum - C[i]) + C[i] * now[i]) % MO;
if(E[i] != F[i]) {
f = i;
break;
}
}
}
if(f == -) {
printf("%d %d\n", f, f);
return;
}
for(int i = ; i <= n; i++) {
if(i == f) {
continue;
}
if(!check(i, f)) {
printf("%d %d\n", std::min(i, f), std::max(i, f));
break;
}
}
return;
} int main() {
srand(time());
int k, x;
scanf("%d%d%d", &n, &d, &k);
MO = k;
bool f = (k == );
for(int i = ; i <= n; i++) {
now[i] = ;
for(int j = ; j <= d; j++) {
scanf("%d", &x);
a[i][j] = x % k;
now[i] += a[i][j] * a[i][j];
}
now[i] %= k;
} f ? solve1() : solve2();
return ;
}
AC代码
题解里还有一种神奇的解法,使用了乘法分配率,每次把一个向量和它上面所有向量的乘积加起来跟(i-1) % MO判断。
分配一下,就是把上面向量的每一维都做前缀和,然后相乘。
这样做其实有一点问题,就是可能有乘积为0的检测不出来。不过上面那种方法也彼此彼此了。
洛谷P1224 向量内积的更多相关文章
- 洛谷 P1224 - [NOI2013] 向量内积(随机化)
洛谷题面传送门 一道很神的随机化. 首先由于我们要求向量点乘 \(\bmod k\) 的值,因此我们可以将所有 \(x_{i,j}\) 都模上 \(k\),显然该操作不影响结果正确性. 注意到这里的 ...
- 洛谷P2520向量
题目传送门 看到数据范围其实就可以确定这是一道结论题. 首先分析,给定你的向量的两个坐标a,b有八种组合方式可以用,但实际上整理一下可以得出实际上只有五种,x/y ±2a,x/y ±2b,x+a,y+ ...
- P1224 [NOI2013]向量内积
传送门 发现这个内积和矩乘有点像,考虑构造一个 $n$ 行 $m$ 列的矩阵 $A$,每一行都是一个题目给定的 $m$ 维向量 设 $B=AA^T$ ,其中 $A^T$ 为 $A$ 的转置矩阵,那么对 ...
- luogu P1224 [NOI2013]向量内积
传送门 挺有意思的一道题 暴力60就是枚举每个向量暴力check,随机选向量就能多骗一些分 然后两个向量内积要模\(k\)为\(0\),那么如果全部不为\(0\)就不合法.先考虑\(k=2\),对于向 ...
- [BZOJ]3243 向量内积(Noi2013)
小C做了之后很有感觉的题目之一,但因为姿势不对调了很久. Description 两个d 维向量A=[a1,a2,...,ad]与B=[b1,b2,...,bd]的内积为其相对应维度的权值的乘积和,即 ...
- 洛谷P3222 [HNOI2012]射箭(计算几何,半平面交,双端队列)
洛谷题目传送门 设抛物线方程为\(y=ax^2+bx(a<0,b>0)\),我们想要求出一组\(a,b\)使得它尽可能满足更多的要求.这个显然可以二分答案. 如何check当前的\(mid ...
- FFT/NTT总结+洛谷P3803 【模板】多项式乘法(FFT)(FFT/NTT)
前言 众所周知,这两个东西都是用来算多项式乘法的. 对于这种常人思维难以理解的东西,就少些理解,多背板子吧! 因此只总结一下思路和代码,什么概念和推式子就靠巨佬们吧 推荐自为风月马前卒巨佬的概念和定理 ...
- 「BZOJ1038」「洛谷P2600」「ZJOI2008」瞭望塔 半平面交+贪心
题目链接 BZOJ/洛谷 题目描述 致力于建设全国示范和谐小村庄的H村村长dadzhi,决定在村中建立一个瞭望塔,以此加强村中的治安. 我们将H村抽象为一维的轮廓.如下图所示: 我们可以用一条山的上方 ...
- 洛谷P4288||bzoj3564 [SHOI2014]信号增幅仪
bzoj3564 洛谷P4288 可以旋转一下坐标轴使得x轴与长轴方向对齐,然后将所有的横坐标变为自身除以放大倍数,然后就做一个最小圆覆盖 #include<cstdio> #includ ...
随机推荐
- 版本控制--git+idea
- Python--文件、文件夹、压缩包、处理模块shutil
高级的 文件.文件夹.压缩包 处理模块 shutil.copyfileobj(fsrc, fdst[, length])将文件内容拷贝到另一个文件中 1 import shutil 2 3 shuti ...
- MyCat数据库中间件 - 分库
MyCat MyCat用于解耦分布式数据库与java,比如分库分表以后,需要查询某条数据时,需要java根据需要查的数据先计算去哪个库查,然而有了Mycat就不用自己计算怎么存储,怎么查询了.MyCa ...
- yield send 的一些使用细节
其实日常中我们使用最多的是 return 很少会使用到 yield 去创造一个生成器.一般就是算算算 算完之后用 return 返回一把. 但是有些情况下 比如需要节约内存不需要一把全部返回,每次使用 ...
- python数据结构与算法第十二天【快速排序】
1. 原理如图所示: 2.代码实现 def quick_sort(alist, start, end): """快速排序""" # 递归的退 ...
- Delphi MDI 子窗体的创建和销毁 [zhuan]
1.如果要创建一个mdi child,先看是否有这个child 存在,如果有,则用它,如果没有再创建 //该函数判断MDI 子窗体是否存在,再进行创建和显示function isInclude(for ...
- vue事件綁定
事件綁定可以是一個句子,一個函數名稱,也可以是一個函數. 事件修飾符,按鍵修飾符.
- hdu-2087(kmp)
题意:模板题,在第一个串中有几个第二个串 解题思路:板子题,拿来练手的: 代码: #include<iostream> #include<algorithm> #include ...
- Nginx websocket反向代理
L:106 现在主流的反向代理,通过长链接可以从服务器推送数据到页面 升级成websocket反向代理必须根据上面的配置做配置 缺点无法多路复用,也就是没办法并行 我们测试下Websocket反向代理 ...
- Nginx 反向代理如何连接上游服务器
L:92 想上游服务器先建立TCP连接 如三次握手 下面指令可以控制握手时间 proxy_next_upstream 指令当出现502可以换个上游服务器 Tcp keepalive 一般都是由进程在 ...