結果
問題 |
No.691 E869120 and Constructing Array 5
|
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
#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; }