终于在抗不住睡觉以前把这道题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;
}
#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;
}