Description
qwb同时也是是之江学院的志愿者,暑期要前往周边地区支教,为了提高小学生的数学水平。她把小学生排成一排,从左至右从1开始依次往上报数。
玩完一轮后,他发现这个游戏太简单了。于是他选了3个不同的数x,y,z;从1依次往上开始报数,遇到x的倍数、y的倍数或z的倍数就跳过。如果x=2,y=3,z=5;第一名小学生报1,第2名得跳过2、3、4、5、6,报7;第3名得跳过8、9、10,报11。
那么问题来了,请你来计算,第N名学生报的数字是多少?
Input
多组测试数据,处理到文件结束。(测试数据数量<=8000)
每个测试例一行,每行有四个整数x,y,z,N。( 2≤x,y,z≤107,1≤N≤1017)。
Output
对于每个测试例,输出第N名学生所报的数字,每个报数占一行。
Sample Input
2 3 5 2
6 2 4 10000
Sample Output
7
19999
Hint
题意
题解:
由于n很大 暴力是不行的 于是考虑二分写法 思路就是找一个很大的数减去这个区间内所有x,y,z的倍数的个数 就是这个很大的数是第几个报数的人啦
x,y,z的倍数个数要用要容斥定理
|A∨B∨C|=|A|+|B|+|C|−|A∧B|−|B∧C|−|A∧C|+|A∧B∧C|
因为是倍数 重复的就是最小公倍数这些啦(这么说有点不严谨- -
AC代码
#include <cstdio>
typedef long long ll;
ll x,y,z,n;
ll gcd(ll a,ll b){
return !b?a:gcd(b,a%b);
}
ll lcm(ll a,ll b){
return (a/gcd(a,b))*b;
}
ll check(ll mid){
ll res = mid - mid/x-mid/y-mid/z+mid/lcm(x,y)+mid/lcm(x,z)+mid/lcm(y,z)-mid/lcm(lcm(x,y),z);
if (res>=n) return 1;
return 0;
}
int main(){
while (~scanf("%lld%lld%lld%lld",&x,&y,&z,&n)){
ll l = 1,r = 1e18;
ll ans,mid;
while (l<=r){
mid = (l+r)>>1;
if (check(mid)){
r = mid-1;
ans = mid;
}else l = mid+1;
}
printf("%lld\n",ans);
}
return 0;
}