結果
| 問題 | No.546 オンリー・ワン | 
| コンテスト | |
| ユーザー |  Manuel1024 | 
| 提出日時 | 2022-03-30 00:06:47 | 
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 3 ms / 2,000 ms | 
| コード長 | 1,041 bytes | 
| コンパイル時間 | 655 ms | 
| コンパイル使用メモリ | 69,504 KB | 
| 実行使用メモリ | 5,248 KB | 
| 最終ジャッジ日時 | 2024-11-08 23:35:18 | 
| 合計ジャッジ時間 | 1,438 ms | 
| ジャッジサーバーID (参考情報) | judge5 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 4 | 
| other | AC * 7 | 
ソースコード
#include <iostream>
#include <vector>
using namespace std;
using ll = long long;
ll mygcd(ll a, ll b){
    while(b != 0){
        ll r = a%b;
        a = b; b = r;
    }
    return a;
}
int main(){
    ll n, l, h; cin >> n >> l >> h;
    vector<ll> c(n);
    for(auto &it: c) cin >> it;
    auto calc = [&](ll x){
        ll res = 0;
        for(int bit = 1; bit < (1 << n); bit++){
            ll m = 1;
            int cnt = 0;
            bool iszero = false;
            for(int i = 0; i < n; i++){
                if((bit >> i) & 1){
                    m *= c[i]/mygcd(m, c[i]);
                    cnt++;
                    if(m > x){
                        iszero = true;
                        break;
                    }
                }
            }
            if(!iszero){
                if(cnt%2 == 0) res -= (x/m)*cnt;
                else res += (x/m)*cnt;
            }
        }
        return res;
    };
    cout << calc(h)-calc(l-1) << endl;
    // cerr << calc(h) << " " << calc(l-1) << endl;
    return 0;
}
            
            
            
        