poj2049Let it Bead

#include<stdio.h>
#include<string.h>
#include<math.h>
#define inta __int64
int gcd(int x,int y){
    int r;
    if(y==0) return x;
    else r=x%y;
    return gcd(y,r);
}
int poww(int a,int b){
    int i,tmp;
    tmp=1;
    if(b==0) return 1;
    if(b==1) return a;
    for(i=1;i<=b;i++)
        tmp*=a;
    return tmp;
}
int main(){
    int n,m;
    int ans;
    while(scanf("%d %d ",&n,&m)!=EOF&&(n!=0||m!=0)){

        ans=0;
        for(int i=1;i<=m;i++){
            ans+=poww(n,gcd(i,m));
        }
        if(m&1){
            ans+=m*poww(n,m/2+1);
        }
        else{
            ans+=m/2*poww(n,m/2)+m/2*poww(n,m/2+1);
        }
        ans/=m*2;
        printf("%d\n",ans);
    }
    return 0;
}