1 #include<cstdio>
  2 #include<iostream>
  3 #include<cstring>
  4 #include<algorithm>
  5 #include<cmath>
  6 #define clr(x) memset(x,0,sizeof(x))
  7 #define mod 1000000007
  8 #define LL long long
  9 using namespace std;
 10 int inf[40010];
 11 int prime[10000];
 12 int mart1[2][2];
 13 int mart2[2][2];
 14 int mart3[2][2];
 15 void is_prime();
 16 int eular(int x);
 17 void exgcd(int a,int b,int &x,int &y);
 18 void quickmul(int x);
 19 int fu(int x);
 20 int main()
 21 {
 22     int n;
 23     int xi,yi;
 24     int ans;
 25     is_prime();
 26     while(scanf("%d",&n)!=EOF)
 27     {
 28         ans=0;
 29         if(n==1)
 30         {
 31             printf("2\n");
 32             continue;
 33         }
 34         for(int i=1;i*i<=n;i++)
 35         {
 36             if(i*i==n)
 37                 ans=(int)(((LL)ans+(LL)eular(n/i)*(LL)fu(i))%mod);
 38             else if(n%i==0)
 39                 ans=(int)(((LL)ans+((LL)eular(n/i)*(LL)fu(i))%mod+((LL)eular(i)*(LL)fu(n/i))%mod)%mod);
 40         }
 41         exgcd(n,mod,xi,yi);
 42         xi=(xi%mod+mod)%mod;
 43         ans=(int)((((LL)ans*(LL)xi)%mod+mod)%mod);
 44         printf("%d\n",ans);
 45     }
 46 }
 47 void exgcd(int a,int b,int &x,int &y)
 48 {
 49     if(b==0)
 50     {
 51         x=1;
 52         y=0;
 53         return ;
 54     }
 55     exgcd(b,a%b,y,x);
 56     y-=x*(a/b);
 57     return ;
 58 }
 59 void is_prime()
 60 {
 61     clr(inf);
 62     clr(prime);
 63     inf[1]=inf[0]==1;
 64     int tot=0;
 65     for(int i=2;i<40000;i++)
 66     {
 67         if(!inf[i]) prime[tot++]=i;
 68         for(int j=0;j<tot;j++)
 69         {
 70             if(i*prime[j]>40000) break;
 71             inf[i*prime[j]]=1;
 72             if(i%prime[j]==0) break;
 73         }
 74     }
 75     return ;
 76 }
 77 int eular(int x)
 78 {
 79     int ans=x;
 80     for(int j=0;prime[j]*prime[j]<=x;j++)
 81     if(x%prime[j]==0)
 82     {
 83         ans-=ans/prime[j];
 84         while(x%prime[j]==0)
 85             x/=prime[j];
 86     }
 87     if(x>1) ans-=ans/x;
 88     return ans;
 89 }
 90 int fu(int x)
 91 {
 92     if(x==1) return 1;
 93     if(x==2) return 3;
 94     clr(mart1);
 95     clr(mart2);
 96     clr(mart3);
 97     mart1[1][1]=mart1[1][0]=mart1[0][1]=mart2[1][1]=mart2[0][0]=1;
 98     quickmul(x-2);
 99     return (int)(((LL)mart2[1][0]*1+(LL)mart2[1][1]*3)%mod);
100 }
101 void quickmul(int x)
102 {
103     int d;
104     while(x)
105     {
106         if(x&1)
107         {
108             for(int i=0;i<2;i++)
109                 for(int j=0;j<2;j++)
110                 {
111                     d=0;
112                     for(int k=0;k<2;k++)
113                         d=(int)(((LL)d+((LL)mart2[i][k]*(LL)mart1[k][j])%mod)%mod);
114                     mart3[i][j]=d;
115                 }
116             for(int i=0;i<2;i++)
117                 for(int j=0;j<2;j++)
118                     mart2[i][j]=mart3[i][j];
119         }
120         x>>=1;
121         for(int i=0;i<2;i++)
122             for(int j=0;j<2;j++)
123             {
124                 d=0;
125                 for(int k=0;k<2;k++)
126                     d=(int)(((LL)d+((LL)mart1[i][k]*(LL)mart1[k][j])%mod)%mod);
127                 mart3[i][j]=d;
128             }
129             for(int i=0;i<2;i++)
130                 for(int j=0;j<2;j++)
131                     mart1[i][j]=mart3[i][j];
132     }
133     return ;
134 }