BJ 集训测试13 平行

http://www.elijahqi.win/archives/2845
二维直角坐标系中有N个点，它们的横、纵坐标的绝对值都不超过10。

一次操作，你可以选择某3个不共线的点A、B、C。然后以AB为对角线、C为另一顶点绘制平行四边形ABCD，然后用点D替换点C（即删去点C保留点A、B、D）。则D是确定并且唯一的。

最初，N个点的坐标互不相同。操作过程中这些点的坐标可以重合。但是，要求所有新点的坐标的绝对值不超过109。
请你判断：能否在2500次操作以内，使得所有N个点都出现在第一象限？即，通过若干次操作使得所有点的横、纵坐标非负。如果可以，请给出一个方案。

考虑横纵坐标绝对值范围只有10，那么如果有三个不共线得点横纵坐标均大于5，那么一定可以把所有点都移到一象限… 考虑叉积若向量b*向量a<0则说明b在a的逆时针方向 以此来做题可以简化代码难度
那么找到一个三角形，然后可以发现他沿着两条边操作两次可以变成平移，然后再考虑他的移动的向量，然后找到两种平移方式就一定可以构成两个单位向量，然后操作即可…

#include<cstdio>
#include<cctype>
#include<algorithm>
#define N 440 
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if(T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
struct node{
    int x,y;
    inline friend node operator +(const node &a,const node &b){
        return (node){a.x+b.x,a.y+b.y};
    }
    inline friend node operator -(const node &a,const node &b){
        return (node){a.x-b.x,a.y-b.y};
    }
}p[N];int n;
struct node1{
    int a,b,c;
};
vector<node1> ans;
inline int cross(const node &a,const node &b){
    return a.x*b.y-a.y*b.x;
}
inline void trans(int &p1,int &p2,int &p3){
    ans.push_back((node1){p2,p3,p1});
    ans.push_back((node1){p1,p2,p3});
    int x1=p1,x2=p2,x3=p3;p3=x1;p2=x3;p1=x2;
}
inline void solve(int p1,int p2,int p3){
    int s1,s2;node a=p[p2]-p[p1],b=p[p3]-p[p1];
    s1=b.x,s2=-a.x;if (s1*a.y+s2*b.y<0) s1=-s1,s2=-s2;
    if (s1>0) for (int i=1;i<=s1*20;++i) trans(p1,p2,p3);
    else for (int i=1;i<=-s1*20;++i) trans(p2,p1,p3);
    if (s2>0) for (int i=1;i<=s2*20;++i) trans(p1,p3,p2);
    else for (int i=1;i<=-s2*20;++i) trans(p3,p1,p2);
    s1=b.y;s2=-a.y;if (s1*a.x+s2*b.x<0) s1=-s1,s2=-s2;
    if (s1>0) for (int i=1;i<=s1*20;++i) trans(p1,p2,p3);
    else for (int i=1;i<=-s1*20;++i) trans(p2,p1,p3);
    if (s2>0) for (int i=1;i<=s2*20;++i) trans(p1,p3,p2);
    else for (int i=1;i<=-s2*20;++i) trans(p3,p1,p2);
    for (int i=0;i<ans.size();++i) p[ans[i].c]=p[ans[i].a]+p[ans[i].b]-p[ans[i].c];
}
int main(){
    freopen("t2.in","r",stdin);
    n=read();int pos=-1;
    for (int i=1;i<=n;++i) p[i].x=read(),p[i].y=read();
    for (int i=1;i<=n;++i)
        if (cross(p[1]-p[i],p[2]-p[i])) {pos=i;break;}
    if(pos==-1) {puts("-1");return 0;}solve(1,2,pos);
    for (int i=3;i<=n;++i){
        if(i!=pos){
            if (cross(p[1]-p[i],p[pos]-p[i])) ans.push_back((node1){1,pos,i});else
            if (cross(p[2]-p[i],p[pos]-p[i])) ans.push_back((node1){2,pos,i});else
            ans.push_back((node1){1,2,i});
        }
    }printf("%d\n",ans.size());
    for (int i=0;i<ans.size();++i) printf("%d %d %d\n",ans[i].a,ans[i].b,ans[i].c);
    return 0;
}