結果

問題 No.3527 Minimum Abs Sum
コンテスト
ユーザー pengin_2000
提出日時 2026-05-04 23:49:27
言語 C
(gcc 15.2.0)
コンパイル:
gcc-15 -O2 -DONLINE_JUDGE -o a.out _filename_ -lm
実行:
./a.out
結果
AC  
実行時間 73 ms / 2,000 ms
コード長 1,917 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 218 ms
コンパイル使用メモリ 39,880 KB
実行使用メモリ 8,188 KB
最終ジャッジ日時 2026-05-04 23:49:44
合計ジャッジ時間 3,989 ms
ジャッジサーバーID
（参考情報） 		judge1_1 / judge2_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<stdio.h>
long long int modpow(long long int a, long long int n, long long int p)
{
	long long int res = 1;
	for (; n > 0; n /= 2, a = a * a % p)
		if (n % 2 > 0)
			res = res * a % p;
	return res;
}
long long int A[200005], B[200005];
long long int h[200005], l;
long long int comp_h(long long int a, long long int b)
{
	a = h[a];
	b = h[b];
	if (B[a] * A[b] > B[b] * A[a])
		return 1;
	else
		return -1;
}
void swap_h(long long int a, long long int b)
{
	long long int f = h[a];
	h[a] = h[b];
	h[b] = f;
	return;
}
void push(long long int ne)
{
	h[l] = ne;
	long long int p = l++;
	for (; p > 0; p = (p - 1) / 2)
		if (comp_h((p - 1) / 2, p) > 0)
			swap_h((p - 1) / 2, p);
	return;
}
long long int pop()
{
	swap_h(0, --l);
	long long int p = 0;
	for (;;)
	{
		if (2 * p + 2 < l)
		{
			if (comp_h(2 * p + 1, 2 * p + 2) > 0)
			{
				if (comp_h(p, 2 * p + 2) > 0)
					swap_h(p, 2 * p + 2);
				p = 2 * p + 2;
			}
			else
			{
				if (comp_h(p, 2 * p + 1) > 0)
					swap_h(p, 2 * p + 1);
				p = 2 * p + 1;
			}
		}
		else if (2 * p + 1 < l)
		{
			if (comp_h(p, 2 * p + 1) > 0)
				swap_h(p, 2 * p + 1);
			p = 2 * p + 1;
		}
		else
			break;
	}
	return h[l];
}
long long int id[200005];
int main()
{
	long long int n;
	scanf("%lld", &n);
	long long int i, j;
	for (i = 0; i < n; i++)
		scanf("%lld", &A[i]);
	for (i = 0; i < n; i++)
		scanf("%lld", &B[i]);
	const long long int p = 1000000007;
	for (j = i = 0; i < n; i++)
	{
		if (A[i] == 0)
			continue;
		if (A[i] < 0)
		{
			A[i] *= -1;
			B[i] *= -1;
		}
		A[j] = A[i];
		B[j] = B[i];
		j++;
	}
	n = j;
	l = 0;
	for (i = 0; i < n; i++)
		push(i);
	for (i = 0; i < n; i++)
		id[i] = pop();
	long long int sum = 0;
	for (i = 0; i < n; i++)
		sum += A[i];
	for (i = 0;; i++)
	{
		sum -= 2 * A[id[i]];
		if (sum <= 0)
			break;
	}
	long long int ans = B[id[i]] * modpow(A[id[i]], p - 2, p) % p;
	ans = (ans + p) % p;
	printf("%lld\n", ans);
	return 0;
}
0