AC answer

#include <bits/stdc++.h>
using namespace std;

int n;
long long a[10000000];

void work()
{
cin>>n;
a[1]=1,a[2]=2,a[3]=4;
for(int i=4;i<=n;i++)
a[i]=(a[i-1]+a[i-2]+a[i-3])%1000000007;
cout<<a[n];
}

int main()
{
work();
return 0;
}