Lucky Number ZOJ - 3233 （容斥）

Lucky Number

 ZOJ - 3233 

题意：给出数组a和b，问区间[L，R]之间有多少个数至少能被数组a中的一个数整除且至少不能被数组b中的一个数整除。

对于数组a，直接容斥

数组b，可以求全都能被乘除的情况，减去就可以了。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
ll a[20],b;
ll n,m,L,R;

ll gcd(ll a,ll b){
    return b?gcd(b,a%b):a;
}
ll lcm(ll a,ll b){
    if(a/gcd(a,b)>R/b) return R+1;
    return a/gcd(a,b)*b;
}
int main(){
    while(scanf("%lld%lld%lld%lld",&n,&m,&L,&R)&&(n||m||L||R)){
        ll q=1;
        for(ll i=0;i<n;i++) scanf("%lld",&a[i]);
        for(ll i=0;i<m;i++) {
            scanf("%lld",&b);
            q=lcm(q,b);
        }
        ll sta=1<<n;
        ll cnt=0,ans=0;
        for(ll i=1;i<sta;i++){
            cnt=0;
            ll temp=1;
            for(ll j=0;j<n&&temp<=R;j++) if((i>>j)&1) cnt++,temp=lcm(temp,a[j]);
            if(temp>R) continue;
            ll s=lcm(q,temp);
            if(cnt&1) ans+=(R/temp-(L-1)/temp)-(R/s-(L-1)/s);
            else ans-=(R/temp-(L-1)/temp)-(R/s-(L-1)/s);
        }
        printf("%lld\n",ans);
    }
}