We will call point a bad if there are different points b and c, not equal to a, from the given set such that angle between vectors and
is acute (i.e. strictly less than
). Otherwise, the point is called good.
The angle between vectors and
in 5-dimensional space is defined as
, where
is the scalar product and
is length of
.
Given the list of points, print the indices of the good points in ascending order.
The first line of input contains a single integer n (1 ≤ n ≤ 103) — the number of points.
The next n lines of input contain five integers ai, bi, ci, di, ei (|ai|, |bi|, |ci|, |di|, |ei| ≤ 103) — the coordinates of the i-th point. All points are distinct.
First, print a single integer k — the number of good points.
Then, print k integers, each on their own line — the indices of the good points in ascending order.
6 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1
1 1
3 0 0 1 2 0 0 0 9 2 0 0 0 5 9 0
0
In the first sample, the first point forms exactly a angle with all other pairs of points, so it is good.
In the second sample, along the cd plane, we can see the points look as follows:
We can see that all angles here are acute, so no points are good.
#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
int a[1010][10];
int abc(int x,int y,int z)
{
//求向量xy和xz夹角正负判断是否为锐角。
int xy=(a[y][0]-a[x][0])*(a[z][0]-a[x][0])+
(a[y][1]-a[x][1])*(a[z][1]-a[x][1])+
(a[y][2]-a[x][2])*(a[z][2]-a[x][2])+
(a[y][3]-a[x][3])*(a[z][3]-a[x][3])+
(a[y][4]-a[x][4])*(a[z][4]-a[x][4]);
if(xy<=0) return 0;
return 1;
}
int sum[1010];
int main()
{
int n;
scanf("%d",&n);
for(int i=0; i<n; i++)
for(int j=0; j<5; j++)
scanf("%d",&a[i][j]);
if(n==1)
printf("%d %d\n",1,1);
else if(n==2)
printf("%d\n%d\n%d\n",2,1, 2);
else if(n>11)
printf("%d\n",0);
else {
int r=0;
for(int i=0; i<n; i++)
{
int flag=0;
for(int j=0; j!=i,j<n; j++)
{
for(int k=j+1; k!=i,k!=j,k<n; k++)
{
if(abc(i,j,k))//<90
{
// printf(" abc=%d ",abc(i,j,k));
// printf("11111111111111\n");
flag=1;
j=n+1;
break;
}
}
}
if(!flag) sum[r++]=i+1;
}
printf("%d\n",r);
if(r)
{
for(int i=0; i<r; i++)
printf("%d ",sum[i]);
}
}
}