#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define popcount __builtin_popcount using namespace std; using ll=long long; typedef pair P; template< int char_size > struct TrieNode { int nxt[char_size]; int exist; vector< int > accept; TrieNode() : exist(0) { memset(nxt, -1, sizeof(nxt)); } }; template< int char_size, int margin > struct Trie { using Node = TrieNode< char_size >; vector< Node > nodes; int root; Trie() : root(0) { nodes.push_back(Node()); } void update_direct(int node, int id) { nodes[node].accept.push_back(id); } void update_child(int node, int child, int id) { ++nodes[node].exist; } void add(const string &str, int str_index, int node_index, int id) { if(str_index == str.size()) { update_direct(node_index, id); } else { const int c = str[str_index] - margin; if(nodes[node_index].nxt[c] == -1) { nodes[node_index].nxt[c] = (int) nodes.size(); nodes.push_back(Node()); } add(str, str_index + 1, nodes[node_index].nxt[c], id); update_child(node_index, nodes[node_index].nxt[c], id); } } void add(const string &str, int id) { add(str, 0, 0, id); } void add(const string &str) { add(str, nodes[0].exist); } void query(const string &str, const function< void(int) > &f, int str_index, int node_index) { for(auto &idx : nodes[node_index].accept) f(idx); if(str_index == str.size()) { return; } else { const int c = str[str_index] - margin; if(nodes[node_index].nxt[c] == -1) return; query(str, f, str_index + 1, nodes[node_index].nxt[c]); } } void query(const string &str, const function< void(int) > &f) { query(str, f, 0, 0); } int count() const { return (nodes[0].exist); } int size() const { return ((int) nodes.size()); } }; /** * @brief Aho-Corasick(エイホ–コラシック法) * @docs docs/aho-corasick.md */ template< int char_size, int margin > struct AhoCorasick : Trie< char_size + 1, margin > { using Trie< char_size + 1, margin >::Trie; const int FAIL = char_size; vector< int > correct; void build(bool heavy = true) { correct.resize(this->size()); for(int i = 0; i < this->size(); i++) { correct[i] = (int) this->nodes[i].accept.size(); } queue< int > que; for(int i = 0; i <= char_size; i++) { if(~this->nodes[0].nxt[i]) { this->nodes[this->nodes[0].nxt[i]].nxt[FAIL] = 0; que.emplace(this->nodes[0].nxt[i]); } else { this->nodes[0].nxt[i] = 0; } } while(!que.empty()) { auto &now = this->nodes[que.front()]; int fail = now.nxt[FAIL]; correct[que.front()] += correct[fail]; que.pop(); for(int i = 0; i < char_size; i++) { if(~now.nxt[i]) { this->nodes[now.nxt[i]].nxt[FAIL] = this->nodes[fail].nxt[i]; if(heavy) { auto &u = this->nodes[now.nxt[i]].accept; auto &v = this->nodes[this->nodes[fail].nxt[i]].accept; vector< int > accept; set_union(begin(u), end(u), begin(v), end(v), back_inserter(accept)); u = accept; } que.emplace(now.nxt[i]); } else { now.nxt[i] = this->nodes[fail].nxt[i]; } } } } map< int, int > match(const string &str, int now = 0) { map< int, int > result; for(auto &c : str) { now = this->nodes[now].nxt[c - margin]; for(auto &v : this->nodes[now].accept) result[v] += 1; } return result; } pair< int64_t, int > move(const char &c, int now = 0) { now = this->nodes[now].nxt[c - margin]; return {correct[now], now}; } pair< int64_t, int > move(const string &str, int now = 0) { int64_t sum = 0; for(auto &c : str) { auto nxt = move(c, now); sum += nxt.first; now = nxt.second; } return {sum, now}; } }; const ll MOD=1e9+7; ll f[99]; int main() { int n; ll l, r; cin>>n>>l>>r; f[0]=1, f[1]=1; for(int i=2; i<90; i++){ f[i]=f[i-1]+f[i-2]; } AhoCorasick<10, '0'> aho; for(int i=1; i<90; i++){ if(l<=f[i] && f[i]<=r){ string s; ll x=f[i]; while(x){ s+='0'+(x%10); x/=10; } reverse(s.begin(), s.end()); aho.add(s); } } aho.build(); int m=aho.nodes.size(); vector dp(m); ll ans=MOD-1; dp[0]=1; for(int i=0; i ndp(m); for(int j=0; j=MOD) ndp[l]-=MOD; } } swap(dp, ndp); } for(int i=0; i