19_05_01校內訓練[polygon]
阿新 • • 發佈:2019-05-01
代碼 its eve pre open ++ 得到 ima pan
題意
把一個邊長為1的正n邊形放到一個正m邊形中,要求m邊形完全覆蓋n邊形,可以有交點,並且中心重合。求正m邊形的最小邊長,至少精確到6位。要求logn計算。
思考
先考慮m|n的情況。
我們知道,正m邊形的邊長與可行區域(即可以完全覆蓋的那些角度)形成單射,當且僅當所有可行區域都成為可數的點時,答案最優。(可以理解為再縮小一點就無解了)
這樣不難證明,把正n邊形的幾條邊剛好卡在正m邊形上是最優的。如n=8,m=4:
這時正m邊形的邊長是容易計算的。相信大家都會初中數學。
這樣再考慮一般情況。由於是中心重合,正n邊形旋轉2π/m度後仍然是能被覆蓋的。
在所有可行的旋轉過程中,將最外圈的點連起來,仍然形成一個正多邊形,且邊數為lcm(n,m)。
例如,n=4,m=6:
用紫線圍出來的正12邊形即為正方形得到的結果。
至於正確性,在於所有的可行區域都是單點。
這樣一來,就可以直接轉化為上一個問題。公式認真推即可。
代碼
1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long int ll; 4 const double pi=acos(-1); 5 ll n,m; 6 ll gcd(ll x,ll y) 7 { 8 return x%y==0?y:gcd(y,x%y); 9View Code} 10 ll lcm(ll x,ll y) 11 { 12 return x/gcd(x,y)*y; 13 } 14 double solve(ll n,ll m) 15 { 16 double len=1/(2*tan(pi/n)); 17 double th=(n/m)*pi/n; 18 return tan(th)*len*2; 19 } 20 int main() 21 { 22 ios::sync_with_stdio(false); 23 cin>>n>>m; 24 double len=1/(2*sin(pi/n)); 25 n=lcm(n,m); 26 double a=sin(pi/n)*len*2; 27 cout<<fixed<<setprecision(9)<<solve(n,m)*a<<endl; 28 return 0; 29 }
19_05_01校內訓練[polygon]