Олимпиадная задача - сумма чисел меньших N, которые делятся на A или на B

#include <iostream>
 
using namespace std;
 
int main(void)
{
    long long n, a, b, sum = 0;
    cin>>n>>a>>b;
 
    for (long long i = a; i < n; i += a) if (i % b != 0) sum += i;
    for (long long i = b; i < n; i += b) sum += i;
 
    cout<<(sum % 1000000007);
    return 0;
}

 int summA,int summB,int summAB,summ;
int qA=n/A;
if(n%A==0)
 qA--;
int lastA =qA*A;
summA=((A+lastA)*qA)/2;
 
int qB=n/B;
if(n%B==0)
 qB--;
int lastB =qB*B;
summB=((B+lastB)*qB)/2;
 
int qAB=n/(A*B);
if(n%(A*B)==0)
 qAB--;
int lastAB =qAB*A*B;
summAB=((A*B+lastAB)*qAB)/2;
 
summ=summA+summB-summAB;

#include "stdafx.h"
#include <iostream>
 
using namespace std;
 
int _tmain(int argc, _TCHAR* argv[])
{
 long long n, a, b, sum = 0;
    cin>>n>>a>>b;
 
    for (long long i = a; i < n; i += a) if (i % b != 0) sum += i;
    for (long long i = b; i < n; i += b) sum += i;
    cout<<endl;
    cout<<(sum % 1000000007)<<endl;
  long long summA,summB,summAB,summ;
 long long qA=n/a;
if(n%a==0)
 qA--;
 long long lastA =qA*a;
summA=((a+lastA)*qA)/2;
 
 long long qB=n/b;
if(n%b==0)
 qB--;
 long long lastB =qB*b;
summB=((b+lastB)*qB)/2;
 
 long long qAB=n/(a*b);
if(n%(a*b)==0)
 qAB--;
 long long lastAB =qAB*a*b;
summAB=((a*b+lastAB)*qAB)/2;
 
summ=summA+summB-summAB;
 
 cout<<(summ % 1000000007)<<endl;   
 return 0;
}

#include <iostream>
 
int foo(long long x, long long n, int modulo)
{
    long long first = x;            //первый член арифметической прогрессии
 
    long long last = (n / x * x);   //последний член арифметической прогресии
 
    long long count  = n / x;       //количество членов в арифметической прогрессии
 
    return ( ( ( (first % modulo + last % modulo) % modulo) * (count % modulo) ) / 2 ) % modulo;
}
 
int main()
{
    const int modulo = 1000000007;
    
    long long n;
 
    std::cin >> n;
    
    long long a, b;
    
    std::cin >> a >> b;
 
    std::cout << (foo(a, n, modulo) - foo(a * b, n, modulo) + foo(b, n, modulo)) % modulo << std::endl; 
}

#include <iostream>
 
int foo(long long x, long long n, int modulo)
{
    long long first = x;            //первый член арифметической прогрессии
 
    long long last = (n / x * x);   //последний член арифметической прогресии
 
    long long count  = n / x;       //количество членов в арифметической прогрессии
 
    return ( ( ( (first % modulo + last % modulo) % modulo) * (count % modulo) ) / 2 ) % modulo;
}
 
template <class T>
T gcd(T a, T b)
{
    while (a ^= b ^= a ^= b %= a);
    return b;
}
 
template <class T>
T lcm(T a, T b)
{
    return a * b / gcd(a, b);
}
 
int main()
{
    const int modulo = 1000000007;
    
    long long n;
 
    std::cin >> n;
    
    long long a, b;
    
    std::cin >> a >> b;
 
    std::cout << (foo(a, n, modulo) - foo(lcm(a, b), n, modulo) + foo(b, n, modulo)) % modulo << std::endl; 
}

#include "stdafx.h"
#include <iostream>
 
using namespace std;
long long NOK(long long A,long long B)
{
if(A%2==0&&B%2==0)
 return 2;
long long c=A;
if(A>B)
  c=B;
for(int i=3;i<c/2;c+=2)
{
 if(A%i==0&&B%i==0)
     return i;
}
return 1;
 
 
}
 
int _tmain(int argc, _TCHAR* argv[])
{
 long long n, a, b, sum = 0;
    cin>>n>>a>>b;
 
    for (long long i = a; i < n; i += a) if (i % b != 0) sum += i;
    for (long long i = b; i < n; i += b) sum += i;
    cout<<endl;
    cout<<(sum % 1000000007)<<endl;
  long long summA,summB,summAB,summ;
  
  
  
  long long qA=n/a;
if(n%a==0)
 qA--;
 long long lastA =qA*a;
summA=((a+lastA)*qA)/2;
 
 long long qB=n/b;
if(n%b==0)
 qB--;
 long long lastB =qB*b;
summB=((b+lastB)*qB)/2;
 
long long nok=NOK(a,b); 
long long first=a*b/nok;
 
long long qAB=n/first;
if(n%first==0)
 qAB--;
long long lastAB =qAB*first;
summAB=((first+lastAB)*qAB)/2;
 
summ=summA+summB-summAB;
 
 cout<<(summ % 1000000007)<<endl;   
 return 0;
}