- Five Dimensional Points(暴力)

ou are given set of n points in 5-dimensional space. The points are labeled from 1 to n. No two points coincide.

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.


Input

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.

Output

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.

Example
Input
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
Output
1
1
Input
3
0 0 1 2 0
0 0 9 2 0
0 0 5 9 0
Output
0
Note

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.

一个点a如果存在b、c点,向量ab与向量ac的夹角小于90度,则说明a是bad;,否则GOOD;


根据公式只需要判断ab与向量ac的乘积是否大于0就可以了。
思路:从2,3维空间超过5,7个点时不存在“good”的点,可以简单推知五维空间内,超过11个点时不存在“good”的点,那么点数小于11时暴力,大于11时输出0.
#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]);
    }
 }
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值