ZJCOJ：qwb与支教（二分）

Problem A: qwb与支教

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 1454   Solved: 266
[ Submit][ Status][ Web Board]

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

思路：二分答案...

# include <bits/stdc++.h>
using namespace std;

typedef long long LL;
const LL INF = 0x3f3f3f3f3f3f3f3f;
LL x, y, z, n;
LL gcd(LL a, LL b)
{
    return b==0?a:gcd(b, a%b);
}

LL lcm(LL a, LL b)
{
    return a/gcd(a, b)*b;
}

bool judge(LL s)
{
    LL res = s - s/x - s/y - s/z + s/lcm(x,y) + s/lcm(x,z) + s/lcm(y,z);
    if(z/gcd(z, lcm(x,y)) <= s/lcm(x,y)) res -= s/lcm(z, lcm(x,y));
    return res >= n;
}
int main()
{
    while(~scanf("%lld%lld%lld%lld",&x,&y,&z,&n))
    {
        LL l=n, r=INF;
        while(l<r)
        {
            LL mid = l+r>>1;
            if(judge(mid)) r = mid;
            else l = mid+1;
        }
        printf("%lld\n",r);
    }
    return 0;
}