結果

問題 No.665 Bernoulli Bernoulli
ユーザー GrenacheGrenache
提出日時 2018-03-12 01:22:12
言語 Java21
(openjdk 21)
結果
AC  
実行時間 407 ms / 2,000 ms
コード長 3,057 bytes
コンパイル時間 4,999 ms
コンパイル使用メモリ 77,740 KB
実行使用メモリ 54,556 KB
最終ジャッジ日時 2024-04-23 15:20:44
合計ジャッジ時間 12,549 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 133 ms
53,984 KB
testcase_01 AC 134 ms
54,104 KB
testcase_02 AC 381 ms
54,220 KB
testcase_03 AC 405 ms
54,344 KB
testcase_04 AC 390 ms
53,984 KB
testcase_05 AC 375 ms
54,324 KB
testcase_06 AC 355 ms
54,108 KB
testcase_07 AC 363 ms
54,244 KB
testcase_08 AC 364 ms
54,296 KB
testcase_09 AC 395 ms
54,368 KB
testcase_10 AC 345 ms
54,276 KB
testcase_11 AC 404 ms
54,164 KB
testcase_12 AC 393 ms
54,004 KB
testcase_13 AC 379 ms
54,156 KB
testcase_14 AC 407 ms
54,128 KB
testcase_15 AC 363 ms
54,420 KB
testcase_16 AC 381 ms
54,556 KB
testcase_17 AC 355 ms
53,996 KB
testcase_18 AC 361 ms
54,316 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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


public class Main_yukicoder665 {

	private static Scanner sc;
	private static Printer pr;

	private static void solve() {
		final int MOD = 1_000_000_007;

		long n = sc.nextLong();
		int nn = (int)(n % MOD);
		int k = sc.nextInt();

//		long stime = System.currentTimeMillis();

		PC pc = new PC(k + 2, MOD);
		long[] inv = new long[k + 2];
		inv[0] = 1;
		for (int i = 1; i <= k + 1; i++) {
			inv[i] = pc.fact[i - 1] * pc.ifact[i] % MOD;
		}

//		pr.println(System.currentTimeMillis() - stime);

		long[] b = new long[k + 1];
		b[0] = 1;
		for (int i = 1; i <= k; i++) {
			long tmp = 0;
//			long C = 1;
			if (i >= 3 && i % 2 == 1) {
				continue;
			}
			for (int j = 0; j < i; j++) {
				if (j >= 3 && j % 2 == 1) {
					continue;
				}
//				tmp += b[j] * pc.C(i + 1, j) % MOD;
	    		long C = pc.fact[i + 1] * pc.ifact[i + 1 - j] % MOD * pc.ifact[j] % MOD;
				tmp += b[j] * C % MOD;
				if (tmp >= MOD) {
					tmp -= MOD;
				}
//				C *= i + 1 - j;
//				C %= MOD;
//				C *= inv[j + 1];
//				C *= MOD;
//				tmp %= MOD;
			}
			b[i] = (-tmp + MOD) % MOD * inv[i + 1] % MOD;
		}

//		pr.println(System.currentTimeMillis() - stime);

		long tmp = 0;
		for (int i = 0; i <= k; i++) {
			tmp += (long)pc.C(k + 1, i) * pc.pow(nn + 1, k + 1 - i) % MOD * b[i] % MOD;
			tmp %= MOD;
		}
		tmp *= inv[k + 1] % MOD;
		tmp %= MOD;

		pr.println(tmp);

//		pr.println(System.currentTimeMillis() - stime);
	}

    private static class PC {
    	// MOD must be prime number.
    	int MOD;
    	// fact[i] : i! % MOD
    	long[] fact;
    	// ifact[i] : 1/i! % MOD
    	long[] ifact;

    	PC(int size, int MOD) {
    		// O(size)
    		// n=sizeまでのnCrを求める。
    		// nHrはn+r-1Crになってしまうので注意

    		this.MOD = MOD;

    		fact = new long[size + 1];
    		fact[0] = 1;
    		for (int i = 1; i <= size; i++) {
    			fact[i] = fact[i - 1] * i % MOD;
    		}

    		ifact = new long[size + 1];

    		int loop = MOD - 2;
    		long x = fact[size];
    		ifact[size] = 1;
    		while (loop > 0) {
    			if (loop % 2 == 1) {
    				ifact[size] = ifact[size] * x % MOD;
    			}
    			x = x * x % MOD;
    			loop /= 2;
    		}

    		for (int i = size - 1; i >= 0; i--) {
    			ifact[i] = ifact[i + 1] * (i + 1) % MOD;
    		}

    	}

    	// 組合せの数
    	long C(int n, int r) {
//    		if (r > n) {
//    			return 0;
//    		}

    		return fact[n] * ifact[n - r] % MOD * ifact[r] % MOD;
    	}

    	long pow(int a, long n) {
    		long loop = n;
    		long ret = 1;
    		long x = a;
    		while (loop > 0) {
    			if (loop % 2 == 1) {
    				ret = ret * x % MOD;
    			}
    			x = x * x % MOD;
    			loop /= 2;
    		}

    		return ret;
    	}
    }

	// ---------------------------------------------------
	public static void main(String[] args) {
		sc = new Scanner(System.in);
		pr = new Printer(System.out);

		solve();

		pr.close();
		sc.close();
	}

	private static class Printer extends PrintWriter {
		Printer(PrintStream out) {
			super(out);
		}
	}
}
0