2017年第0届浙江工业大学之江学院程序设计竞赛决赛 Ａ．qwb与支教【二分＋简单容斥】

Problem A: qwb与支教

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 1361   Solved: 225
[ 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

思路：

我们设定ans是第n个数的话，那么对应有：

n=ans/x+ans/y+ans/z-ans/lcm(x,y)-ans/lcm(x,z)-ans/lcm(y,z)+ans/lcm(x,y,z)；

而且这里包含单调性问题，我们直接二分答案去check即可。

Ac代码：

#include<stdio.h>
#include<string.h>
using namespace std;
#define ll long long int
ll gcd(ll x,ll y)
{
    return y==0?x:gcd(y,x%y);
}
ll x,y,z,n,xy,xz,yz,xyz;
ll Slove(ll ans)
{
    ll tmp=ans/x+ans/y+ans/z;
    tmp-=ans/xy;
    tmp-=ans/xz;
    tmp-=ans/yz;
    tmp+=ans/xyz;
    return tmp;
}
int main()
{
    while(~scanf("%lld%lld%lld%lld",&x,&y,&z,&n))
    {
        xy=x*y/gcd(x,y);
        xz=x*z/gcd(x,z);
        yz=y*z/gcd(y,z);
        xyz=xy*z/gcd(xy,z);
        ll l=0;
        ll r=2000000000000000000;
        ll ans=-1;
        while(r-l>=0)
        {
            ll mid=(l+r)/2;
            //printf("%lld %lld\n",mid,mid-Slove(mid));
            if(mid-Slove(mid)>n)
            {
                r=mid-1;
            }
            else if(mid-Slove(mid)<n)
            {
                l=mid+1;
            }
            else
            {
                ans=mid;
                r=mid-1;
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}