結果

問題 No.1973 Divisor Sequence
ユーザー tenten
提出日時 2023-07-24 17:34:52
言語 Java
(openjdk 23)
結果
AC  
実行時間 86 ms / 2,000 ms
コード長 3,388 bytes
コンパイル時間 2,271 ms
コンパイル使用メモリ 90,124 KB
実行使用メモリ 39,324 KB
最終ジャッジ日時 2024-10-01 12:42:46
合計ジャッジ時間 4,313 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.io.*;
import java.util.*;
import java.util.stream.*;

public class Main {
    static final int MOD = 1000000007;
    public static void main(String[] args) throws Exception {
        Scanner sc = new Scanner();
        int n = sc.nextInt() - 1;
        long m = sc.nextLong();
        HashMap<Long, Integer> counts = new HashMap<>();
        for (long i = 2; i <= Math.sqrt(m); i++) {
            while (m % i == 0) {
                counts.put(i, counts.getOrDefault(i, 0) + 1);
                m /= i;
            }
        }
        if (m > 1) {
            counts.put(m, 1);
        }
        long result = 1;
        for (int x : counts.values()) {
            long[][][] matrix = new long[20][x + 1][x + 1];
            for (int i = 0; i <= x; i++) {
                for (int j = 0; j + i <= x; j++) {
                    matrix[0][i][j] = 1;
                }
            }
            for (int i = 1; i < 20; i++) {
                for (int a = 0; a <= x; a++) {
                    for (int b = 0; b <= x; b++) {
                        for (int c = 0; c <= x; c++) {
                            matrix[i][a][c] += matrix[i - 1][a][b] * matrix[i - 1][b][c] % MOD;
                            matrix[i][a][c] %= MOD;
                        }
                    }
                }
            }
            long[] ans = new long[x + 1];
            int current = n;
            Arrays.fill(ans, 1);
            for (int i = 19; i >= 0; i--) {
                if (current < (1 << i)) {
                    continue;
                }
                current -= (1 << i);
                long[] next = new long[x + 1];
                for (int a = 0; a <= x; a++) {
                    for (int b = 0; b <= x; b++) {
                        next[a] += matrix[i][a][b] * ans[b] % MOD;
                        next[a] %= MOD;
                    }
                }
                ans = next;
            }
            long tmp = 0;
            for (long y : ans) {
                tmp += y;
                tmp %= MOD;
            }
            result *= tmp;
            result %= MOD;
        }
        System.out.println(result);
    }
    
}
class Utilities {
    static String arrayToLineString(Object[] arr) {
        return Arrays.stream(arr).map(x -> x.toString()).collect(Collectors.joining("\n"));
    }
    
    static String arrayToLineString(int[] arr) {
        return String.join("\n", Arrays.stream(arr).mapToObj(String::valueOf).toArray(String[]::new));
    }
}
class Scanner {
    BufferedReader br = new BufferedReader(new InputStreamReader(System.in));
    StringTokenizer st = new StringTokenizer("");
    StringBuilder sb = new StringBuilder();
    
    public Scanner() throws Exception {
        
    }
    
    public int nextInt() throws Exception {
        return Integer.parseInt(next());
    }
    
    public long nextLong() throws Exception {
        return Long.parseLong(next());
    }
    
    public double nextDouble() throws Exception {
        return Double.parseDouble(next());
    }
    
    public int[] nextIntArray() throws Exception {
        return Stream.of(br.readLine().split(" ")).mapToInt(Integer::parseInt).toArray();
    }
    
    public String next() throws Exception {
        while (!st.hasMoreTokens()) {
            st = new StringTokenizer(br.readLine());
        }
        return st.nextToken();
    }
    
}
0