2021-ICPD昆明站-I Mr. Main and Windmills

传送门

题意：给定一个起点坐标s一个终点坐标t，然后给定n(1~1000)个点的坐标，这些点都在s-t的一侧（而且不会在s-t上）。

一个点旋转一次，比如最开始a在b的视角左边，后来变成了不是右边（这里默认包含在一条视线上）。

官方图：

然后有m次查询，每次查询编号为h的点旋转k次的时候点ans的坐标（ans从s到t滑行）。

题解：

1）以原始点为基础，记录该点“贡献”（比如pi-pj与s-t相交于点a，那么g[i].push_back(a),g[j].push_back(a)）过的所有与s-t的交点。每遇到一个交点表示“旋转”一次。

2）求交点：其中一条线的两点为pi,pj，另一条线上的两点为s,t。那么交点为s+(t-s)*x。若a为交点，则x表示sa.len()/(sa.len()+st.len())。（简单推一下就好了）

3）x怎么求？？叉积表示平行四边形的面积，即底*高，如果底相等，那么面积之比则为两个平行四边形的底到另一条平行线之比，然后利用相似。

emmmmmmm不太直观，总之先知道叉积表示的是平行四边形的面积，然后还有利用底*高，以及两个高的比（相似）得到sa.len()与st.len()的比例。

建议直接看图。（假装好理解一些）

总结：

1）其实也不难想，但是感觉如果不是做过一些计算几何的题目，发到这种题还是会发怵的（ps:精度问题搞人心态，而且计算几何没咋写过的话代码也不是那么容易搞）。——比如说我。

2）自己开始向的时候，也想到了记录所有可能在s,t上的交点，然后一个一个记录是哪些点共线的。以s-t上的交点为基础记录贡献点的下标，最差的空间复杂度1e6*1e3（牛客竟然能提交，不过一直在61.54%的正确率）。

3）无奈看别人代码，着实震惊了，太，简洁了（我今后遇到计算几何题也要尽量自己推，不抄模板）。

4）xcpc最不能少的就是思考！！！如果将时间浪费在写代码上就不好了——这一题，写了这么多还是没有结果。如果能再多想一下，不接受1e6*1e3的空间复杂度，那还是比较容易想出优化方案的：换点的下标为基础就好了，然后记录产生的所有贡献点。

 

代码：

#include<bits/stdc++.h>

#define pb push_back
#define dbg(x) cout<<#x<<"===="<<x<<endl;
#define mst(a,x) memset(a,x,sizeof(a))
using namespace std;
void read(int &x){
    int res=0,f=1;
    char c=getchar();
    while(c<'0'||c>'9'){
        if(c=='-') f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9') res=(res<<3)+(res<<1)+(c-'0'),c=getchar();
    x=res*f;
}
void OUT(int x){
    if(x<0) putchar('-'),x=-x;
    if(x>9) OUT(x/10);
    putchar(x%10+'0');
}
void print(int a,char c){OUT(a),putchar(c);}

const int maxn=1e3+10;
const double eps=1e-8;
const double pi=acos(-1.0);
int sgn(double x){
    if(fabs(x)<eps) return 0;
    if(x<0) return -1;
    else return 1;
}

int n,m,h,k;
struct Point{
    double x,y;
    Point(){}
    Point(double _x,double _y){x=_x,y=_y;}
    void input(){scanf("%lf%lf",&x,&y);}
    void output(){printf("%.10lf %.10lf\n",x,y);}
    Point operator-(const Point &b)const{return Point(x-b.x,y-b.y);}
    double operator^(const Point &b)const{return x*b.y-y*b.x;}
    Point operator*(const double &k)const{return Point(x*k,y*k);}
    Point operator+(const Point &b)const{return Point(x+b.x,y+b.y);}
}p[maxn],s,t;
vector<double> g[maxn];

signed main(){
    read(n),read(m);
    s.input(),t.input();
    // if(sgn(s.x-t.x)>0) swap(s,t);//毒瘤代码，不知道发生了这个，有它就wa。反正我是不想纠结了
    for(int i=1;i<=n;i++) p[i].input();
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){//这样写起来更容易看懂很多，从i+1开始也是可以的
            if(i==j) continue;
            double a1=(p[i]-p[j])^(s-p[j]);//面积s1
            double a2=(p[i]-p[j])^(t-p[j]);//面积s2
            if(a1*a2<=eps){//s,t不能再pi,pj的同一边（在同一边就与线段s-t没有交点了）
                a2=-a2;//使a1,a2同号
                g[i].pb(a1/(a1+a2));//比值x=a1/(a1+a2)
            }
        }
    }
    for(int i=1;i<=n;i++) sort(g[i].begin(),g[i].end());//排序，按顺序贡献交点
    while(m--){
        read(h),read(k);
        if(g[h].size()<k) puts("-1");
        else{
            Point ans=s+(t-s)*g[h][k-1];
            ans.output();//偷个懒
        }
    }
    return 0;
}