-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathcf422.cpp
More file actions
89 lines (82 loc) · 2.01 KB
/
cf422.cpp
File metadata and controls
89 lines (82 loc) · 2.01 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define Debug(x) cout << "> " << #x << " : " << x << "\n";
#define DebugArr(a,n) For(i,0,n) cout << i << " : " << a[i] << "\n";
#define For(i,a,b) for(int i=(a);i<(b);i++)
#define ForD(i,a,b) for(int i=(a);i>(b);i--)
#define tests int t; cin>>t; while(t--)
#define Imax INT_MAX
#define Imin INT_MIN
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef pair<int,int> pii;
#define sz(a) int((a).size())
#define pb push_back
#define it iterator
#define all(c) (c).begin(),(c).end()
#define tr(c,i) for(typeof((c).begin() i = (c).begin(); i != (c).end(); i++)
#define present(c,x) ((c).find(x) != (c).end())
#define cpresent(c,x) (find(all(c),x) != (c).end())
#define mod 1000000007
/////////////////////////////////////////////////////////////////////
std::vector<int> v(5001000, 0);
long long int mpow(long long int x,long long int n,long long int p) {
long long int res = 1;
x %= p;
while(n > 0) {
if(n & 1) res = (res*x*1LL) % p; // if y is odd
n = n>>1; // n/2
x = (x*x*1LL) % p;
}
return res;
}
void pre(int r){
v[1] = 0;v[2] = 1;v[3] = 3;
long long int best,tp;
int numg;
For(i,4,r+1){
int up = sqrt(i);up += 2;
best = i*(i-1);
best = best/2;//base
For(j,2,up){
if(i%j == 0){
//form group of size j
numg = i/j;
tp = i*(j-1);
tp = tp/2;
tp += v[numg];
best = min(best,tp);
//form group of size i/j
tp = i*(numg-1);
tp = tp/2;
tp += v[j];
best = min(best,tp);
}
}
v[i] = best;
// Debug(v[i])
}
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("input.txt", "r", stdin);
freopen("output.txt", "w", stdout);
#endif
ios_base::sync_with_stdio(false);
cin.tie(NULL);cout.tie(NULL);
long long int t,l,r,ans=0,step;
cin >> t >> l >> r;
pre(r);
For(i,l,r+1){
step = v[i];
step = step%mod;
step *= mpow(t,i-l,mod);
step = step%mod;
ans += step;
ans = ans%mod;
}
cout << ans;
return 0;
}