这个地方其实也不难。。不过还是看了别人的博客
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;
}