ZOJ-1041-Transmitters

终于在抗不住睡觉以前把这道题AC了，对测试数据出正确答案以后，我就开始提交，结果WA了4次，崩溃了都快，怎么改都有问题，可是算法推理出来又确实没有问题，最后找呀找，随手把所有的非计数非标志的数据类型改成double。。。你猜怎么着，就这样AC了。。。OJ不能老这么折磨我啊，每次都是出这种问题。

这道题猛一看很变态，尤其一想怎么让那个半圆转起来，直接就晕了，加上那份流传的ZOJ分类上是这样说这道题的“1041 Transmitters 变态题……好烦人……现在绝对没有信心再做一遍了…… ”，让人心里平添几分恐惧。。。不过它的AC率还是蛮高的，其实算起来困难也并不太大，我直接遍历就过了，连记录已计算的结果避免重复计算的步骤都没加

主要算法是筛选出半径内的点，然后检验每一点在和其他点的夹角，当满足顺时针小于pi时，计数器加1。。。就是这样。。。

C++ 00:00.01 400K

 

#include<stdio.h>
#include<math.h>

#define MAXN 151
#define pi 3.1415926

class POINT
...{
public:
    double x,y;
};

POINT e[MAXN];//检验点
int n;//检验点数
POINT o;//原点f
double r;//半径
bool check[MAXN];

double dist(POINT p1) // 返回两点之间欧氏距离
...{
    return(sqrt(double((p1.x-o.x)*(p1.x-o.x)+(p1.y-o.y)*(p1.y-o.y))));
}

/**//* 返回顶角在o点，起始边为os，终止边为oe的夹角(单位：弧度)
角度小于pi，返回正值
角度大于pi，返回负值
*/
double angle(POINT s,POINT e)
...{
    double cosfi,fi,norm;
    double dsx = s.x - o.x;
    double dsy = s.y - o.y;
    double dex = e.x - o.x;
    double dey = e.y - o.y;
    
    cosfi=dsx*dex+dsy*dey;
    norm=(dsx*dsx+dsy*dsy)*(dex*dex+dey*dey);
    cosfi /= sqrt( norm );

    if (cosfi >= 1.0 ) 
        return 0;
    else if (cosfi <= -1.0 ) 
        return 0;
    
    fi=acos(cosfi);

    if (dsx*dey-dsy*dex>0) return fi; // 说明矢量os 在矢量 oe的顺时针方向
    return -fi;

}

int solve()
...{
    int i,j,count1,key=0;
    double temp,dis;
    for(i=1;i<=n;i++)
    ...{
        dis = dist(e[i]);
        if(dis<=r)
            check[i] = 1;
        else
            check[i] = 0;
    }
    
    for(i=1;i<=n;i++)
    ...{
        if(check[i])
        ...{                        
            count1=1;
            for(j=1;j<=n;j++)
            ...{
                if(i!=j && check[j])
                ...{
                    temp = angle(e[i],e[j]);
                    if(temp>=0.0)
                        count1++;
                }
            }            
            if(count1>key)
                key = count1;

        }
    }
    return key;
}

int main()
...{    
    int i;
//    freopen("1041.txt","r",stdin);
    scanf("%lf %lf %lf",&o.x,&o.y,&r);    
    while(r>0.0)
    ...{
        scanf("%d",&n);
        for(i=1;i<=n;i++)
        ...{
            scanf("%lf %lf",&e[i].x,&e[i].y);
        }
        for(i=0;i<MAXN;i++)
            check[i] = 0;
        printf("%d/n",solve());
        scanf("%lf %lf %lf",&o.x,&o.y,&r);
    }
//    fclose(stdin);
    return 0;
}