• [DLX反复覆盖] hdu 3656 Fire station


    题意:

    N个点。再点上建M个消防站。

    问消防站到每一个点的最大距离的最小是多少。

    思路:

    DLX直接二分推断TLE了。

    这时候一个非常巧妙的思路

    我们求的距离一定是两个点之间的距离

    因此我们把距离都求一遍排序一下。

    然后用下标二分  这样就AC了。

    代码:

    #include"stdio.h"
    #include"algorithm"
    #include"string.h"
    #include"iostream"
    #include"cmath"
    #include"queue"
    #include"map"
    #include"vector"
    #include"string"
    using namespace std;
    #define eps 1e-6
    #define N 55*55
    #define RN 55
    #define CN 55
    double dist[N];
    int kx;
    struct DLX
    {
        int n,m,C;
        int U[N],D[N],L[N],R[N],Row[N],Col[N];
        int H[RN],S[CN],cnt,ans[RN];
        void init(int _n,int _m)
        {
            n=_n;
            m=_m;
            for(int i=0; i<=m; i++)
            {
                S[i]=0;
                U[i]=D[i]=i;
                L[i]=(i==0?m:i-1);
                R[i]=(i==m?

    0:i+1); } C=m; for(int i=1; i<=n; i++) H[i]=-1; } void link(int x,int y) { C++; Row[C]=x; Col[C]=y; S[y]++; U[C]=U[y]; D[C]=y; D[U[y]]=C; U[y]=C; if(H[x]==-1) H[x]=L[C]=R[C]=C; else { L[C]=L[H[x]]; R[C]=H[x]; R[L[H[x]]]=C; L[H[x]]=C; } } void del(int x) { for(int i=D[x]; i!=x; i=D[i]) { R[L[i]]=R[i]; L[R[i]]=L[i]; } } void rec(int x) { for(int i=U[x]; i!=x; i=U[i]) { R[L[i]]=i; L[R[i]]=i; } } int used[CN]; int h() { int sum=0; for(int i=R[0]; i!=0; i=R[i]) used[i]=0; for(int i=R[0]; i!=0; i=R[i]) { if(used[i]==0) { sum++; used[i]=1; for(int j=D[i]; j!=i; j=D[j]) for(int k=R[j]; k!=j; k=R[k]) used[Col[k]]=1; } } return sum; } int dance(int x) { if(x+h()>=cnt || x+h()>kx) return 0; if(R[0]==0) { cnt=min(cnt,x); if(cnt<=kx) return 1; return 0; } int now=R[0]; for(int i=R[0]; i!=0; i=R[i]) { if(S[i]<S[now]) now=i; } for(int i=D[now]; i!=now; i=D[i]) { del(i); for(int j=R[i]; j!=i; j=R[j]) del(j); if (dance(x+1)) return 1; for(int j=L[i]; j!=i; j=L[j]) rec(j); rec(i); } return 0; } } dlx; struct node { double x,y; } dian[55]; double dis(double x1,double y1,double x2,double y2) { return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)); } int main() { int t; cin>>t; while(t--) { int n; scanf("%d%d",&n,&kx); for(int i=1;i<=n;i++) scanf("%lf%lf",&dian[i].x,&dian[i].y); int cnt=0; for(int i=1;i<=n;i++) { for(int j=i;j<=n;j++) { dist[cnt++]=dis(dian[i].x,dian[i].y,dian[j].x,dian[j].y); } } sort(dist,dist+cnt); int l=0,r=cnt-1; double ans; while(l<=r) { int mid=(l+r)/2; dlx.init(n,n); for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { double tep=dis(dian[i].x,dian[i].y,dian[j].x,dian[j].y); if(tep>dist[mid]) continue; dlx.link(i,j); } } dlx.cnt=999; if(dlx.dance(0)) { ans=dist[mid]; r=mid-1; } else l=mid+1; } printf("%.6f ",ans); } return 0; }



  • 相关阅读:
    我的Win32开发抉择,Delphi老将复出
    【月入41万】Mono For Android中使用百度地图SDK
    “开源”是什么,能吃吗?
    (转)oracle中用户删除不了,ORA-01940提示 “无法删除当前已连接用户”
    Oracle 去除两边空格
    (转)C#中键值对类型Hashtable与Dictionary比较和相关用法
    (转)对Oracle导出文件错误和DMP文件结构的分析,EXP-00008: 遇到 ORACLE 错误 904 ORA-00904: "MAXSIZE": invalid identifier
    (转)ajax.dll,ajaxpro.dll的区别和用法
    npoi实现 从固定的行读取数据作为表头并返回datable
    字符串格式化测试
  • 原文地址:https://www.cnblogs.com/yxwkf/p/5403079.html
Copyright © 2020-2023  润新知