結果

問題 No.2954 Calculation of Exponentiation
コンテスト
ユーザー tobisatis
提出日時 2024-11-08 22:23:07
言語 C#
(.NET 8.0.404)
結果
AC  
実行時間 55 ms / 2,000 ms
コード長 3,079 bytes
コンパイル時間 9,336 ms
コンパイル使用メモリ 167,452 KB
実行使用メモリ 188,248 KB
最終ジャッジ日時 2024-11-08 22:23:21
合計ジャッジ時間 10,175 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (89 ms)。
MSBuild のバージョン 17.9.6+a4ecab324 (.NET)
  main -> /home/judge/data/code/bin/Release/net8.0/main.dll
  main -> /home/judge/data/code/bin/Release/net8.0/publish/

ソースコード

diff # 

namespace AtCoder;

#nullable enable

using System.Numerics;

static class Extensions
{
    public static T[] Repeat<T>(this int time, Func<T> F) => Enumerable.Range(0, time).Select(_ => F()).ToArray();

    public static long Root(this long x, int root = 2)
    {
        var limits = new long[] { 0, 0, 3037000500L, 2097152, 55109, 6209, 1449, 512, 235, 128, 79, 53, 39 };
        if (x < 0 || root < 2 || root >= limits.Length) throw new Exception("not supported");
        var pass = 0L;
        var fail = limits[root];
        while (fail - pass > 1)
        {
            var middle = (fail + pass) >> 1;
            var next = 1L;
            for (var i = 0; i < root; i++) next *= middle;
            if (next <= x) pass = middle; else fail = middle;
        }
        return pass;
    }
}

class AtCoder
{
    object? Solve()
    {
        var d = 10000L;
        static long GCD(long x, long y) => y == 0 ? x : GCD(y, x % y);
        var ad = Input<decimal>();
        var bd = Input<decimal>();
        if (bd == 0) return true;
        var ap = (long)(ad * d);
        var aq = d;
        if (bd < 0) (ap, aq) = (aq, ap);
        if (aq != GCD(ap, aq)) return false;
        ap /= aq;
        var bp = (long)(bd * d);
        var div = GCD(Math.Abs(bp), d);
        var bq = d / div;
        if (bq == 1) return true;
        if (bq > 1000) return false;
        if (bq < 10)
        {
            var r = Extensions.Root(ap, (int)bq);
            var o = 1L;
            for (var i = 0; i < bq; i++) o *= r;
            return o == ap;
        }
        {
            var r = 1L;
            while (true)
            {
                Int128 o = 1L;
                for (var i = 0; i < bq; i++)
                {
                    if (o * r > long.MaxValue) return false;
                    o *= r;
                }
                if ((long)o == ap) return true;
                r++;
            }
        }
    }

    public static void Main() => new AtCoder().Run();
    public void Run()
    {
        var res = Solve();
        if (res != null)
        {
            if (res is bool yes) res = yes ? "Yes" : "No";
            sw.WriteLine(res);
        }
        sw.Flush();
    }

    string[] input = Array.Empty<string>();
    int iter = 0;
    readonly StreamWriter sw = new(Console.OpenStandardOutput()) { AutoFlush = false };

    string String()
    {
        while (iter >= input.Length) (input, iter) = (Console.ReadLine()!.Trim().Split(' '), 0);
        return input[iter++];
    }
    T Input<T>() where T : IParsable<T> => T.Parse(String(), null);
    int Int() => Input<int>();
    void Out(object? x, string? separator = null)
    {
        separator ??= Environment.NewLine;
        if (x is System.Collections.IEnumerable obj and not string)
        {
            var firstLine = true;
            foreach (var item in obj)
            {
                if (!firstLine) sw.Write(separator);
                firstLine = false;
                sw.Write(item);
            }
        }
        else sw.Write(x);
        sw.WriteLine();
    }
}
0