結果

問題 No.1973 Divisor Sequence
コンテスト
ユーザー tenten
提出日時 2023-07-24 17:13:47
言語 Java
(openjdk 23)
結果
MLE  
実行時間 -
コード長 3,320 bytes
コンパイル時間 2,373 ms
コンパイル使用メモリ 85,352 KB
実行使用メモリ 738,292 KB
最終ジャッジ日時 2024-10-01 12:16:56
合計ジャッジ時間 7,996 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21 MLE * 1
権限があれば一括ダウンロードができます

ソースコード

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();
        TreeSet<Long> each = new TreeSet<>();
        for (long i = 1; i <= Math.sqrt(m); i++) {
            if (m % i == 0) {
                each.add(i);
                each.add(m / i);
            }
        }
        int idx = 0;
        int size = each.size();
        long[] values = new long[size];
        for (long x : each) {
            values[idx++] = x;
        }
        long[][][] matrix = new long[20][size][size];
        for (int i = 0; i < size; i++) {
            for (int j = 0; j < size; j++) {
                if (m / values[i] < values[j]) {
                    break;
                }
                if (each.contains(values[i] * values[j])) {
                    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];
                        matrix[i][a][c] %= MOD;
                    }
                }
            }
        }
        long[] ans = new long[size];
        Arrays.fill(ans, 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;
        }
        long result = 0;
        for (long x : ans) {
            result += x;
            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