LA4356 Fire-Control System

题目

LA4356 Fire-Control System

题解

我开始的想法：半径最多n个，夹角（与x轴正半轴夹角）最多n个，离散一下然后按半径排序，对于每种半径在那一圈上用前缀和搞一搞然后用滑动窗口确定最小面积，然后发现这样不仅不好写，而且精度误差对离散有很大影响，就算保存半径的平方，夹角也没办法避免误差，所以GG。
做法：按夹角排序，枚举每个点i，设点i到原点的距离为dis[i],找到满足$dis[j] <= dis[i]$的所有j，成一个序列，因为是环不好操作所以断环，变为2倍长度，然后枚举是哪个区间的被计算，因为这里一个下标对应一个点，所以枚举其实位置可以方便的知道结束位置，然后计算扇形夹角的时候注意+=2*pi。

注：之前因为手残一直狂RE不止。

代码

//QWsin
#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=5000+10;
const double pi=acos(-1.0);

int n,k,x[maxn],y[maxn];
struct Node{
    double jiao;int dis;
    Node(double jiao=0,int dis=0):jiao(jiao),dis(dis){}
    bool operator < (const Node &rhs)const{
        return jiao < rhs.jiao;
    }
}a[maxn];

set<int>vis;

int t[maxn*2],kase=0;
inline void print(double ans){printf("Case #%d: %.2f\n",++kase,ans);}
inline void solve()
{
    for(int i=1;i<=n;++i) 
    {
        scanf("%d%d",x+i,y+i);
        a[i].jiao=atan2(y[i],x[i]);
        if(a[i].jiao<0) a[i].jiao+=2*pi;
        a[i].dis=x[i]*x[i]+y[i]*y[i];
    }
    if(!k) {print(0);return ;}
    vis.clear();
    sort(a+1,a+n+1);

    double ans=1e10;
    for(int i=1;i<=n;++i) if(!vis.count(a[i].dis))
    {
        int p=0;vis.insert(a[i].dis);
        for(int j=1;j<=n;++j) if(a[j].dis <= a[i].dis) t[++p]=j;
        if(p<k) continue;
        for(int j=1;j<p;++j) t[p+j]=t[j];

        for(int j=1;j<=p;++j) 
        {
            double dj=a[t[j+k-1]].jiao-a[t[j]].jiao;//之前a数组的下标写的是j+k-1,23333
            if(j+k-1>p) dj+=2*pi;
            ans=min(ans,dj/2*a[i].dis);
        }
    }

    print(ans);
}

int main()
{
    while(scanf("%d%d",&n,&k)==2){
        if(!n&&!k) break; solve();
    }
    return 0;
}