链接:https://ac.nowcoder.com/acm/contest/329/A
来源:牛客网
众所周知,处女座是数学大师。他定义了k维空间里的处女座点。
对于给出的k维度空间上N个点,处女座点满足:
对于这个点P和空间里任意其他两个点P1、P2,有dot(→PP1,→PP2)<0;
现在给你一个k维空间和这N个点,请求出这里面所有的处女座点。
Hint: dot(→A,→B)dot(A→,B→)表示两个向量的点乘(内积)。两个向量→A={a1,a2,⋯,an}和→B={b1,b2,⋯,bn}的点乘(内积)定义为:dot(→A,→B)=a1∗b1+a2∗b2+⋯+an∗bn
输入描述:
输入数据第一行包含一个整数T,表示用例组数。 对于每组用例,第一行两个整数 N,k。 接下来 N 行,每行k个整数,表示每个点的坐标。
输出描述:
对于每一组数据,先输出一个整数ans,表示处女座点的个数。 接下来ans行,每行k个数,表示每个处女座点的坐标。 如果有多个点满足条件,则按照输入数据中出现的顺序进行输出。
示例1
输入
1 3 2 0 0 1 1 -1 -1
输出
1 0 0
备注:
3≤N≤ 1≤k≤10 T≤100
对于每一个点,保证每一维坐标的绝对值在以内
1.处女座点的数量最多为 1
证明:
2. 若K维空间一个点集{P1,P2,⋯,PnP1,P2,⋯,Pn}存在处女座点,则𝑵≤𝑲+𝟐≤𝟐𝑲+𝟏
严格的证明需要利用高等代数。但是利用数学归纳法或线性代数方法证明同条件下证明𝑁≤2𝐾+1是显然的。 (提示:将处女座点的定义弱化为向量夹角为钝角或直角后即可证明)
综上所述,若𝑁>𝐾+2则直接输出0,否则暴力检查即可。由于𝐾≤10,用𝑁>2𝐾+1判断也是可以的。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e4+5;
int main()
{
ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
int T,n,kk;
cin>>T;
while(T--)
{
ll a[maxn][15];
cin>>n>>kk;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=kk;j++)
cin>>a[i][j];
}
if(n>kk+2)
{
cout<<0<<endl;
continue;
}
else
{
int x;
for(int i=1;i<=n;i++){
int f=1;
for(int j=1;j<=n;j++) {
for(int k=j+1;k<=n;k++) {
ll ans=0;
if(j==i||k==i)continue;
for(int l=1;l<=kk;l++)
ans+=(a[j][l]-a[i][l])*(a[k][l]-a[i][l]);
if(ans>=0){
f=0;
break;
}
}
}
if(f)
{
cout<<1<<endl;
for(int lll=1;lll<=kk;lll++)
cout<<a[i][lll]<<" ";
cout<<endl;
break;
}
}
}
}
return 0;
}