ZZULIOJ 1882: 蛤玮的魔法 （计算几何）

1882: 蛤玮的魔法

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 313   Solved: 95

Submit Status Web Board

Description

为了成为魔法少女,蛤玮正在学习画魔法阵,他首先画了一个正n边形,查阅魔法书后蛤玮发现书上要求魔法阵的面积不能超过L,他很头疼,因为用尺规作这个正n边形花了他好大经历,他不想重新画一边,于是他想了个好主意,把每条边的中点依次连起来,就能得到一个缩小的正n边行.现在蛤玮想知道他需要按着个方法缩小多少次才能合乎魔法书上的要求.

Input

T(1<=T<=100),表示数据组数.

每组数据三个整数n(3<=n<=10),a(1<=a<=100),L(1<=L<=1000),其中a表示蛤玮画的正n边行的边长,n,L如题中描述.

数据保证[L-1e-5,L+1e-5]内答案唯一.

Output

每组数据输出一个整数,表示蛤玮操作的次数.

Sample Input

1
4 2 3

Sample Output

1

HINT

Source

按照题目要进行一步步计算，要用到三角函数正弦定理 a/sinA=b/sinB=c/sinC 求边，以及用三角函数求三角形面积s=0.5*a*b*sinC;

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#define PI acos(-1.0)
using namespace std;
int main()
{
	int t;
	double n,a,l,l1,x1,x2,x3,x4,s1,a1,s;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lf%lf%lf",&n,&a,&l);
		x1=(PI/n)*2;
		x2=(PI-x1)/2;
		x3=x2*2;
		x4=(PI-x3)/2;
		l1=a*sin(x2)/sin(x1);
		s=0.5*sin(x1)*l1*l1*n;
		int cnt=0;
		while(s-l>1e-6)
		{
			a1=a/2;
			s1=0.5*a1*a1*sin(x3);
			s=s-s1*n;
			a=a1*sin(x3)/sin(x4);
			cnt++;
		}
		printf("%d\n",cnt);
	}
	return 0;
}