Число Фибоначчи 10​^18

#include <iostream>
 
int main(int k, int f)
{
    unsigned long long  Fi, F0, F1;
    int n, i;
    std::cin >> n;
        
    F0 = 0;
    F1 = 1;
    for (i = 2; i <= n; i++)
    {
        Fi = F1 + F0;
        F0 = F1;
        F1 = Fi;
    }
        
        Fi = Fi % 1000000007;
        std::cout << Fi << "\n";
    
 
    return 0;
}

#include <cstddef>
#include <iostream>
 
template <class T, T mod>
class mat22
{
public:
    mat22(T m00, T m01, T m10, T m11) : m_{m00 % mod, m01 % mod, m10 % mod, m11 % mod} { };
    explicit mat22(T n = 0) : mat22{n, 0, 0, n} { }
    T& operator()(std::size_t row, std::size_t col) { return m_[row][col]; }
    const T& operator()(std::size_t row, std::size_t col) const { return m_[row][col]; }
    mat22& operator*=(const mat22& rhs);
private:
    T m_[2][2];
};
 
template <class T, T mod>
mat22<T, mod> operator*(const mat22<T, mod>& lhs, const mat22<T, mod>& rhs)
{
    return {
        lhs(0, 0) * rhs(0, 0) + lhs(0, 1) * rhs(1, 0),
        lhs(0, 0) * rhs(0, 1) + lhs(0, 1) * rhs(1, 1),
        lhs(1, 0) * rhs(0, 0) + lhs(1, 1) * rhs(1, 0),
        lhs(1, 0) * rhs(0, 1) + lhs(1, 1) * rhs(1, 1),
        };
}
 
template <class T, T mod>
mat22<T, mod>& mat22<T, mod>::operator*=(const mat22<T, mod>& rhs)
{
    return *this = *this * rhs;
}
 
template <class T, class U>
T power(T base, U exp)
{
    T result{1};
    while (true /*exp != 0*/)
    {
        if (exp % 2 != 0)
            result *= base;
        exp /= 2;
        if (exp == 0)
            break;
        base *= base;
    }
    return result;
}
 
int main()
{
    unsigned long long n;
    std::cin >> n;
    std::cout << power(mat22<unsigned long long, 1000000007>{1, 1, 1, 0}, n)(0, 1);
}