結果

問題 No.2555 Intriguing Triangle
ユーザー 👑 binapbinap
提出日時 2023-12-01 01:26:33
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 2,219 bytes
コンパイル時間 4,432 ms
コンパイル使用メモリ 262,884 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-26 15:27:11
合計ジャッジ時間 5,580 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#include<atcoder/all>
#define rep(i,n) for(int i=0;i<n;i++)
using namespace std;
using namespace atcoder;
typedef long long ll;
typedef vector<int> vi;
typedef vector<long long> vl;
typedef vector<vector<int>> vvi;
typedef vector<vector<long long>> vvl;
typedef long double ld;
typedef pair<int, int> P;

ostream& operator<<(ostream& os, const modint& a) {os << a.val(); return os;}
template <int m> ostream& operator<<(ostream& os, const static_modint<m>& a) {os << a.val(); return os;}
template<typename T> istream& operator>>(istream& is, vector<T>& v){int n = v.size(); assert(n > 0); rep(i, n) is >> v[i]; return is;}
template<typename U, typename T> ostream& operator<<(ostream& os, const pair<U, T>& p){os << p.first << ' ' << p.second; return os;}
template<typename T> ostream& operator<<(ostream& os, const vector<T>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : " "); return os;}
template <typename T> ostream& operator<<(ostream& os, const vector<vector<T>>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : ""); return os;}

const long double delta = 0.000000001;

int main(){
	int a, b, c;
	cin >> a >> b >> c;
	for(int x = 1; x <= 200; x++){
		for(int y = 1; y <= 200; y++){
			if(x + a + y <= abs(b - c) or x + a + y >= b + c) continue;
			long double s = (long double)((x + a + y) + b + c) / 2;
			long double area = sqrtl(s * (s - (x + a + y)) * (s - b) * (s - c));
			
			long double sin_beta = 2 * area / (x + a + y) / b;
			long double cos_beta = sqrtl(1 - sin_beta * sin_beta);
			if(c * c > b * b + (x + a + y) * (x + a + y)) cos_beta *= -1;
			
			long double sin_gamma = 2 * area / (x + a + y) / c;
			long double cos_gamma = sqrtl(1 - sin_gamma * sin_gamma);
			if(b * b > c * c + (x + a + y) * (x + a + y)) cos_gamma *= -1;
			
			long double p = b * b + (x + a) * (x + a) - 2 * b * (x + a) * cos_beta;
			long double q = c * c + (y + a) * (y + a) - 2 * c * (y + a) * cos_gamma;
			
			long double d = (x + a) * (x + a) * sin_beta * sin_beta / p - (y + a) * (y + a) * sin_gamma * sin_gamma / q;
			if(abs(d) < delta){
				cout << "Yes\n";
//				cout << x << ' ' << y << "\n";
				return 0;
			}
		}
	}
	cout << "No\n";
	return 0;
}
0