codefoeces B. Friends and Presents

本文链接：https://blog.csdn.net/wust_ZJX/article/details/40475829

You have two friends. You want to present each of them several positive integers. You want to present cnt₁ numbers to the first friend and cnt₂ numbers to the second friend. Moreover, you want all presented numbers to be distinct, that also means that no number should be presented to both friends.

In addition, the first friend does not like the numbers that are divisible without remainder by prime number x. The second one does not like the numbers that are divisible without remainder by prime number y. Of course, you're not going to present your friends numbers they don't like.

Your task is to find such minimum number v, that you can form presents using numbers from a set 1, 2, ..., v. Of course you may choose not to present some numbers at all.

A positive integer number greater than 1 is called prime if it has no positive divisors other than 1 and itself.

Input

The only line contains four positive integers cnt₁, cnt₂, x, y (1 ≤ cnt₁, cnt₂ < 10⁹; cnt₁ + cnt₂ ≤ 10⁹; 2 ≤ x < y ≤ 3·10⁴) — the numbers that are described in the statement. It is guaranteed that numbers x, y are prime.

Output

Print a single integer — the answer to the problem.

Sample test(s)

Input

3 1 2 3

Output

Input

1 3 2 3

Output

Note

In the first sample you give the set of numbers {1, 3, 5} to the first friend and the set of numbers {2} to the second friend. Note that if you give set {1, 3, 5} to the first friend, then we cannot give any of the numbers 1, 3, 5 to the second friend.

In the second sample you give the set of numbers {3} to the first friend, and the set of numbers {1, 2, 4} to the second friend. Thus, the answer to the problem is 4.

题意：求出数字v ,使得能在1——v之间，分别选择cnt1不能被x整除和 cnt2个数不能被y整除，并且cnt1和cnt2两组数中不能有相同的。

思路：二分加容斥原理

#include <iostream>
#include <stdio.h>
#include <string>
#include <cstring>
#include <algorithm>
#include <cmath>
#define ll __int64
#define N 2000000000
//容斥原理加二分。。
using namespace std;
ll cnt1,cnt2,x,y;
ll gcd(ll a,ll b)
{
    return b==0?a:gcd(b,a%b);
}

ll check(ll n)
{
    ll t1,t2,t3,temp1,temp2,temp3;

     t1=n/x;//被x整除的个数
     t2=n/y;//被y整除的个数
     t3=n/(x*y/gcd(x,y));//被x和y最小公倍数整除的个数

    if(n-t1>=cnt1&&n-t2>=cnt2&&n>=cnt1+cnt2+t3)
    return true;
    else
    return false;
}

int main()
{

    while(~scanf("%I64d%I64d%I64d%I64d",&cnt1,&cnt2,&x,&y))
    {
        ll ri=N,le=1;
        ll mid;
        ll ans=0;
        while(le<ri)
        {
            mid=(ri+le)/2;
            if(check(mid))
            {
                ans=mid;
                ri=mid;
            }
            else
            le=mid+1;
        }
        cout<<ans<<endl;
    }
    return 0;
}