Pagini recente » Cod sursa (job #1671115) | Cod sursa (job #2711228) | Cod sursa (job #1905214) | Cod sursa (job #886297) | Cod sursa (job #1820692)
#include <fstream>
#include <vector>
#include <bitset>
#define mod 0x26F5
#define MXP 0xF424A
std::vector <int> prim;
std::ofstream t ("ssnd.out");
class parser{
public:
parser() {}
parser(const char *file_name){
input_file.open(file_name);
cursor=0;
input_file.read(buffer,SIZE);}
inline parser &operator >>(int64_t &n){
while(buffer[cursor]<'0' or buffer[cursor]>'9') inc();
n=0;
while('0'<=buffer[cursor] and buffer[cursor]<='9')
n=n*10+buffer[cursor]-'0',inc();
return *this;}
private:
std::ifstream input_file;
static const int SIZE=0x8800;
char buffer[SIZE];
int cursor=0;
inline void inc(){
if(++cursor==SIZE)
cursor=0,input_file.read(buffer,SIZE);
}
};
int64_t pow(int64_t x,int64_t y)
{
x%=mod;
int64_t m=1;
while(y){
if(y&1)
m=(m*x)%mod;
x=(x*x)%mod;
y>>=1;
}
return m;
}
void ciuruieste()
{
std::bitset <MXP> prime;
for (int i=1;((i*i)<<1)+(i<<1)<MXP;++i)
if(!prime[i])
for(int j=((i*i)<<1)+(i<<1);(j<<1)+1<MXP;j+=(i<<1)+1)
prime[j]=true;
prim.push_back(2);
for(int i=1;(i<<1)+1<MXP;++i)
if(!prime[i])
prim.push_back((i<<1)+1);
}
void solve(int64_t target){
int how_many=1,sum=1,aux;
for (std::vector<int>::iterator i=prim.begin();(*i)*(*i)<=target and i!=prim.end();advance(i,1)){
if (target%(*i)==0){aux=0;
while (target%(*i)==0){target/=(*i);++aux;}
how_many*=(1+aux);
sum=(((sum*((pow((*i),1+aux)-1)%mod))%mod)*(pow((*i)-1,mod-2)%mod))%mod;}
}
if (target>1){
how_many<<=1;
sum=(sum*(1+target))%mod;}
t<<how_many<<" "<<sum<<'\n';
}
int main()
{
parser f ("ssnd.in");
int64_t aux,q;
ciuruieste();
f>>q;
for (;q;--q)
f>>aux,solve(aux);
return 0;
}