HDU4116 Fruit Ninja (2011 Asia ChengDu Regional Contest)

被这道题玩死了要。。。

对每个圆枚举得到和其他圆切线斜率，两个切线之间的斜率范围是可以经过的。得到的许多斜率范围可以转化为区间覆盖，nlogn排序得解，总复杂度O(n^2logn)。

要特别处理圆相交、包含的情况。

切线的斜率范围，是按 -pi~pi 算的，算晕了很久。。

然后就是恐怖的TLE，最后发现交C++TLE，交G++就过了，额滴个神，累死了，调一晚上。。

还有，这题用sort比用qsort快了5s。。。

  1 #include<stdio.h>
  2 #include<string.h>
  3 #include<stdlib.h>
  4 #include<math.h>
  5 #include<algorithm>
  6 #include<iostream>
  7 using namespace std;
  8 const int maxn = 1 << 10 | 1;
  9 const double pi = acos(-1.0);
 10 const double eps = 1e-6;
 11 typedef struct
 12 {
 13     double x, y, r;
 14 } Cir;
 15 typedef struct
 16 {
 17     double site;
 18     bool se;
 19 } Seg;
 20 Cir ra[maxn];
 21 Seg s[maxn << 3];
 22 int n, stp, ans;
 23 inline double max(double a, double b)
 24 {
 25     return a > b ? a : b;
 26 }
 27 inline double min(double a, double b)
 28 {
 29     return a < b ? a : b;
 30 }
 31 inline double CalDis(double a, double b)
 32 {
 33     return sqrt(a * a + b * b);
 34 }
 35 inline bool IinJ(int i, int j, double ijdis)
 36 {
 37     return ra[i].r + ijdis <= ra[j].r;
 38 }
 39 inline bool IcutJ(int i, int j, double ijdis)
 40 {
 41     return ijdis <= ra[i].r + ra[j].r;
 42 }
 43 void AddSeg(double start, double end)
 44 {
 45     if(start + pi < -eps) start += 2 * pi;
 46     if(start - pi > -eps) start -= 2 * pi;
 47     if(end + pi < -eps) end += 2 * pi;
 48     if(end - pi > -eps) end -= 2 * pi;
 49     if(start - end > eps) AddSeg(start, pi - eps * 2), AddSeg(-pi, end);
 50     else
 51     {
 52         s[stp].site = start, s[stp].se = 1;
 53         ++ stp;
 54         s[stp].site = end, s[stp].se = 0;
 55         ++ stp;
 56     }
 57 }
 58 int comp(const void *a, const void *b)
 59 {
 60     if(fabs((*(Seg *)a).site - (*(Seg *)b).site) < eps)
 61         return -((*(Seg *)a).se && !(*(Seg *)b).se);
 62     return (*(Seg *)a).site > (*(Seg *)b).site ? 1 : -1;
 63 }
 64 /*
 65 //sort
 66 bool comp(const Seg &a, const Seg &b)
 67 {
 68       if(fabs(a.site - b.site) < eps)
 69             return a.se && !b.se;
 70       return a.site < b.site;
 71 }
 72 */
 73 void Cal(int i)
 74 {
 75     int sum, j;
 76     for(j = sum = stp = 0; j < n; ++ j)
 77     {
 78         if(i != j)
 79         {
 80             double ijdis = CalDis(ra[i].x - ra[j].x, ra[i].y - ra[j].y);
 81             double xlij = atan2(ra[i].y - ra[j].y, ra[i].x - ra[j].x);
 82             double xlji = atan2(ra[j].y - ra[i].y, ra[j].x - ra[i].x);
 83             double asipj = asin((ra[i].r + ra[j].r) / ijdis);
 84             double asimj = asin((ra[i].r - ra[j].r) / ijdis);
 85             if(IinJ(i, j, ijdis)) ++ sum;
 86             else if(IinJ(j, i, ijdis));
 87             else if(IcutJ(i, j, ijdis)) AddSeg(xlij + asimj, xlji - asimj);
 88             else
 89             {
 90                 AddSeg(xlji - asipj, xlji - asimj);
 91                 AddSeg(xlij + asimj, xlij + asipj);
 92             }
 93         }
 94     }
 95     qsort(s, stp, sizeof(Seg), comp);
 96     //  sort(s, s + stp, comp);
 97     for(i = 0; i < stp; ++ i)
 98     {
 99         if(s[i].se) ++ sum;
100         else -- sum;
101         if(sum > ans) ans = sum;
102     }
103     if(sum > ans) ans = sum;
104 }
105 int main()
106 {
107     int t, ca, i;
108     for(scanf("%d", &t), ca = 1; ca <= t; ++ ca)
109     {
110         scanf("%d", &n);
111         for(i = 0; i < n; ++ i)
112             scanf("%lf%lf%lf", &ra[i].x, &ra[i].y, &ra[i].r);
113         for(i = ans = 0; i < n; ++ i) Cal(i);
114         printf("Case #%d: %d\n", ca, ans + 1);
115     }
116     return 0;
117 }

 

 

 
posted @ 2012-08-27 20:42  CSGrandeur  阅读(653)  评论(0)  编辑  收藏  举报