Cod sursa(job #2569950)

Utilizator RazvanPanaiteRazvan Panaite RazvanPanaite Data 4 martie 2020 14:26:01
Problema Suma si numarul divizorilor Scor 70
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.95 kb
#include <bits/stdc++.h>
#define pb push_back

using namespace std;

ifstream fin("ssnd.in");
ofstream fout("ssnd.out");

void debug_out() { cerr << '\n'; }
template <typename Head, typename... Tail> void debug_out(Head H, Tail... T) { cerr << " " << H; debug_out(T...);}
#define dbg(...) cerr << #__VA_ARGS__ << " ->", debug_out(__VA_ARGS__)
#define dbg_v(x, n) do{cerr<<#x"[]: ";for(int _=0;_<n;++_)cerr<<x[_]<<" ";cerr<<'\n';}while(0)
#define dbg_ok cerr<<"OK!\n"

typedef pair<int,int> pii;
typedef long long int ll;
typedef long double ld;

const int DMAX = 1e6+10;
const int MOD = 9973;

vector <int> prime;

bool uz[DMAX];

ll n;
ll diviz;
int ans;

void ciur();
int explog(int x,int n);
int invmod(int x);

int main(){
    ios_base::sync_with_stdio(false);
    cin.tie(0);

    int t,i,j;
    int cnt;

    ciur();

    fin>>t;
    while(t--){
        fin>>n;
        diviz=ans=1;
        for(i=0;i<prime.size() && 1LL*prime[i]*prime[i] <= n;i++){
            if(n%prime[i] % 2 == 0){
                cnt=0;
                while(n%prime[i] == 0){
                    cnt++;
                    n/=prime[i];
                }
                diviz*=(cnt+1);
                ans=ans*((explog(prime[i]%MOD,cnt+1)+MOD-1)%MOD)%MOD;
                ans=(ans*invmod(prime[i]-1))%MOD;
            }
        }
        if(n > 1){
            diviz*=2;
            ans=(ans*(n+1))%MOD;
        }
        fout<<diviz<<' '<<ans<<'\n';
    }

    return 0;
}
void ciur(){
    int i,j;
    for(i=2;i<=(int)1e6;i++)
        if(!uz[i]){
            prime.pb(i);
            if(1LL*i*i > (int)1e6)
                continue;
            for(j=i*i;j<=(int)1e6;j+=i)
                uz[j]=true;
        }
}
int explog(int x,int n){
    if(!n)
        return 1;
    if(n & 1)
        return 1LL*x*explog(1LL*x*x%MOD,n/2)%MOD;
    return explog(1LL*x*x%MOD,n/2);
}
int invmod(int x){
    return explog(x%MOD,MOD-2);
}