qwb与支教 二分加容斥

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\vee B\vee C| = |A|+|B|+|C|-|A\wedge B|-|B\wedge C|-|A\wedge C|+|A\wedge B \wedge 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;
}