#pragma GCC optimize("Ofast")
#include <bits/stdc++.h>
#include <iostream>
#include <iomanip>
#include <cassert>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <vector> 
#include <chrono>
#include <algorithm>
using namespace std;
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace __gnu_pbds;
template <typename T>
using ordered_set = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
// order_of_key(val): returns the number of values less than val
// find_by_order(k): returns the iterator to kth largest element.(0-based)
#define int long long
#define pb push_back
// #define pb emplace
#define endl "\n"
#define all(con) con.begin(), con.end()
#define fi first
#define se second
#define F first
#define S second
#define MP make_pair
#define sz(x) (int)x.size()
#define rp(i,n) for(int i=0;i<n;i++)
#define fr(i,n) for(int i=1;i<=n;i++)
#define ff(i,a,b) for(int i=a;i<=b;i++)
#define fff(i,a,b) for(int i=a;i<b;i++)
#define TRACE
#ifndef ONLINE_JUDGE
#define trace(...) __f(#__VA_ARGS__, __VA_ARGS__)
template <typename Arg1>
void __f(const char *name, Arg1 &&arg1){cerr<<name<<" : "<<arg1<<endl;}
template <typename Arg1, typename... Args>
void __f(const char *names, Arg1 &&arg1, Args &&... args){const char *comma = strchr(names + 1, ',');cerr.write(names, comma - names) << " : " << arg1 << " | ";__f(comma + 1, args...);}
#else
#define trace(...)
#endif
typedef long double ld;
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<string> vs;
typedef vector<char> vc;
typedef vector<pair<int,int>> vpi;
const int mod = 1e9 + 7;
const int mod1 = 998244353;
const int inf = 2e18;
// const int MOD = 998244353;
const int MOD = 1e9+7;
int gcd(int a,int b){return __gcd(a,b);}
int lcm(int a,int b){return (a*b)/(__gcd(a,b));}
int hcf(int a,int b){if(a==0)return b;return hcf(b%a, a);}
int modpow(int a,int b,int m){int ans=1;while(b){if(b&1)ans=(ans*a)%m;b/=2;a=(a*a)%m;}return ans;}
bool isprime(int n){if(n<2)return false;if(n==2)return true;if(n%2==0)return false;for(int i=3;i*i<=n;i+=2){if(n%i==0)return false;}return true;}
const int MXX=1e5+69;
int fact[MXX],invfact[MXX];
int modinv(int k){return modpow(k,MOD-2,MOD);}
void combinatorics(){fact[0]=fact[1]=1;for(int i=2;i<MXX;i++){fact[i]=fact[i-1]*i;fact[i]%=MOD;}invfact[MXX-1]=modinv(fact[MXX-1]);for(int i=MXX-2;i>= 0;i--){invfact[i]=invfact[i+1]*(i+1);invfact[i] %= MOD;}}
int nCr(int x, int y){if(y>x)return 0;int num=fact[x];num*=invfact[y];num%=MOD;num*=invfact[x-y];num%=MOD;return num;}
string lwr(string s){transform(all(s),s.begin(),::tolower);return s;}
string upr(string s){transform(all(s),s.begin(),::toupper);return s;}
struct DSU{vector<int> par,sz;DSU(int n){par.resize(n);sz.resize(n);for(int i=0;i<n;i++)par[i]=i,sz[i]=1;}
int find(int x){if(x == par[x])return x;else return (par[x]=find(par[x]));}
void unite(int x,int y){x=find(x);y=find(y);if(x==y)return;if(sz[x]<sz[y])swap(x,y);sz[x]+=sz[y];par[y]=x;}
bool Same(int x, int y){return find(x)==find(y);}int size(int x){return sz[find(x)];}};
template<typename T> void prev(vector<T> &v){for(auto t:v)cout<<t<<" ";cout<<endl;}
void prar(int a[],int n){for(int i=0;i<n;i++){cout<<a[i]<<" ";}cout<<endl;}
template <typename T> std::pair<bool, int > findInVector(const std::vector<T>  & vecOfElements, const T  & element)
{std::pair<bool, int > result;auto it = std::find(vecOfElements.begin(), vecOfElements.end(), element);if (it != vecOfElements.end()){result.second = distance(vecOfElements.begin(), it);result.first = true;}
else{result.first = false;result.second = -1;}return result;}
bool isvowel(char ch){string s="aeiou";for(char c:s){if(c==ch){return true;}}return false;}
void dream();
int32_t main(){
ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
srand(chrono::high_resolution_clock::now().time_since_epoch().count());
cout << setprecision(8);
// #ifndef ONLINE_JUDGE // Comment this if the question demands input-output through a file
// freopen("in.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
// #endif
// combinatorics(); // Take care of MOD
int ttt = 1;
// cin>>ttt;
for(int i=1;i<=ttt;i++) dream();
return 0;
}
int f(int a,int b,int c){
    return gcd(a,gcd(b,c));
}
void dream(){
    int n;
    cin>>n;
    int ans=0;
    fr(i,n){
        fr(j,n){
            fr(k,n){
                ans += f(i,j,k);
            }
        }
    }
    cout<<ans<<endl;
	// you should actually read the stuff at the bottom
}
/* stuff you should look for
	* int overflow, array bounds
	* special cases (n=1?)
    * boundary cases
	* do something instead of nothing and stay organized
*/
/* some basic stuff that you should learn:
        * Sets
            * begin() - returns an iterator to first element of set
            * end() - Returns an iterator to the theoretical element that follows last element in the set.
            * size() - Returns the number of elements in the set.
            * max_size() - Returns the maximum number of elements that the set can hold.
            * empty() - Returns whether the set is empty.
            * find(x) - Returns an iterator to the element ‘x’ in the set if found, else returns the iterator to end.
            * erase(x)– Removes the value ‘x’ from the set.
            * count(x) - Returns 1 or 0 based on the element ‘x’ is present in the set or not.
            * lower_bound(x) – Returns an iterator to the first element that is equivalent to ‘x’ or 
                               definitely will not go before the element ‘x’ in the set.
            * Return reverse iterator to reverse beginning: (s.rbegin())
            * Set to be sorted in decreasing order: set<int, greater<int> > 
            * s.count(x) : check if an element exists in a set #returns 0 or 1 accordingly
            * emplace() : If the function successfully inserts the element (because no 
                    equivalent element existed already in the set), the function returns a 
                    pair of an iterator to the newly inserted element and a value of true.
            * const bool is_in = s.find(x) != s.end(); check if an element is in the set. 
            * myset.erase(x) #directly by value. logarithmic complexity
            * myset.erase(iteratortoanelement,s.end())
            * for_each(all(s),do_something);
            * set<int> :: iterator it = myset.find(60); 
        
        * Vectors 
            * begin() - Returns an iterator pointing to the first element in the vector
            * end() - Returns an iterator pointing to the theoretical element that follows the last element in the vector
            * rbegin() - Returns a reverse iterator pointing to the last element in the vector (reverse beginning).
            * rend() - Returns a reverse iterator pointing to the theoretical element preceding the first element in the vector
            * cbegin(),cend(),crbegin(),crend()
            * size() - Returns the number of elements in the vector.
            * max_size() -  Returns the maximum number of elements that the vector can hold.
            * capacity() - Returns the size of the storage space currently allocated to the vector expressed as number of elements.
            * resize() - Resizes the container so that it contains ‘n’ elements.
            * empty() - Returns whether the container is empty.
            * shrink_to_fit() - Reduces the capacity of the container to fit its size and destroys all elements beyond the capacity. 
            * reserve() - Requests that the vector capacity be at least enough to contain n elements.
            * [i] - Returns a reference to the element at position ‘i’ in the vector
            * at(i) - Returns a reference to the element at position ‘i’ in the vector
            * front() - Returns a reference to the first element in the vector
            * back() - Returns a reference to the last element in the vector
            * data() - Returns a direct pointer to the memory array used internally by the vector to store its owned elements.
            * assign() - It assigns new value to the vector elements by replacing old ones
            * push_back() - It push the elements into a vector from the back
            * pop_back() - It is used to pop or remove elements from a vector from the back.
            * insert() - It inserts new elements before the element at the specified position
            * erase() - It is used to remove elements from a container from the specified position or range
            * swap() - It is used to swap the contents of one vector with another vector of same type. Sizes may differ.
            * clear() - It is used to remove all the elements of the vector container
            * emplace() - It extends the container by inserting new element at position
            * emplace_back() - It is used to insert a new element into the vector container, the new element is added to the end of the vector
            * for (auto it = myvector.begin(); it != myvector.end(); ++it) 
        
        * Queue
            * empty() - Returns whether the queue is empty.
            * size() - Returns the size of the queue.
            * emplace() - Insert a new element into the queue container, the new element is added to the end of the queue.
            * front() - front() function returns a reference to the first element of the queue.
            * back() - returns a reference to the last element of the queue.
            * push(x) - adds the element ‘x’ at the end of the queue
            * pop() - deletes the first element of the queue
            
        * Priority_Queue
            * priority_queue <int, vector<int>, greater<int>>
            * empty() - returns whether the queue is empty.
            * size() -  returns the size of the queue.
            * top() -  Returns a reference to the top most element of the queue
            * push(x) - adds the element ‘x’ at the end of the queue.
            * pop() -  deletes the first element of the queue.
            * swap() -  swap the contents of one priority queue with another priority queue of same type and size.
            * emplace() - insert a new element into the priority queue container, the new element is added to the top of the priority queue.
        
        * Deque
            Double ended queues are sequence containers with the feature of expansion and contraction on both the ends.
            They are similar to vectors, but are more efficient in case of insertion and deletion of elements. 
            Unlike vectors, contiguous storage allocation may not be guaranteed.
            * insert()
            * rbegin()
            * rend()
            * assign()
            * resize()
            * push_front()
            * push_back()
            * pop_front()
            * pop_back()
            * front()
            * back()
*/

Measured memory_bytes: 3632