#include #define REP_(i, a_, b_, a, b, ...) for (int i = (a), END_##i = (b); i < END_##i; ++i) #define REP(i, ...) REP_(i, __VA_ARGS__, __VA_ARGS__, 0, __VA_ARGS__) #define ALL(x) std::begin(x), std::end(x) using Int = long long; using Uint = unsigned long long; using Real = long double; template inline bool chmax(T &a, U b) { return a < b and ((a = b), true); } template inline bool chmin(T &a, U b) { return a > b and ((a = b), true); } template inline int ssize(const T &a) { return (int) a.size(); } template constexpr T kBigVal = std::numeric_limits::max() / 2; template inline std::ostream &print_one(const T &x, char endc) { if constexpr (std::is_same::value) { return std::cout << (x ? "Yes" : "No") << endc; } else { return std::cout << x << endc; } } template inline std::ostream &print(const T &x) { return print_one(x, '\n'); } template std::ostream &print(const T &head, Ts... tail) { return print_one(head, ' '), print(tail...); } inline std::ostream &print() { return std::cout << '\n'; } template std::ostream &print_seq(const Container &seq, const char *sep = " ", const char *ends = "\n", std::ostream &os = std::cout) { const auto itl = std::begin(seq), itr = std::end(seq); for (auto it = itl; it != itr; ++it) { if (it != itl) os << sep; os << *it; } return os << ends; } struct CastInput { template operator T() const { T x; std::cin >> x; return x; } struct Sized { int n; template operator T() const { T xs(n); for (auto &x: xs) std::cin >> x; return xs; } }; Sized operator()(int n) const { return {n}; } } in; #ifdef MY_DEBUG #include "debug_dump.hpp" #include "backward.hpp" backward::SignalHandling kSignalHandling; #else #define DUMP(...) #define cerr if(false)cerr #endif using namespace std; template inline bool has_bit(const T &x, int i) { return (x >> i) & 1; } inline int popcount(unsigned x) { return __builtin_popcount(x); } inline int popcount(Uint x) { return __builtin_popcountll(x); } // Saturating multiplication template T sat_mul(T x, T y) { static_assert(std::is_integral::value); static constexpr T kMin = std::numeric_limits::lowest(); static constexpr T kMax = std::numeric_limits::max(); static_assert(kMax != 0); if (T res; not __builtin_mul_overflow(x, y, &res)) { return res; } else if constexpr (not std::is_signed::value) { return kMax; } else { return ((x ^ y) < 0) ? kMin : kMax; } } auto solve() { int n = in; Int L = in, H = in; vector C = in(n); Int ret = 0; REP(bits, 1, 1 << n) { Int blcm = -1; bool over = false; REP(i, n) { if (has_bit(bits, i)) { if (blcm == -1) { blcm = C[i]; } else { Int g = std::gcd(blcm, C[i]); blcm = sat_mul(blcm / g, C[i]); if (blcm > H) { over = true; break; } } } } if (over) continue; int bsize = popcount((unsigned) bits); Int sign = (bsize & 1) ? 1 : -1; Int co = sign * (H / blcm - ((L - 1) / blcm)) * bsize; ret += co; } print(ret); } int main() { std::ios::sync_with_stdio(false), cin.tie(nullptr); cout << std::fixed << std::setprecision(18); const int T = 1;//in; REP(t, T) { (solve()); } }