結果

問題 No.187 中華風 (Hard)
ユーザー EmKjpEmKjp
提出日時 2020-06-25 20:44:47
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 233 ms / 3,000 ms
コード長 8,283 bytes
コンパイル時間 3,746 ms
コンパイル使用メモリ 111,280 KB
実行使用メモリ 23,284 KB
最終ジャッジ日時 2023-09-16 22:58:41
合計ジャッジ時間 8,160 ms
ジャッジサーバーID
(参考情報)
judge13 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 51 ms
19,512 KB
testcase_01 AC 49 ms
19,156 KB
testcase_02 AC 183 ms
21,216 KB
testcase_03 AC 182 ms
21,180 KB
testcase_04 AC 230 ms
23,284 KB
testcase_05 AC 233 ms
21,204 KB
testcase_06 AC 221 ms
19,308 KB
testcase_07 AC 225 ms
21,440 KB
testcase_08 AC 170 ms
21,236 KB
testcase_09 AC 171 ms
21,340 KB
testcase_10 AC 174 ms
21,220 KB
testcase_11 AC 226 ms
21,200 KB
testcase_12 AC 226 ms
21,152 KB
testcase_13 AC 95 ms
21,284 KB
testcase_14 AC 93 ms
21,196 KB
testcase_15 AC 177 ms
21,288 KB
testcase_16 AC 180 ms
21,416 KB
testcase_17 AC 52 ms
21,188 KB
testcase_18 AC 52 ms
21,212 KB
testcase_19 AC 51 ms
21,232 KB
testcase_20 AC 183 ms
21,204 KB
testcase_21 AC 51 ms
21,224 KB
testcase_22 AC 233 ms
21,340 KB
testcase_23 AC 52 ms
21,208 KB
testcase_24 AC 52 ms
23,272 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff #

using System;
using System.Collections;
using System.Collections.Generic;
using System.Diagnostics;
using System.Globalization;
using System.IO;
using System.Text;
using System.Linq;

using E = System.Linq.Enumerable;

internal partial class Solver {

    public static long Gcd(long s, long t) {
        var r = s % t;
        while (r > 0) {
            s = t;
            t = r;
            r = s % t;
        }
        return t;
    }

    private static long ExtendedGcd(long a, long b, out long x, out long y) {
        long s = a, sx = 1, sy = 0;
        long t = b, tx = 0, ty = 1;
        while (true) {
            var d = s / t;
            var u = s - t * d;
            if (u == 0) break;
            var ux = sx - tx * d;
            var uy = sy - ty * d;
            s = t; sx = tx; sy = ty;
            t = u; tx = ux; ty = uy;
        }
        x = tx; y = ty;
        return t;
    }

    private static long Inverse(long a, long mod) {
        long x = 0, y = 0;
        if (ExtendedGcd(a, mod, out x, out y) == 1) {
            return (x + mod) % mod;
        } else {
            throw new Exception("Invalid inverse " + a + " " + mod);
        }
    }

    /// <summary>
    /// m 同士がすべて互いに素になるように変換する
    /// </summary>
    public static bool PreGarner(long[] b, long[] m) {
        for (int i = 0; i < b.Length; i++) {
            for (int j = 0; j < i; j++) {
                var gcd = Gcd(m[i], m[j]);
                if (b[i] % gcd != b[j] % gcd) return false;
                var gi = m[i] / gcd;
                var lcm = gi * m[j];
                while (true) {
                    var g = Gcd(gi * gi, gcd);
                    if (g == 1) break;
                    gi *= g;
                    gcd /= g;
                };
                m[i] = gi;
                m[j] = lcm / gi;
            }
        }
        for (int i = 0; i < b.Length; i++) {
            b[i] %= m[i];
        }
        return true;
    }

    /// <summary>
    /// X = b[i] (mod m[i]) for all i
    /// を満たす最小の X % mod を返す
    /// m はすべて互いに素でなければならない
    /// </summary>
    public static long Garner(long[] b, long[] m, long mod) {
        // x = c1 + c2 m1 + c3 m1 m2 + c4 m1 m2 m3 ...
        var c = new long[m.Length];
        long sum = 0;
        for (int i = 0; i <= m.Length; i++) {
            long coef = 1; // m1*m2*m3* ...
            var modInternal = i == m.Length ? mod : m[i];
            sum = 0;
            for (int j = 0; j < i; j++) {
                sum += coef * c[j]; sum %= modInternal;
                coef *= m[j]; coef %= modInternal;
            }
            if (i == m.Length) break;
            c[i] = (b[i] - sum) * Inverse(coef, modInternal) % modInternal;
            if (c[i] < 0) c[i] += modInternal;
        }
        return sum;
    }

