結果

問題 No.2954 Calculation of Exponentiation
ユーザー ygussanyygussany
提出日時 2024-11-08 21:50:33
言語 C
(gcc 12.3.0)
結果
WA  
実行時間 -
コード長 2,079 bytes
コンパイル時間 233 ms
コンパイル使用メモリ 32,508 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-08 21:50:36
合計ジャッジ時間 1,071 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 1 ms
5,248 KB
testcase_02 AC 1 ms
5,248 KB
testcase_03 AC 1 ms
5,248 KB
testcase_04 AC 1 ms
5,248 KB
testcase_05 AC 1 ms
5,248 KB
testcase_06 AC 1 ms
5,248 KB
testcase_07 AC 1 ms
5,248 KB
testcase_08 AC 1 ms
5,248 KB
testcase_09 AC 1 ms
5,248 KB
testcase_10 AC 1 ms
5,248 KB
testcase_11 AC 1 ms
5,248 KB
testcase_12 AC 0 ms
5,248 KB
testcase_13 AC 1 ms
5,248 KB
testcase_14 WA -
testcase_15 AC 1 ms
5,248 KB
testcase_16 AC 1 ms
5,248 KB
testcase_17 AC 1 ms
5,248 KB
testcase_18 AC 1 ms
5,248 KB
testcase_19 AC 1 ms
5,248 KB
testcase_20 WA -
testcase_21 AC 1 ms
5,248 KB
testcase_22 WA -
testcase_23 AC 0 ms
5,248 KB
testcase_24 AC 1 ms
5,248 KB
testcase_25 AC 1 ms
5,248 KB
testcase_26 AC 1 ms
5,248 KB
testcase_27 AC 1 ms
5,248 KB
testcase_28 AC 1 ms
5,248 KB
testcase_29 AC 1 ms
5,248 KB
testcase_30 AC 1 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <stdlib.h>

long long gcd(long long a, long long b)
{
	if (a == 0) return b;
	else return gcd(b % a, a);
}

typedef struct {
	long long numer, denom;
} ratio;

void reduce_ratio(ratio *a)
{
	int g = gcd(abs(a->numer), abs(a->denom));
	a->numer /= g;
	a->denom /= g;
	if (a->denom < 0) {
		a->numer *= -1;
		a->denom *= -1;
	}
}

int is_integer_ratio(ratio a)
{
	reduce_ratio(&a);
	if (a.denom == 1) return 1;
	else return 0;
}

int main()
{
	char A[15], B[15];
	scanf("%s", A);
	scanf("%s", B);
	
	int i;
	ratio a, b;
	for (i = 0, a.numer = 0; A[i] != '.'; i++) a.numer = a.numer * 10 + A[i] - '0';
	for (i++; A[i] != 0; i++) a.numer = a.numer * 10 + A[i] - '0';
	a.denom = 10000;
	for (i = (B[0] == '-')? 1: 0, b.numer = 0; B[i] != '.'; i++) b.numer = b.numer * 10 + B[i] - '0';
	for (i++; B[i] != 0; i++) b.numer = b.numer * 10 + B[i] - '0';
	b.denom = 10000;
	if (B[0] == '-') {
		a.denom ^= a.numer;
		a.numer ^= a.denom;
		a.denom ^= a.numer;
	}
	reduce_ratio(&a);
	reduce_ratio(&b);
	if (b.numer == 0) {
		printf("Yes\n");
		return 0;
	}

	__int128 l, r, m;
	while (b.denom % 2 == 0) {
		b.denom /= 2;
		l = 1;
		r = a.numer;
		while (l < r) {
			m = (l + r) / 2;
			if (m * m < a.numer) l = m + 1;
			else r = m;
		}
		if (l * l == a.numer) a.numer = l;
		else {
			printf("No\n");
			return 0;
		}
		
		l = 1;
		r = a.denom;
		while (l < r) {
			m = (l + r) / 2;
			if (m * m < a.denom) l = m + 1;
			else r = m;
		}
		if (l * l == a.denom) a.denom = l;
		else {
			printf("No\n");
			return 0;
		}
	}
	while (b.denom % 5 == 0) {
		b.denom /= 5;
		l = 1;
		r = a.numer;
		while (l < r) {
			m = (l + r) / 2;
			if (m * m * m * m * m < a.numer) l = m + 1;
			else r = m;
		}
		if (l * l * l * l * l == a.numer) a.numer = l;
		else {
			printf("No\n");
			return 0;
		}
		
		l = 1;
		r = a.denom;
		while (l < r) {
			m = (l + r) / 2;
			if (m * m * m * m * m < a.denom) l = m + 1;
			else r = m;
		}
		if (l * l * l * l * l == a.denom) a.denom = l;
		else {
			printf("No\n");
			return 0;
		}
	}
	printf("Yes\n");
	fflush(stdout);
	return 0;
}
0