이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#pragma warning(disable:4996)
#include<stdio.h>
#include<algorithm>
using namespace std;
int n, m, w[101000];
long long Mod = 1000000007, F[201000], G[201000], Res, invF[201000];
bool chk[101000];
long long Inv(long long a){
long long cnt = Mod - 2, r = 1;
while (cnt){
if (cnt % 2){
r = r*a%Mod;
}
cnt /= 2;
a = a*a%Mod;
}
return r;
}
long long Gap(int a){
if (a == 0)return 0;
return F[a + n - 1] * invF[n] % Mod * invF[a - 1] % Mod;
}
void Do(int a){
G[a] = (G[a] + Gap(a) - Gap(a-1) + Mod) % Mod;
int i, ck = chk[a];
for (i = a * 2; i <= 100000; i += a){
G[i] = (G[i] + Mod - G[a]) % Mod;
if (chk[i])ck = true;
}
if (ck)Res = (Res + G[a]) % Mod;
}
int main()
{
int i;
scanf("%d%d", &n, &m);
for (i = 1; i <= m; i++){
scanf("%d", &w[i]);
chk[w[i]] = true;
}
F[0] = 1;
for (i = 1; i <= 200010; i++){
F[i] = F[i - 1] * i%Mod;
invF[i] = Inv(F[i]);
}
for (i = 1; i <= 100000; i++){
Do(i);
}
printf("%lld\n", Res);
}
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |