# | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
---|---|---|---|---|---|---|---|
21418 | gs14004 | Actual visible points (kriii2_AC) | C++11 | 86 ms | 6708 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
using namespace std;
typedef long long lint;
typedef pair<int, int> pi;
const int mod = 1e9 + 7;
lint ipow(lint x, lint p){
lint ret = 1, piv = x % mod;
while(p){
if(p&1) ret *= piv;
piv *= piv;
ret %= mod;
piv %= mod;
p >>= 1;
}
return ret % mod;
}
lint fact[200005], invf[200005];
lint bino(int x, int y){
return fact[x] * (invf[x-y] * invf[y] % mod) % mod;
}
int n, q, chk[100005];
lint dp[100005];
int mob[100005];
int main(){
fact[0] = invf[0] = 1;
for(int i=1; i<=200000; i++){
fact[i] = fact[i-1] * i % mod;
invf[i] = ipow(fact[i], mod-2);
}
cin >> n >> q;
for(int i=0; i<q; i++){
int x;
scanf("%d",&x);
chk[x] = 1;
}
for(int i=1; i<=100000; i++){
dp[i] += bino(n + i - 2, n - 1);
dp[i] %= mod;
for(int j=2*i; j<=100000; j+=i){
dp[j] += mod - dp[i];
}
}
lint ret = 0;
for(int i=1; i<=100000; i++){
for(int j=i; j<=100000; j+=i){
if(chk[j]){
chk[i] = 1;
break;
}
}
if(chk[i]) ret += dp[i];
}
ret %= mod;
cout << ret << endl;
}
컴파일 시 표준 에러 (stderr) 메시지
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |