結果

問題 No.691 E869120 and Constructing Array 5
ユーザー mcfx0mcfx0
提出日時 2019-06-29 00:44:21
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 333 ms / 1,000 ms
コード長 4,412 bytes
コンパイル時間 1,733 ms
コンパイル使用メモリ 171,180 KB
実行使用メモリ 4,504 KB
最終ジャッジ日時 2023-09-14 22:56:03
合計ジャッジ時間 12,444 ms
ジャッジサーバーID
(参考情報)
judge13 / judge15
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 333 ms
4,376 KB
testcase_01 AC 329 ms
4,380 KB
testcase_02 AC 316 ms
4,376 KB
testcase_03 AC 316 ms
4,376 KB
testcase_04 AC 316 ms
4,380 KB
testcase_05 AC 317 ms
4,376 KB
testcase_06 AC 329 ms
4,380 KB
testcase_07 AC 325 ms
4,376 KB
testcase_08 AC 321 ms
4,380 KB
testcase_09 AC 317 ms
4,380 KB
testcase_10 AC 317 ms
4,380 KB
testcase_11 AC 320 ms
4,380 KB
testcase_12 AC 321 ms
4,380 KB
testcase_13 AC 316 ms
4,380 KB
testcase_14 AC 323 ms
4,380 KB
testcase_15 AC 324 ms
4,376 KB
testcase_16 AC 313 ms
4,380 KB
testcase_17 AC 315 ms
4,376 KB
testcase_18 AC 319 ms
4,380 KB
testcase_19 AC 322 ms
4,380 KB
testcase_20 AC 317 ms
4,504 KB
testcase_21 AC 322 ms
4,380 KB
testcase_22 AC 329 ms
4,380 KB
testcase_23 AC 316 ms
4,376 KB
testcase_24 AC 323 ms
4,380 KB
testcase_25 AC 323 ms
4,380 KB
testcase_26 AC 320 ms
4,380 KB
testcase_27 AC 3 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>

typedef unsigned int uint;
typedef long long ll;
typedef unsigned long long ull;
typedef double lf;
typedef long double llf;
typedef std::pair<int,int> pii;

#define xx first
#define yy second

template<typename T> inline T max(T a,T b){return a>b?a:b;}
template<typename T> inline T min(T a,T b){return a<b?a:b;}
template<typename T> inline T abs(T a){return a>0?a:-a;}
template<typename T> inline bool repr(T &a,T b){return a<b?a=b,1:0;}
template<typename T> inline bool repl(T &a,T b){return a>b?a=b,1:0;}
template<typename T> inline T gcd(T a,T b){T t;if(a<b){while(a){t=a;a=b%a;b=t;}return b;}else{while(b){t=b;b=a%b;a=t;}return a;}}
template<typename T> inline T sqr(T x){return x*x;}
#define mp(a,b) std::make_pair(a,b)
#define pb push_back
#define I __attribute__((always_inline))inline
#define mset(a,b) memset(a,b,sizeof(a))
#define mcpy(a,b) memcpy(a,b,sizeof(a))

#define fo0(i,n) for(int i=0,i##end=n;i<i##end;i++)
#define fo1(i,n) for(int i=1,i##end=n;i<=i##end;i++)
#define fo(i,a,b) for(int i=a,i##end=b;i<=i##end;i++)
#define fd0(i,n) for(int i=(n)-1;~i;i--)
#define fd1(i,n) for(int i=n;i;i--)
#define fd(i,a,b) for(int i=a,i##end=b;i>=i##end;i--)
#define foe(i,x)for(__typeof((x).end())i=(x).begin();i!=(x).end();++i)
#define fre(i,x)for(__typeof((x).rend())i=(x).rbegin();i!=(x).rend();++i)

struct Cg{I char operator()(){return getchar();}};
struct Cp{I void operator()(char x){putchar(x);}};
#define OP operator
#define RT return *this;
#define UC unsigned char
#define RX x=0;UC t=P();while((t<'0'||t>'9')&&t!='-')t=P();bool f=0;\
if(t=='-')t=P(),f=1;x=t-'0';for(t=P();t>='0'&&t<='9';t=P())x=x*10+t-'0'
#define RL if(t=='.'){lf u=0.1;for(t=P();t>='0'&&t<='9';t=P(),u*=0.1)x+=u*(t-'0');}if(f)x=-x
#define RU x=0;UC t=P();while(t<'0'||t>'9')t=P();x=t-'0';for(t=P();t>='0'&&t<='9';t=P())x=x*10+t-'0'
#define TR *this,x;return x;
I bool IS(char x){return x==10||x==13||x==' ';}template<typename T>struct Fr{T P;I Fr&OP,(int&x)
{RX;if(f)x=-x;RT}I OP int(){int x;TR}I Fr&OP,(ll &x){RX;if(f)x=-x;RT}I OP ll(){ll x;TR}I Fr&OP,(char&x)
{for(x=P();IS(x);x=P());RT}I OP char(){char x;TR}I Fr&OP,(char*x){char t=P();for(;IS(t);t=P());if(~t){for(;!IS
(t)&&~t;t=P())*x++=t;}*x++=0;RT}I Fr&OP,(lf&x){RX;RL;RT}I OP lf(){lf x;TR}I Fr&OP,(llf&x){RX;RL;RT}I OP llf()
{llf x;TR}I Fr&OP,(uint&x){RU;RT}I OP uint(){uint x;TR}I Fr&OP,(ull&x){RU;RT}I OP ull(){ull x;TR}};Fr<Cg>in;
#define WI(S) if(x){if(x<0)P('-'),x=-x;UC s[S],c=0;while(x)s[c++]=x%10+'0',x/=10;while(c--)P(s[c]);}else P('0')
#define WL if(y){lf t=0.5;for(int i=y;i--;)t*=0.1;if(x>=0)x+=t;else x-=t,P('-');*this,(ll)(abs(x));P('.');if(x<0)\
x=-x;while(y--){x*=10;x-=floor(x*0.1)*10;P(((int)x)%10+'0');}}else if(x>=0)*this,(ll)(x+0.5);else *this,(ll)(x-0.5);
#define WU(S) if(x){UC s[S],c=0;while(x)s[c++]=x%10+'0',x/=10;while(c--)P(s[c]);}else P('0')
template<typename T>struct Fw{T P;I Fw&OP,(int x){WI(10);RT}I Fw&OP()(int x){WI(10);RT}I Fw&OP,(uint x){WU(10);RT}
I Fw&OP()(uint x){WU(10);RT}I Fw&OP,(ll x){WI(19);RT}I Fw&OP()(ll x){WI(19);RT}I Fw&OP,(ull x){WU(20);RT}I Fw&OP()
(ull x){WU(20);RT}I Fw&OP,(char x){P(x);RT}I Fw&OP()(char x){P(x);RT}I Fw&OP,(const char*x){while(*x)P(*x++);RT}
I Fw&OP()(const char*x){while(*x)P(*x++);RT}I Fw&OP()(lf x,int y){WL;RT}I Fw&OP()(llf x,int y){WL;RT}};Fw<Cp>out;

const int N=3;

std::map<llf,std::vector<int>>u[N];
llf req[N],eps[N];

inline std::vector<int> merge(const std::vector<int>&a,const std::vector<int>&b)
{
	std::vector<int>r;
	r.resize(a.size()+b.size());
	fo0(i,a.size())r[i]=a[i];
	fo0(i,b.size())r[i+a.size()]=b[i];
	return r;
}

inline void add(llf x,const std::vector<int>&s,int i)
{
	//printf("add:%.15Lf %d %.15Lf %.15Lf\n",x,i,abs(x-req[i]),eps[i]);
	if(abs(x-req[i])<eps[i])
	{
		u[i][x]=s;
		if(i==N-1)return;
	}else return;
	llf t=req[i+1]-x;
	auto p=u[i].lower_bound(t);
	if(p!=u[i].end())add(p->xx+x,merge(p->yy,s),i+1);
	if(p!=u[i].begin())--p,add(p->xx+x,merge(p->yy,s),i+1);
}

int main()
{
	eps[N-1]=9.9e-11l;
	fd0(i,N-1)eps[i]=sqrtl(eps[i+1]);
	for(int T=in;T--;)
	{
		fo0(i,N)u[i].clear();
		in,req[2];
		fd0(i,N-1)req[i]=req[i+1]/2;
		int t=sqr(req[0]);
		while(!u[N-1].size())
		{
			int v=rand()%t;
			llf sv=sqrtl(v);
			int u=sqr(req[0]-sv);
			add(sv+sqrtl(u),{u,v},0);
			add(sv+sqrtl(u+1),{u+1,v},0);
		}
		//printf("%.15Lf\n",u[N-1].begin()->xx);
		auto&tmp=u[N-1].begin()->yy;
		out,(int)tmp.size();foe(i,tmp)out,' ',*i;out,'\n';
	}
}
0