Cod sursa(job #2854653)

Utilizator LionMan101Achim-Panescu Silvian LionMan101 Data 21 februarie 2022 16:38:05
Problema Suma si numarul divizorilor Scor 70
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 2.94 kb
#include <bits/stdc++.h>

using namespace std;

using ll = long long;
const int INF=1e9+7;
const ll INFF=1e18+7;
#define nl '\n'
void __print(int x) { cout << x; }
void __print(long x) { cout << x; }
void __print(long long x) { cout << x; }
void __print(unsigned x) { cout << x; }
void __print(unsigned long x) { cout << x; }
void __print(unsigned long long x) { cout << x; }
void __print(float x) { cout << x; }
void __print(double x) { cout << x; }
void __print(long double x) { cout << x; }
void __print(char x) { cout << '\'' << x << '\''; }
void __print(const char* x) { cout << '\"' << x << '\"'; }
void __print(const string& x) { cout << '\"' << x << '\"'; }
void __print(bool x) { cout << (x ? "true" : "false"); }
template<typename T, typename V>
void __print(const pair<T, V>& x) { cout << '{'; __print(x.first); cout << ','; __print(x.second); cout << '}'; }
template<typename T>
void __print(const T& x) { int f = 0; cout << '{'; for (auto& i : x) cout << (f++ ? "," : ""), __print(i); cout << "}"; }
void _print() { cout << "]\n"; }
template <typename T, typename... V>
void _print(T t, V... v) { __print(t); if (sizeof...(v)) cout << ", "; _print(v...); }
#ifndef ONLINE_JUDGE
#define debug(x...) cout << "[" << #x << "] = ["; _print(x)
#else
#define debug(x...)
#endif

const ll mod=9973;

ll power(ll x, ll p){
    if(p==0) return 1;
    if(p&1){
        ll nr=power(x,p/2);
        return (((nr*nr)%mod)*x)%mod;
    }
    else{
        ll nr=power(x,p/2);
        return (nr*nr)%mod;
    }
}

void solve()
{
    int n; cin>>n;
    ll maxx=0;
    vector<ll> a(n);
    for(int i=0; i<n; i++){
        cin>>a[i];
        maxx=max(maxx,a[i]);
    }
    ll m=1e6+5;
    vector<bool> ciur(m);
    vector<ll> prim;  
    for(ll i=2; i*i<=maxx; i++){
        if(ciur[i]) continue;
        prim.emplace_back(i);
        ll p=i+i;
        while(p<=m){
            ciur[p]=true;
            p+=i;
        }
    }
    
    for(int i=0; i<n; i++){
        int it=0;
        ll sum=1, div=1;
        while(prim[it]*prim[it]<=a[i]){
            ll nr=0;
            while(a[i]%prim[it]==0 && a[i]>1){
                a[i]/=prim[it];
                nr++;
            }
            ll p1 = (power(prim[it], nr+1) - 1) % mod;
            ll p2 = power(prim[it]-1, mod-2) % mod;
            sum = (((sum * p1) % mod) * p2) % mod;
            div*=nr+1;
            it++;
        }
        if(a[i]>1){
            sum = (sum*(a[i]+1)) % mod;
            div*=2;
        }
        sum%=mod;
        cout<<div<<" "<<sum<<nl;
    }
}

int main()
{
    // freopen("input.txt", "r", stdin);
    // freopen("output.txt", "w", stdout);

    freopen("ssnd.in", "r", stdin);
    freopen("ssnd.out", "w", stdout);

    ios::sync_with_stdio(0);
    cin.tie(0);
    int t=1;
    // int t;
    // cin >> t;
    for(int tt=1; tt<=t; tt++){
        // cout<<"#Case "<<t<<nl;
        solve();
    }
    return 0;
}