結果

問題 No.691 E869120 and Constructing Array 5
ユーザー square1001
提出日時 2018-05-18 23:33:53
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 158 ms / 1,000 ms
コード長 1,249 bytes
コンパイル時間 874 ms
コンパイル使用メモリ 84,220 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-28 14:48:28
合計ジャッジ時間 5,164 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cmath>
#include <vector>
#include <iomanip>
#include <iostream>
using namespace std;
int q; long double 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++) {
			long double bs = sqrt((long double)(j));
			long double rem = sqrt((long double)(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((long double)(j - m)) - sqrt((long double)(j + m))) l = m;
					else r = m;
				}
				rem -= 2 * bs - sqrt((long double)(j - l)) - sqrt((long double)(j + l));
				delta[k] = l;
			}
			if (rem < 1.0e-10L) {
				vector<long double> ret(2 * sz);
				for (int k = 0; k < sz; k++) {
					ret[2 * k] = j - delta[k];
					ret[2 * k + 1] = j + delta[k];
				}
				long double 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