    public void Run() {
        var n = ni();
        var x = new long[n];
        var y = new long[n];
        var mod = 1000000007;
        for (int i = 0; i < n; i++) {
            x[i] = nl();
            y[i] = nl();
        }

        if (!PreGarner(x, y)) {
            cout.WriteLine(-1);
        } else if (x.All(a => a == 0)) {
            var ans = y.Aggregate((a, b) => a * b % mod);
            cout.WriteLine(ans);
        } else {
            var g = Garner(x, y, mod);
            cout.WriteLine(g);
        }
    }
}

// PREWRITEN CODE BEGINS FROM HERE

static public class StringExtensions {
    static public string JoinToString<T>(this IEnumerable<T> source, string separator = " ") {
        return string.Join(separator, source);
    }
}

internal partial class Solver : Scanner {
    public static void Main() {
#if LOCAL
        byte[] inputBuffer = new byte[1000000];
        var inputStream = Console.OpenStandardInput(inputBuffer.Length);
        using (var reader = new StreamReader(inputStream, Console.InputEncoding, false, inputBuffer.Length)) {
            Console.SetIn(reader);
            new Solver(Console.In, Console.Out).Run();
        }
#else
        Console.SetOut(new StreamWriter(Console.OpenStandardOutput()) { AutoFlush = false });
        new Solver(Console.In, Console.Out).Run();
        Console.Out.Flush();
#endif
    }

#pragma warning disable IDE0052
    private readonly TextReader cin;
    private readonly TextWriter cout;
    private readonly TextWriter cerr;
#pragma warning restore IDE0052

    public Solver(TextReader reader, TextWriter writer)
        : base(reader) {
        cin = reader;
        cout = writer;
        cerr = Console.Error;
    }

    public Solver(string input, TextWriter writer)
        : this(new StringReader(input), writer) {
    }

#pragma warning disable IDE1006
#pragma warning disable IDE0051
    private int ni() { return NextInt(); }
    private int[] ni(int n) { return NextIntArray(n); }
    private long nl() { return NextLong(); }
    private long[] nl(int n) { return NextLongArray(n); }
    private double nd() { return NextDouble(); }
    private double[] nd(int n) { return NextDoubleArray(n); }
    private string ns() { return Next(); }
    private string[] ns(int n) { return NextArray(n); }
#pragma warning restore IDE1006
#pragma warning restore IDE0051
}

internal static class LinqPadExtension {
    [Conditional("DEBUG")]
    public static void Dump<T>(this T obj) {
#if DEBUG
        LINQPad.Extensions.Dump(obj);
#endif
    }
}

public class Scanner {
    private readonly TextReader Reader;
    private readonly CultureInfo ci = CultureInfo.InvariantCulture;

    private readonly char[] buffer = new char[2 * 1024];
    private int cursor = 0, length = 0;
    private string Token;
    private readonly StringBuilder sb = new StringBuilder(1024);

    public Scanner()
        : this(Console.In) {
    }

    public Scanner(TextReader reader) {
        Reader = reader;
    }

    public int NextInt() { return checked((int)NextLong()); }
    public long NextLong() {
        var s = Next();
        long r = 0;
        int i = 0;
        bool negative = false;
        if (s[i] == '-') {
            negative = true;
            i++;
        }
        for (; i < s.Length; i++) {
            r = r * 10 + (s[i] - '0');
#if DEBUG
            if (!char.IsDigit(s[i])) throw new FormatException();
#endif
        }
        return negative ? -r : r;
    }
    public double NextDouble() { return double.Parse(Next(), ci); }
    public string[] NextArray(int size) {
        string[] array = new string[size];
        for (int i = 0; i < size; i++) {
            array[i] = Next();
        }

        return array;
    }
    public int[] NextIntArray(int size) {
        int[] array = new int[size];
        for (int i = 0; i < size; i++) {
            array[i] = NextInt();
        }

        return array;
    }

    public long[] NextLongArray(int size) {
        long[] array = new long[size];
        for (int i = 0; i < size; i++) {
            array[i] = NextLong();
        }

        return array;
    }

    public double[] NextDoubleArray(int size) {
        double[] array = new double[size];
        for (int i = 0; i < size; i++) {
            array[i] = NextDouble();
        }

        return array;
    }

    public string Next() {
        if (Token == null) {
            if (!StockToken()) {
                throw new Exception();
            }
        }
        var token = Token;
        Token = null;
        return token;
    }

    public bool HasNext() {
        if (Token != null) {
            return true;
        }

        return StockToken();
    }

    private bool StockToken() {
        while (true) {
            sb.Clear();
            while (true) {
                if (cursor >= length) {
                    cursor = 0;
                    if ((length = Reader.Read(buffer, 0, buffer.Length)) <= 0) {
                        break;
                    }
                }
                var c = buffer[cursor++];
                if (33 <= c && c <= 126) {
                    sb.Append(c);
                } else {
                    if (sb.Length > 0) break;
                }
            }

            if (sb.Length > 0) {
                Token = sb.ToString();
                return true;
            }

            return false;
        }
    }
}
0