結果

問題 No.546 オンリー・ワン
ユーザー GOTKAKO
提出日時 2025-07-22 10:13:34
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 940 bytes
コンパイル時間 2,226 ms
コンパイル使用メモリ 198,848 KB
実行使用メモリ 7,716 KB
最終ジャッジ日時 2025-07-22 10:13:38
合計ジャッジ時間 3,449 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 7
権限があれば一括ダウンロードができます

ソースコード

diff #

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

template<typename T>
void FastMobius(vector<T> &A){
    int n = A.size();
    for(int i=1; i<n; i<<=1){
        for(int k=0; k<n; k++){
            if(i&k) A.at(k) -= A.at(k-i);
            else k += i-1;
        }
    }
}

int main(){
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    long long N,L,R; cin >> N >> L >> R;
    vector<long long> C(N);
    for(auto &c : C) cin >> c;

    int n2 = 1<<N;
    vector<long long> A(n2);
    for(int i=0; i<N; i++) A.at(1<<i) = 1;
    FastMobius(A);

    vector<long long> LCM(n2,1);
    for(int i=1; i<n2; i++) for(int k=0; k<N; k++) if(i&(1<<k)){
        long long lcm = LCM.at(i-(1<<k));
        LCM.at(i) = lcm*C.at(k)/gcd(lcm,C.at(k));
    }

    long long answer = 0,inf = 1001001001;
    for(int i=1; i<n2; i++){
        long long lcm = LCM.at(i);
        answer += (R/lcm-(L-1)/lcm)*A.at(i);
    }
    cout << answer << endl;
}
0