結果

問題 No.890 移調の限られた旋法
コンテスト
ユーザー kakel-san
提出日時 2026-02-15 23:24:38
言語 C#
(.NET 10.0.102)
結果
AC  
実行時間 97 ms / 2,000 ms
コード長 3,863 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 17,957 ms
コンパイル使用メモリ 172,232 KB
実行使用メモリ 203,128 KB
最終ジャッジ日時 2026-02-15 23:25:00
合計ジャッジ時間 15,256 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 32
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (138 ミリ秒)。
  main -> /home/judge/data/code/bin/Release/net10.0/main.dll
  main -> /home/judge/data/code/bin/Release/net10.0/publish/

ソースコード

diff #
raw source code 

using System;
using static System.Console;
using System.Linq;
using System.Collections.Generic;

class Program
{
    static int NN => int.Parse(ReadLine());
    static int[] NList => ReadLine().Split().Select(int.Parse).ToArray();
    static string[] SList(long n) => Enumerable.Repeat(0, (int)n).Select(_ => ReadLine()).ToArray();
    public static void Main()
    {
        Solve();
    }
    static void Solve()
    {
        var c = NList;
        var (n, k) = (c[0], c[1]);
        if (n == k)
        {
            WriteLine(n == 1 ? 0 : 1);
            return;
        }
        var dlist = DivList(n);
        var ans = 0L;
        var mod = 1_000_000_007;
        var ncr = new NCR(n, mod);
        var mdic = new Dictionary<int, long>();
        foreach (var d in dlist)
        {
            // 長さN/dのパターン
            if (k % d != 0) continue;
            var mul = 1L;
            if (mdic.ContainsKey(n / d)) mul = (1 - mdic[n / d] + mod) % mod;
            ans = (ans + ncr.Calc(n / d, k / d) * mul % mod) % mod;
            var ddlist = DivList(n / d);
            foreach (var dd in ddlist)
            {
                if (k / d % dd != 0) continue;
                // 長さN/dのパターンを数えたときに含まれる長さN/d/ddのパターン数
                var len = n / d / dd;
                if (mdic.ContainsKey(len)) mdic[len] = (mdic[len] + mul) % mod;
                else mdic[len] = mul;
            }
        }
        WriteLine(ans);
    }
    static List<int> DivList(int n)
    {
        var list = new List<int>();
        var rev = new List<int>();
        for (var i = 2; i * i <= n; ++i)
        {
            if (n % i == 0)
            {
                list.Add(i);
                if (i * i < n) rev.Add(n / i);
            }
        }
        rev.Reverse();
        list.AddRange(rev);
        return list;
    }
    class NCR
    {
        int[] facts;
        int[] revFacts;
        int mod;
        public NCR(int n, int mod)
        {
            facts = new int[n + 1];
            revFacts = new int[n + 1];
            this.mod = mod;
            facts[0] = 1;
            var tmp = 1L;
            for (var i = 1; i <= n; ++i)
            {
                tmp = tmp * i % mod;
                facts[i] = (int)tmp;
            }
            tmp = Exp(facts[n], mod - 2);
            revFacts[n] = (int)tmp;
            for (var i = n; i > 1; --i)
            {
                tmp = tmp * i % mod;
                revFacts[i - 1] = (int)tmp;
            }
            revFacts[0] = 1;
        }
        public long Exp(long n, long k)
        {
            n = n % mod;
            if (k == 0) return 1;
            if (k == 1) return n;
            var half = Exp(n, k / 2);
            var result = half * half % mod;
            return ((k % 2) == 0) ? result : (result * n % mod);
        }
        public long Calc(int n, int r)
        {
            return (long)facts[n] * revFacts[r] % mod * revFacts[n - r] % mod;
        }
        /// <summary>nが大きくrが小さい場合の計算</summary>
        public long Calc2(int n, int r)
        {
            var tmp = 1L;
            for (var i = 0; i < r; ++i)
            {
                tmp = tmp * (n - i) % mod;
            }
            return tmp * revFacts[r] % mod;
        }
        public long NPR(int n, int r)
        {
            return (long)facts[n] * revFacts[n - r] % mod;
        }
        public long Fact(int n)
        {
            return facts[n];
        }
        public long RevFact(int n)
        {
            return revFacts[n];
        }
        public long Inverse(int n)
        {
            return (long)revFacts[n] * facts[n - 1] % mod;
        }
    }
}
0