POJ1905解题报告

题目链接

 

题目的那点在于如何求解角度，这里采用的方法是二分求解：

一般地，对于函数f(x),如果存在实数c,当x=c时f(c)=0,那么把x=c叫做函数f(x)的零点。 

  解方程即要求f(x)的所有零点。 

  先找到a、b，使f(a)，f(b)异号，说明在区间(a,b)内一定有零点，然后求f[(a+b)/2], 

  现在假设f(a)<0,f(b)>0,a<b 

  ①如果f[(a+b)/2]=0，该点就是零点， 

  如果f[(a+b)/2]<0,则在区间（(a+b)/2，b)内有零点，(a+b)/2=>a，从①开始继续使用

  中点函数值判断。

  如果f[(a+b)/2]>0，则在区间(a,(a+b)/2)内有零点，(a+b)/2=>b，从①开始继续使用

  中点函数值判断。

  这样就可以不断接近零点。

  通过每次把f(x)的零点所在小区间收缩一半的方法，使区间的两个端点逐步迫近函数的零点，以求得零点的近似值，这种方法叫做二分法。

（以上摘自百度）



因此，可以用二分迭代的方式求解该题的角度。

当两点之间的距离小于一个阈值时，停止迭代计算，输出结果即可。

代码如下：

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const double MyError = 10e-12;
const double PI = atan(1.0)*4;
double alpha;
double f(const double &x)
{
	return alpha*sin(x)-x;
}

double GetSita(double b,double e)
{
	//recursion
// 	double now = e-b;
// 	if(now< MyError && now >0)
// 		return e;
// 	else now = (b+e)/2;
// 	if(f(now)>0.0)
// 		return GetSita(now,e);
// 	else return GetSita(b,now);
	double now = e-b;
	while (now >MyError)
	{
		now = (b+e)/2;
		if( f(now) > 0.0)
			b = now;
		else e = now;
		now = e-b;
	}
	return e;
}

int main()
{
 	double n,L,c;
   	while(cin>>L>>n>>c)
   	{
   		if(n == -1 && L == -1 && c == -1)
   			break;
   
		if(n<MyError || c< MyError)
		{
  			cout<<"0.000"<<endl;
			continue;
		}
   		alpha =1.0+ n*c;
  		double sita = GetSita(double(1e-5),10.0);
  		printf("%.3f/n",L/2/tan( (PI-sita)/2 ));
   	}
	return 0;
}
 



这里还要说的一点就是不能走意识流的路线，看着题目中的n和L是整数就以为它们真的是整数？F**K，那是POJ逗你玩。。。白WA了两次。。。