Cod sursa(job #2568953)

Utilizator RazvanPanaiteRazvan Panaite RazvanPanaite Data 4 martie 2020 10:35:52
Problema Suma si numarul divizorilor Scor 70
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.93 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 ans,diviz;

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*=(explog(prime[i],cnt+1)+MOD-1)%MOD;
                ans%=MOD;
                ans*=invmod(prime[i]-1);
                ans%=MOD;
            }
        if(n > 1){
            diviz*=2;
            ans=(ans*(explog(n,2)-1))%MOD;
            ans=ans*invmod(n-1)%MOD;
        }
        fout<<diviz<<' '<<ans<<'\n';
    }

    return 0;
}
void ciur(){
    ll i,j;
    for(i=2;i<=1e6;i++)
        if(!uz[i]){
            prime.pb(i);
            for(j=i*i;j<=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-2);
}