結果

問題 No.3074 Divide Points Fairly
ユーザー nouka28
提出日時 2025-03-28 23:28:28
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 1,572 bytes
コンパイル時間 5,967 ms
コンパイル使用メモリ 335,080 KB
実行使用メモリ 7,328 KB
最終ジャッジ日時 2025-03-28 23:28:49
合計ジャッジ時間 19,380 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30 RE * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

#include<atcoder/all>
using namespace atcoder;
using mint=atcoder::modint998244353;

#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#define rep(i,n) for(int i=0;i<(n);i++)
#define rng(i,l,r) for(int i=(l);i<(r);i++)
#define rrep(i,n) for(int i=(n)-1;i>=0;i--)
#define rrng(i,l,r) for(int i=(r)-1;i>=(l);i--)

#define ALL(x) (x).begin(),(x).end()

#define int long long
#define fi first
#define se second

struct fast_io{fast_io(){cin.tie(nullptr)->sync_with_stdio(false);}}_;

void solve();
signed main(){
	int T=1;
	// cin>>T;
	while(T--)solve();
}

//乱数系
random_device rnd; 
mt19937 mt(rnd());
const long long MT_MAX=1e18;
uniform_int_distribution<long long> rd(0,MT_MAX);
double randd(){
    return 1.0*rd(mt)/MT_MAX;
}
long long randint(long long a,long long b){
    // [a,b]の乱数を生成
    return a+rd(mt)%(b-a+1);
}


void solve(){
	int N;cin>>N;
	vector<pair<int,int>> v(2*N);
	for(auto&&[x,y]:v)cin>>x>>y;

	auto get=[&](double deg)->pair<int,int> {
		double a=tan(deg/180*acos(-1)),b=1;
		// cout<<"deg : "<<a<<" "<<b<<endl;
		int p=100000/max(abs(a),abs(b));
		int A=a*p,B=b*p;
		return {A,B};
	};
	
	for(int q=-1000;q<=1000;q+=200){

		rep(i,500){
			auto[A,B]=get(double(i)*180/500.0);
			// cout<<i<<" : "<<A<<" "<<B<<endl;
			int cnt1=0,cnt2=0;
			for(auto&&[x,y]:v){
				if(A*x-B*y-A*q<0)cnt1++;
				if(A*x-B*y-A*q==0)cnt2++;
			}

			if(cnt2)continue;
			// cout<<i<<" : "<<cnt1<<endl;
			if(cnt1==N){
				cout<<A<<" "<<-B<<" "<<-A*q<<endl;
				return;
			}
		}
	}

	assert(0);
}
0