Pagini recente » Cod sursa (job #373039) | Cod sursa (job #3285241) | Cod sursa (job #336412) | Cod sursa (job #3122795) | Cod sursa (job #3250385)
#include <fstream>
#define NMAX 1000001
using namespace std;
ifstream f("mins.in");
ofstream g("mins.out");
long long c,d,mn,mx,rez,i,j,k,fact,p,x,m,s,nrc;
long long v[NMAX];
int main()
{
f>>c>>d;
c--; d--;
mn=min(c,d); mx=max(c,d);
rez=mx; // toate numerele sunt prime cu 1
for(k=2;k<=mn;k++)
{
// se construieste vectorul cu factorii primi ai lui k (v cu m elem)
x=k;
m=0;
fact=2;
while (x>1)
{
if(x%fact==0) v[++m]=fact;
while(x%fact==0) x=x/fact;
i++;
if(fact*fact>x) fact=x;
}
s=mx; // s=nr. de numere prime cu i care sunt <= mx (scad din mx numerele neprime cu i)
for(i=1; i < (1LL << m); i++)
{
nrc = 0; // numarul de elemente din submultimea curenta
p = 1; // produsul elementelor din submultimea curenta
for(j = 0;j < m; j++)
if(i & (1LL << j))
{
nrc++;
p = p * v[j + 1];
}
if(nrc % 2 == 1) s = s - mx / p;
else s = s + mx / p;
}
rez+=s;
}
g<<rez<<'\n';
return 0;
}