結果

問題 No.691 E869120 and Constructing Array 5
ユーザー vjudge1
提出日時 2025-09-21 17:06:58
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 158 ms / 1,000 ms
コード長 987 bytes
コンパイル時間 1,780 ms
コンパイル使用メモリ 167,764 KB
実行使用メモリ 7,716 KB
最終ジャッジ日時 2025-09-21 17:07:05
合計ジャッジ時間 6,465 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 1
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
using ld=long double;
int q;ld x;
int main() {
	cin>>q;
	for(int i=0;i<q;i++) {
		cin>>x;
		int sz=10;
		int lim=x*x /(4*sz*sz)+1.01L;
		for(int j=lim;;j++) {
			ld bs=sqrt((ld)(j));
			ld rem=sqrt((ld)(j))*(2*sz)-x;
			vector<int> delta(sz);
			for(int k=0;k<sz;k++) {
				int l=0,r=j;
				while(r-l>1) {
					int m=(l+r)>>1;
					if(rem>2*bs-sqrt((ld)(j-m))-sqrt((ld)(j+m))) l=m;
					else r=m;
				}
				rem -= 2*bs-sqrt((ld)(j-l))-sqrt((ld)(j+l));
				delta[k]=l;
			}
			if(rem<1.0e-10L) {
				vector<ld> ret(2*sz);
				for(int k=0;k<sz;k++) {
					ret[2*k]=j-delta[k];
					ret[2*k+1]=j+delta[k];
				}
				ld val=0.0;
				for(int k=0;k<ret.size();k++) {
					val += sqrt(ret[k]);
				}
				cerr << fixed << setprecision(15) << j << ' ' << lim << ' ' << val << ' ' << abs(val-x) << endl;
				cout << ret.size();
				for(int k=0;k<ret.size();k++) {
					cout << ' ' << ret[k];
				}
				cout << '\n';
				break;
			}
		}
	}
	return 0;
}
0