#include #include #include #define MOD 1000000007 typedef unsigned long long llu; typedef signed long long lld; llu n; void mtrMultiple( /* const */ llu* left[], /* const */ llu* right[], int sizeA, int sizeB, int sizeC, llu* resultVar[]){ int i, j, k; llu resultdata[sizeA][sizeC]; for (i = 0; i < sizeA; i++) { for (j = 0; j < sizeC; j++) { llu cell = 0; for (k = 0; k < sizeB; k++) { cell += (left[i][k] * right[k][j]) % MOD; cell %= MOD; } resultdata[i][j] = cell; } } for (i = 0; i < sizeA; i++) { for (j = 0; j < sizeC; j++) { resultVar[i][j] = resultdata[i][j]; } } } void powmtr(/* const */ llu* mtr[], llu n, int size, llu* resultVar[]){ int i, j; llu prevmtr[size][size]; llu resultdata[size][size]; llu* argMtr[size]; llu* argPrm[size]; llu* argResd[size]; for (i = 0; i < size; i++) { argMtr[i] = mtr[i]; argPrm[i] = prevmtr[i]; argResd[i] = resultdata[i]; } for (i = 0; i < size; i++) { for (j = 0; j < size; j++) { resultdata[i][j] = (i == j ? 1 : 0); prevmtr[i][j] = mtr[i][j]; } } while (n) { if (n & 1) { mtrMultiple((llu**)argResd, (llu**)argPrm, size, size, size, (llu**)argResd); } n = (n >> 1); mtrMultiple((llu**)argPrm, (llu**)argPrm, size, size, size, (llu**)argPrm); } for (i = 0; i < size; i++) { for (j = 0; j < size; j++) { resultVar[i][j] = resultdata[i][j]; } } } void solve(){ int i; llu leftMatrix[1][2] = { {1, 0} }; llu rightMatrix[2][2] = { {1, 1}, {1, 0} }; llu resultMatrix[1][2]; llu answer; llu* argLeftm[1] = {leftMatrix[0]}; llu* argRightm[2] = {rightMatrix[0], rightMatrix[1]}; llu* argResm[1] = {resultMatrix[0]}; powmtr((llu**)argRightm, n, 2, (llu**)argRightm); mtrMultiple((llu**)argLeftm, (llu**)argRightm, 1, 2, 2, (llu**)argResm); answer = resultMatrix[0][0] * resultMatrix[0][1] % MOD; printf("%llu\n", answer); return; } int main(void){ scanf("%llu", &n); solve(); return 0; }