#include #include #include #include using namespace std; const int mod = 1e9 + 7; struct Modint { int n; Modint(int n = 0) : n(n) {} }; Modint operator+(Modint a, Modint b) { return Modint((a.n += b.n) >= mod ? a.n - mod : a.n); } Modint operator-(Modint a, Modint b) { return Modint((a.n -= b.n) < 0 ? a.n + mod : a.n); } Modint operator*(Modint a, Modint b) { return Modint(1LL * a.n * b.n % mod); } Modint &operator+=(Modint &a, Modint b) { return a = a + b; } Modint &operator-=(Modint &a, Modint b) { return a = a - b; } Modint &operator*=(Modint &a, Modint b) { return a = a * b; } int main() { int n; long long L, R; cin >> n >> L >> R; vector> dp(n + 1, vector(63)); Modint bad; dp[0][0] = 1; for (int i = 0; i < n; i++) { for (int j = 0; j <= 62; j++) { // update dp[i + 1][min(61, j + 1)] += dp[i][j]; // keep if (i != n - 1) { dp[i + 1][j] += dp[i][j] * i; } else { bad += dp[i][j] * i; } } } Modint ans; if (L == 0) { ans += bad; } for (int i = 1; i <= 61; i++) { long long way = 1LL << (i - 1); if (L <= way && way <= R) { ans += dp[n][i]; } } cout << ans.n << endl; }