結果

問題 No.125 悪の花弁
コンテスト
ユーザー 古寺いろは
提出日時 2015-04-06 03:13:05
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,312 bytes
コンパイル時間 1,443 ms
コンパイル使用メモリ 165,652 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-04 03:17:24
合計ジャッジ時間 2,289 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2 WA * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;

long long mod = (long long)1e9 + 7;

int gcd(int a, int b){
	if (b == 0) return a;
	return gcd(b, a%b);
}

long long powmod(long long a, int p){
	if (p == 0) return 1;
	if (p % 2 == 1){
		return (powmod(a, p - 1) * a) % mod;
	}
	long long mid = powmod(a, p / 2);
	return (mid * mid) % mod;
}

long rev(long long a){
	return powmod(a, mod - 2);
}

int main() {
	int K;
	cin >> K;
	vector<int> C(K);
	int sum = 0;
	for (int i = 0; i < K; i++)
	{
		cin >> C[i];
		sum += C[i];
	}
	int g = C[0];
	for (int i = 0; i < K; i++)
	{
		g = gcd(g, C[i]);
	}

	long long ans = 0;

	vector<int> factor;
	vector<long long> num;

	for (int i = g; i >= 1; i--)
	{
		if (g%i != 0) continue;
		int remain = sum / i - 1;

		long long mul = 1;
		long long div = 1;
		for (int j = 0; j < K; j++)
		{
			int start = 0;
			if (j == 0) start = 1;
			for (int l = start; l < C[j] / i; l++)
			{
				mul *= remain--;
				mul %= mod;
				div *= l + 1 - start;
				div %= mod;
			}
		}
		mul *= rev(div);
		for (int j = 0; j < factor.size(); j++)
		{
			if (factor[j] % i != 0) continue;
			mul -= num[j];
		}
		mul %= mod;
		mul += mod;
		mul %= mod;


		mul *= rev(C[0] / i);
		mul %= mod;
		factor.push_back(i);
		num.push_back(mul);
		ans += mul;
		ans %= mod;

	}

	cout << ans << endl;
}
0