Sunday, October 22, 2017

10392 - Factoring Large Numbers

/***

            Bismillahir Rahmanir Rahim
            Read the name of Allah, who created you!!!
            Author : Shah Newaj Rabbi Shishir
            Department of CSE, City University, Bangladesh.

***/

#include <bits/stdc++.h>
using namespace std;

#define sf scanf
#define pf printf
#define scase sf ("%d",&tc)
#define sn sf ("%d",&n)
#define whilecase while (tc--)
#define eof while (cin >> n)
#define forloop for (pos=1; pos<=tc; pos++)
#define arrayloop (i=0; i<n; i++)
#define cinstr cin >> str
#define getstr getline (cin,str)
#define pcase pf ("Case %d: ",pos)
#define pb push_back
#define in insert
#define llu unsigned long long
#define lld long long
#define u unsigned int

const int MAX = 1000005;
vector <int> v;
bool prime[MAX];

void sieve ()
{
    int i,j;

    v.pb(2);

    prime[0] = prime[1] = true;

    for (i=4; i<=MAX; i+=2)
        prime[i] = true;

    for (i=3; i*i<=MAX; i+=2)
    {
        if (!prime[i])
        {
            v.pb(i);

            for (j=i*i; j<=MAX; j+=2*i)
                prime[j] = true;
        }
    }

    for (i=1001; i<=MAX; i+=2)
        if (!prime[i])
            v.pb(i);

    return;
}

int main (void)
{
    sieve ();

    lld n,i,len=v.size();

    while (sf ("%lld",&n) && n != -1)
    {
        for (i=0; i<len; i++)
        {
            if (n % v[i] == 0)
            {
                while (n % v[i] == 0)
                {
                    pf ("    %d\n",v[i]);
                    n /= v[i];
                }
            }
        }

        if (n > 1)
            pf ("    %lld\n",n);

        pf ("\n");
    }

    return 0;
}

No comments: