int gcd(int x,int y)
{
if(y==0)
return x;
return gcd(y,x%y);
}
ll rolya()
{
int i;
ll sum=0;
repf(i,0,n-1)
sum+=(ll)pow(1.0*col,gcd(n,i));
if(n%2==0)
{
sum+=(ll)n/2*pow(1.0*col,n/2);
sum+=(ll)n/2*pow(1.0*col,n/2+1);
}
else
sum+=(ll)n*pow(1.0*col,n/2+1);
return sum/2/n;
}