結果

問題 No.546 オンリー・ワン
コンテスト
ユーザー beet
提出日時 2018-11-29 17:55:08
言語 C++17
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 987 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,652 ms
コンパイル使用メモリ 196,968 KB
最終ジャッジ日時 2025-01-06 17:46:02
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 7
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
using namespace std;
using Int = long long;
template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;}
template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;}

//INSERT ABOVE HERE
signed main(){
  Int n,l,h;
  cin>>n>>l>>h;
  vector<Int> c(n);
  for(Int i=0;i<n;i++) cin>>c[i];

  const Int MAX = 1e9+100;
  auto lcm=[&](Int x,Int y){
             if(max(x,y)>=MAX) return MAX;
             Int l=x/__gcd(x,y)*y;
             return min(MAX,l);
           };

  auto calc=
    [&](Int x){
      Int s=1<<n;
      vector<Int> dp(s);
      for(Int b=0;b<s;b++){
        Int l=1;
        for(Int i=0;i<n;i++)
          if((b>>i)&1) l=lcm(l,c[i]);
        dp[b]=x/l;
      }

      for(Int i=0;i<n;i++)
        for(Int b=0;b<s;b++)
          if((~b>>i)&1) dp[b]-=dp[b|(1<<i)];
      
      Int res=0;
      for(Int i=0;i<n;i++)
        res+=dp[1<<i];
      return res;
    };
  
  cout<<calc(h)-calc(l-1)<<endl;
  return 0;
}
0