結果

問題 No.1731 Product of Subsequence
コンテスト
ユーザー platinum
提出日時 2021-11-05 22:01:03
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,150 ms / 2,000 ms
コード長 1,691 bytes
コンパイル時間 3,252 ms
コンパイル使用メモリ 204,976 KB
最終ジャッジ日時 2025-01-25 12:18:39
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define rep(i,n) for(int i = 0; i < (int)(n); i++)

using namespace std;
using LL = long long;
using P = pair<int,int>;
using vv = vector<vector<int>>;
const int INF = (int)1e9;
const LL LINF = (LL)1e18;

long long const mod = 1000000007;
struct mint{
	long long val;
	mint(long long val = 0): val(val % mod) {}
	mint& operator += (const mint n){
		val += n.val;
		if(val >= mod) val -= mod;
		return *this;
	}
	mint& operator -= (const mint n){
		val -= n.val;
		if(val < 0) val += mod;
		return *this;
	}
	mint& operator *= (const mint n){
		val *= n.val;
		val %= mod;
		if(val < 0) val += mod;
		return *this;
	}
	mint operator + (const mint n) const{
		mint res(*this);
		return res += n;
	}
	mint operator - (const mint n) const{
		mint res(*this);
		return res -= n;
	}
	mint operator * (const mint n) const{
		mint res(*this);
		return res *= n;
	}
	mint pow(long long n) const{
		if(n == 0) return 1;
		mint m = pow(n >> 1);
		m *= m;
		if(n & 1) m *= *this;
		return m;
	}
	// mint division for prime mod
	mint inv() const{
		return pow(mod - 2);
	}
	mint& operator /= (const mint n){
		return (*this) *= n.inv();
	}
	mint operator / (const mint n) const{
		mint res(*this);
		return res /= n;
	}
};

int main(){
	int N;
	LL K;
	cin >> N >> K;
	vector<LL> A(N);
	rep(i,N) cin >> A[i];
	if(K == 1){
		mint base = 2;
		mint ans = base.pow(N);
		ans -= 1;
		cout << ans.val << endl;
		return 0;
	}
	vector<map<LL,mint>> dp(N+1);
	dp[0][K] = 1;
	rep(i,N){
		for(auto itr : dp[i]){
			LL res = itr.first;
			mint cnt = itr.second;
			dp[i+1][res] += cnt;
			LL g = gcd(res, A[i]);
			dp[i+1][res / g] += cnt;
		}
	}
	cout << dp[N][1].val << endl;

	return 0;
}
0