【油滴扩展】dfs/STL next_permutation

洛谷P1378 油滴扩展
题目描述
在一个长方形框子里，最多有N（0≤N≤6）个相异的点，在其中任何一个点上放一个很小的油滴，那么这个油滴会一直扩展，直到接触到其他油滴或者框子的边界。必须等一个油滴扩展完毕才能放置下一个油滴。那么应该按照怎样的顺序在这N个点上放置油滴，才能使放置完毕后所有油滴占据的总体积最大呢？（不同的油滴不会相互融合）

注：圆的面积公式V=pirr，其中r为圆的半径。

输入格式
第1行一个整数N。

第2行为长方形边框一个顶点及其对角顶点的坐标，x,y,x’,y’。

接下去N行，每行两个整数xi,yi，表示盒子的N个点的坐标。

以上所有的数据都在[-1000，1000]内。

输出格式
一行，一个整数，长方形盒子剩余的最小空间（结果四舍五入输出）

输入输出样例
输入 #1复制
2
20 0 10 10
13 3
17 7
输出 #1复制
50

这道题蒟蒻博主写了两天。。。圆周率一开始写3.14直接wa声一片，看了题解才知道最少七位，最好调用头文件<math.h>中的反余弦double pi = acos(-1),非常好记，妈妈再也不用担心我卡圆周率了。
进入正题，当时是搜dfs算法找到的这道题，但自己dfs写全排列太麻烦了，不如用next_permutation多爽，结果整整被这个卡了两天。。。
卡了两处，第一个是我用两个数组分别存坐标的x和y，然后dowhile里面一起next_permutation，debug了很久才发现跑着跑着不对应了。。。第二个是我醒悟后用结构体存坐标，然后只有样例4wa了，很不巧样例4数据还很大。。。我插桩看了半天才发现next_permutation之前如果不给结构体排序的话，会少前面几种排列。。。天啊，我太天真了。

附上AC代码。

#include <iostream>
#include<bits/stdc++.h>
using namespace std;
#define pi 3.14159265358979323846
int n;
double ans,sum;
int x,y,xx,yy;
int xi[10],yi[10];
double r[10];//存每个油滴的半径
bool vis[10];
struct node
{
    int xi,yi;
}s[10];
bool cmp(node a,node b)
{
    return a.xi!=b.xi?a.xi<b.xi:a.yi<b.yi;
}
int main()
{
    cin>>n;
    cin>>x>>y>>xx>>yy;
    sum=abs(x-xx)*abs(y-yy);
    ans=sum;
    for (int i = 0; i<n; i++)
    {
        cin>>s[i].xi>>s[i].yi;
    }
    sort(s,s+n,cmp);
    do//先根据边界初始化油滴半径，vis标记用过的油滴，之后如果碰到，就取min
    {
        double re=sum;
        memset(vis,0,sizeof(vis));
        for(int i = 0; i<n; i++)
        {
            r[i] = min(min(abs(x-s[i].xi),abs(xx-s[i].xi)),min(abs(y-s[i].yi),abs(yy-s[i].yi)));//根据边界初始化油滴半径
            vis[i]=1;
            for(int j =0; j<n; j++)
            {
                if(j==i)continue;//跳过自己
                if(vis[j])//j油滴用过了，需要i重新判断
                {
                    double d = sqrt(abs(s[i].xi-s[j%n].xi)*abs(s[i].xi-s[j%n].xi)+abs(s[i].yi-s[j%n].yi)*abs(s[i].yi-s[j%n].yi));//两油滴距离
                    if(d<=r[j])//两油滴距离小于等于已有油滴的半径
                    {
                        r[i]=0;
                        break;
                    }else
                    {
                        r[i] = min(r[i],d-r[j]);
                    }

                }
            }
            re-=pi*r[i]*r[i];
        }
        ans = min(ans,re);

    }
    while(next_permutation(s,s+n,cmp));
    cout<<fixed<<setprecision(0)<<ans<<endl;
    return 0;
}
/*
样例3
-98 5 30 30
-42 11
-51 17
-11 22

2547
*/
/*
样例5 
-1 67 972 591 
626 266 
322 316 
32 334 
361 183 
480 345 

300682
*/