結果

問題 No.498 ワープクリスタル (給料日編)
ユーザー 37zigen
提出日時 2016-07-14 00:35:14
言語 Java
(openjdk 23)
結果
AC  
実行時間 309 ms / 2,000 ms
コード長 1,862 bytes
コンパイル時間 2,027 ms
コンパイル使用メモリ 77,160 KB
実行使用メモリ 77,176 KB
最終ジャッジ日時 2024-10-15 20:00:37
合計ジャッジ時間 8,041 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.util.Scanner;

public class Main {
	public static void main(String[] args) {
		solver();
	}

	static long MOD = 1_000_000_000 + 7;

	static void solver() {
		Scanner sc = new Scanner(System.in);
		int gx = sc.nextInt();
		int gy = sc.nextInt();
		int k = sc.nextInt();// the number of kinds of crystals

		P np = new P(-1, -1);
		P[] ps = new P[1 << (4 * k)];
		for (int i = 1; i < ps.length; i++) {
			ps[i] = np;
		}
		ps[0] = new P(0, 0);

		for (int i = 0; i < k; i++) {
			int x = sc.nextInt();
			int y = sc.nextInt();
			P v = new P(x, y);
			int n = sc.nextInt();

			for (int j = 0; n > 0; j++) {
				int mul = Math.min(1 << j, n);
				n -= mul;
				P u = v.multiply(mul);
				for (int l = (1 << (4 * i + j)) - 1 ; l >= 0; l--) {
					if (ps[l] == np)
						continue;
					ps[l + (mul << (4 * i))] = ps[l].add(u);
				}
			}
		}
		long ans = 0;
		for (int i = 0; i < (1 << (4 * k)); i++) {
			if (ps[i] == np)
				continue;
			if (ps[i].x != gx || ps[i].y != gy)
				continue;
			int s = i;
			long total = 0;
			long denomi = 1;
			for (; s > 0; s >>= 4) {
				int d = s & 15;
				total += d;
				denomi = denomi * inv(fact(d), MOD) % MOD;
			}
			ans = (ans + fact(total) * denomi) % MOD;
		}
		System.out.println(ans);
	}

	static class P {
		int x;
		int y;

		P(int x, int y) {
			this.x = x;
			this.y = y;
		}

		P multiply(int m) {
			return new P(x * m, y * m);
		}

		P add(P p) {
			return new P(x + p.x, y + p.y);
		}
	}

	static long fact(long nume) {
		long ans = 1;
		for (int i = 1; i <= nume; i++) {
			ans *= i;
			ans %= MOD;
		}
		return ans;
	}

	static long inv(long a, long mod) {
		a = a % mod;
		long b = mod;
		long p = 1, q = 0;

		while (b > 1) {
			long c = b / a;
			b = b % a;
			q = q - p * c;
			long d = b;
			b = a;
			a = d;
			d = p;
			p = q;
			q = d;
		}
		while (q < 0)
			q += mod;
		return q;
	}
}
0