結果

問題 No.1365 [Cherry 1st Tune] Whose Fault?
ユーザー publfl
提出日時 2021-01-22 21:52:06
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 342 ms / 2,000 ms
コード長 1,373 bytes
コンパイル時間 276 ms
コンパイル使用メモリ 34,176 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-29 09:35:02
合計ジャッジ時間 17,825 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 46
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>

int next[100010];
int find(int k)
{
	if(next[k]==k) return k;
	else return next[k] = find(next[k]);
}

int A[100010],B[100010];
double P[100010];
int check[100010];

int main()
{
	int a;
	scanf("%d",&a);
	for(int i=1;i<=a;i++) scanf("%d",&A[i]);
	for(int i=1;i<=a;i++) scanf("%d",&B[i]);
	for(int i=1;i<=a;i++)
	{
		int c;
		scanf("%d",&c);
		if(c==0) check[i] = 1;
		else P[i] = (double)1/(2*c);
	}
	
	double sum = 0;
	for(int i=1;i<=a;i++) next[i] = i;
	for(int i=1;i<=a;i++)
	{
		if(check[i]==0)
		{
			long long int lambda = (A[i]-B[i]);
			sum += (lambda*lambda)/(2*P[i]);
		}
	}
	
	int b;
	scanf("%d",&b);
	for(int i=1;i<=b;i++)
	{
		int c,d;
		scanf("%d%d",&c,&d);
		c = find(c);
		d = find(d);
		if(c!=d)
		{
			long long int lambda;
			if(check[c]==1&&check[d]==1);
			else if(check[c]==1)
			{
				lambda = (A[d]-B[d]);
				sum -= (lambda * lambda)/(2*P[d]);
				next[c] = d;
				check[d] = 1;
			}
			else if(check[d]==1)
			{
				lambda = (A[c]-B[c]);
				sum -= (lambda * lambda)/(2*P[c]);
				next[c] = d;
			}
			else
			{
				lambda = (A[c]-B[c]);
				sum -= (lambda * lambda)/(2*P[c]);
				lambda = (A[d]-B[d]);
				sum -= (lambda * lambda)/(2*P[d]);
				
				next[c] = d;
				A[d] += A[c];
				B[d] += B[c];
				P[d] += P[c];
				
				lambda = (A[d]-B[d]);
				sum += (lambda * lambda)/(2*P[d]);
			}
		}
		printf("%.12lf\n",sum);
	}
}
0