PIPIOJ1341: PIPI学容斥

#include<bits/stdc++.h> 
using namespace std;
const double PI=acos(-1.0);
#define _max -0x7fffffff
long long int  gcd(long long int a,long long int b) //求a,b的最大公约数 
{
	if(a<b)
	{
	int temp=a;
	a=b;
	b=temp;
	}
	if(a%b==0)
		return b;
	else
		return gcd(b,a%b);	
} 
long long int ab(long long int a,long long int b)  //求a,b的最小公倍数 
{
	return a*b/gcd(a,b);
}
int main()
{
	long long int n,a,b,c;
	while(scanf("%lld%lld%lld%lld",&n,&a,&b,&c)!=EOF)
	{
		long long int temp=ab(a,b);
		printf("%lld\n",n/a+n/b+n/c-n/ab(a,b)-n/ab(a,c)-n/ab(b,c)+n/ab(temp,c));	
	}
	return 0;	
}