結果

問題 No.1973 Divisor Sequence
ユーザー tenten
提出日時 2024-04-02 14:08:50
言語 Java
(openjdk 23)
結果
MLE  
実行時間 -
コード長 2,710 bytes
コンパイル時間 2,436 ms
コンパイル使用メモリ 79,196 KB
実行使用メモリ 742,580 KB
最終ジャッジ日時 2024-09-30 23:13:18
合計ジャッジ時間 8,924 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 21 MLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #

import java.io.*;
import java.util.*;
import java.util.function.*;
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();
        List<Long> divs = new ArrayList<>();
        for (long i = 1; i * i <= m; i++) {
            if (m % i == 0) {
                divs.add(i);
                if (i * i < m) {
                    divs.add(m / i);
                }
            }
        }
        int size = divs.size();
        long[][][] matrix = new long[20][size][size];
        for (int i = 0; i < size; i++) {
            long current = m / divs.get(i);
            for (int j = 0; j < size; j++) {
                if (current % divs.get(j) == 0) {
                    matrix[0][i][j] = 1;
                }
            }
        }
        for (int i = 1; i < 20; i++) {
            for (int a = 0; a < size; a++) {
                for (int b = 0; b < size; b++) {
                    for (int c = 0; c < size; 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[size];
        ans[0] = 1;
        for (int i = 19; i >= 0; i--) {
            if (n < (1 << i)) {
                continue;
            }
            n -= (1 << i);
            long[] next = new long[size];
            for (int a = 0; a < size; a++) {
                for (int b = 0; b <size; b++) {
                    next[a] += matrix[i][a][b] * ans[b] % MOD;
                    next[a] %= MOD;
                }
            }
            ans = next;
        }
        System.out.println(ans[0]);
    }
}
class Scanner {
    BufferedReader br;
    StringTokenizer st = new StringTokenizer("");
    StringBuilder sb = new StringBuilder();
    
    public Scanner() {
        try {
            br = new BufferedReader(new InputStreamReader(System.in));
        } catch (Exception e) {
            
        }
    }
    
    public int nextInt() {
        return Integer.parseInt(next());
    }
    
    public long nextLong() {
        return Long.parseLong(next());
    }
    
    public double nextDouble() {
        return Double.parseDouble(next());
    }
    
    public String next() {
        try {
            while (!st.hasMoreTokens()) {
                st = new StringTokenizer(br.readLine());
            }
        } catch (Exception e) {
            e.printStackTrace();
        } finally {
            return st.nextToken();
        }
    }
    
}
0