【NOI2013】向量内积

最新推荐文章于 2021-08-08 08:53:47 发布

weixin_33888907

最新推荐文章于 2021-08-08 08:53:47 发布

阅读量169

点赞数

原文链接：http://www.cnblogs.com/shenben/p/6636679.html

版权

#121. 【NOI2013】向量内积

统计

两个

⟨ A, B ⟩ = \sum i = 1 d a i b i = a 1 b 1 + a 2 b 2 + \dots + a d b d

现在有

输入格式

第一行包含

接下来

输出格式

包含两个整数，用空格隔开。

如果存在两个向量

若不存在这样的向量组合，则输出两个

样例一

input

output

2 3

explanation

限制与约定

测试点编号
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

时间限制：

空间限制：

下载

样例数据下载

#include<cstdio>
#define set(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
#define EF if(ch==EOF) return x;
using namespace std;
typedef long long ll;
const int N=1e5+5;
int n,d,K,flag;
ll ans,a[N][101];
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;EF;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void work2(){
    for(int i=1;i<=n;i++){
        for(int j=1;j<i;j++){
            ans=0;
            for(int k=0;k<d;k++){
                ans+=a[i][k]*a[j][k];
            }
            if(!(ans&1)){printf("%d %d\n",j,i);flag=1;return ;}
        }
    }
}
void work3(){
    for(int i=1;i<=n;i++){
        for(int j=1;j<i;j++){
            ans=0;
            for(int k=0;k<d;k++){
                ans+=a[i][k]*a[j][k];
            }
            if(!(ans%K)){printf("%d %d\n",j,i);flag=1;return ;}
        }
    }
}
int main(){
    set(meow);
    n=read();d=read();K=read();
    for(int i=1;i<=n;i++){
        for(int j=0;j<d;j++){
            a[i][j]=read();
        }
    }
    if(K==2) 
        work2();
    else 
        work3();
    if(!flag) puts("-1 -1");
    return 0;
}

75分乱搞

——转自网络：

先考虑k=2的情况
以下运算全部在mod 2的环境下运行。
将n个向量压进一个[Math Processing Error]的矩阵A1里。
设一个[Math Processing Error]的矩阵A2满足[Math Processing Error]

[Math Processing Error]

设矩阵[Math Processing Error]
如果存在一对向量的内积为0
那么等价于矩阵P在非对角线位置上存在一个元素0。

对角线上的0元素干扰了我们的判断。
所以再设一个[Math Processing Error]的矩阵F
满足[Math Processing Error] 其他元素都为0。
由于只有对角线上的元素是非零的，所以可以用个数组存下F。

再设一个[Math Processing Error]的矩阵G，元素全部为1。

再设一个矩阵T=G-F-P

至此，存在一对向量的内积为0已经等价于矩阵T中存在一个元素1。

但同时，计算这些矩阵的时间、空间复杂度都是巨大的。
我们也无法判断矩阵T是否存在元素1，更无法判断他在哪里。
我们一步步来。

随机一个[Math Processing Error]的矩阵X
如果我们能求出矩阵L=[Math Processing Error]，那么如果[Math Processing Error]
那么T矩阵的第i行上就一定存在一个元素1
剩下的事情就很容易了。

但问题是T是一个[Math Processing Error]的矩阵，太巨大了,既存不下，也算不得，怎么办？

我们观察到: [Math Processing Error]
什么意思？分配律！
[Math Processing Error]

G是一个元素全部为1的[Math Processing Error]的矩阵，O(n)算出[Math Processing Error]

F只有对角线上的元素非0，O(n)算出[Math Processing Error]

而求[Math Processing Error]可以利用结合律
求[Math Processing Error]的复杂度是O(dn)
再求[Math Processing Error],复杂度也是O(dn)

完美的避开了高复杂度的运算。

由于矩阵[Math Processing Error]算出来没有1的概率是很小的，所以随机一定的次数大概就没问题了。

对于k=3的情况。由于运算出来的结果除了0,1，还有2
所以上述方法就不能继续沿用了

但我们再进一步观察: [Math Processing Error]
也就是说，如果我们求的是内积平方，那么问题就解决了。
但问题是如果求内积平方就不能用矩阵来做了。

所以我们再观察一下所谓的“内积平方”。
假设d=3吧。
向量A=[a1,a2,a3],向量B=[b1,b2,b3]

[Math Processing Error]

貌似就变成了:
向量A[a1a1,a1a2,a1a3,a2a1,a2a2,a2a3,a3a1,a3a2,a3a3]
向量B[b1b1,b1b2,b1b3,b2b1,b2b2,b2b3,b3b1,b3b2,b3b3]
两者的内积。
也就是说只要把d维向量变成上述[Math Processing Error]维的向量
问题就和k=2的情况几乎完全相同了。
时间复杂度：
对于k=2：O(CND)
对于k=3：O(CND^2)
期望得分：100

#include<ctime>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=1e5+10;
const int Z=110;
int n,m,p,v[N][Z],r[N],b[Z],c[Z][Z];
struct buf{
    char z[1<<25],*s;
    buf():s(z){
        fread(z,1,1<<25,stdin);
    }
    operator int(){
        int x=0,f=1;
        for(;*s<'0'||*s>'9';s++)if(*s=='-')f=-1;
        for(;*s>='0'&&*s<='9';s++)x=x*10+*s-'0';
        return x*f;
    }
}R;
int sol(int i){
    int s=0;
    if(p==2)
        for(int j=1;j<=m;j++) s^=b[j]&v[i][j],b[j]^=v[i][j];
    else
        for(int j=1;j<=m;j++)
            for(int k=1;k<=m;k++)
                s+=c[j][k]*v[i][j]*v[i][k],c[j][k]+=v[i][j]*v[i][k];
    return s%p;
}
void sol(int i,int j){
    int t(0);
    for(int k=1;k<=m;k++) t+=v[i][k]*v[j][k];
    if(!(t%p)){
        if(i>j) swap(i,j);
        printf("%d %d\n",i,j);
        exit(0);
    }
}
int main(){
    n=R;m=R;p=R;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            v[i][j]=R%p;
    for(int i=1;i<=n;i++) r[i]=i;
    srand(20000329);
    random_shuffle(r+1,r+n+1);
    for(int i=1;i<=n;i++)
        if(sol(r[i])!=(i-1)%p)
            for(int j=1;;sol(r[i],r[j++]));
    puts("-1 -1");
    return 0;
}