题目地址:
思路:
列出方程 2*R*x=L‘
2*R*sin(x)=L
两式相除即得 x/sin(x) = 1+n*c 前提x!=0 就是 n*c>0
答案就是 L/2* (1/sin(x)-1/tan(x)) 三角函数化简为 L/2*tan(x/2) x 在0~ PI/2 于是关于x是单增的, 反函数也是单增的,x/sin(x) 也是单增的
一开始的思路是先用 二分把 x/sin(x) = 1+n*c 里面的x解出来,然后带进 L/2*tan(x/2) 计算 虽然样例过了,但是还是wa,看来是精度问题。
所以直接二分答案,ans= L/2*tan(x/2) => x=2*arctan(ans*2/L) 这个关于ans是单增的,带进 x/sin(x) 也是单增的 可以二分了
代码:
#include#include #include const long double PI=acos(-1.0);long double L,n,c;long double f(long double x){ double xx=2*atan(2*x/L); return xx/sin(xx);}using namespace std;int main(){ long double l,r,mid; while(cin>>L>>n>>c) { if(L==-1) break; l=0; //r=PI/2; r=L; double ans; if(n*c>0) { while(r-l>1e-4) { mid=(l+r)/2; if(f(mid)>1+n*c) r=mid; else l=mid; } //ans=1.0/sin(l)-1.0/tan(l); ans=l; } else ans=0; printf("%.3lf\n",ans); }}