Lucky Number ZOJ - 3233 容斥原理 +两个条件

这个地方其实也不难。。不过还是看了别人的博客
http://blog.csdn.net/acm_cxlove/article/details/7878893
http://blog.csdn.net/acm_baihuzi/article/details/49304911
而且自己的代码不知道哪里错了。可能什么地方溢出了
只能借别人的,,可能就是lcm这个地方我的溢出了?

ll gcd(ll a,ll b) {
    return b==0?a:gcd(b,a%b);
}

ll lcm(ll a,ll b) {
    if(a==-1||b==-1) return -1;
    ll gd=gcd(a,b);
    if(inf/b<a/gd)   return -1;
    return a/gd*b;
}
void solve() {
    ll sum=0;
    for(int i=1; i<(1<<n); i++) {
        ll lc=1;
        int cnt=0;
        ll temp=i;
        for(int j=0; j<n; j++) {
            if(temp&1)cnt++,lc=lcm(lc,a[j]);
            temp>>=1;
        }
        if(lc==-1)lc=inf;
        ll tmp=lcm(bm,lc);
        if(tmp==-1) tmp=inf;
        ll num=(r/lc-r/tmp)-(l/lc-l/tmp);
        if(cnt&1)sum+=num;
        else     sum-=num;
    }
    printf("%lld\n",sum);
}
int main() {
    //freopen("test.in","r",stdin);
    while(~scanf("%d%d%lld%lld",&n,&m,&l,&r)&&n) {
        l--;
        for(int i=0; i<n; i++)scanf("%d",&a[i]);
        bm=1;
        for(int i=0; i<m; i++) {
            ll x;
            scanf("%lld",&x);
            bm=lcm(bm,x);
        }
        solve();
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值