计算
S
n
=
∑
k
=
1
n
[
(
3
k
+
6
)
!
+
1
3
k
+
7
−
⌊
(
3
k
+
6
)
!
3
k
+
7
⌋
]
S_n=\sum^{n}_{k=1}[\frac{(3k+6)!+1}{3k+7}-\lfloor\frac{(3k+6)!}{3k+7}\rfloor]
Sn=k=1∑n[3k+7(3k+6)!+1−⌊3k+7(3k+6)!⌋]
其中
⌊
x
⌋
\lfloor x\rfloor
⌊x⌋表示不大于
x
x
x的最大整数。
威尔逊定理:当且仅当p为素数时,
(
p
−
1
)
!
≡
−
1
m
o
d
  
p
( p − 1)! \ ≡−1\mod p
(p−1)! ≡−1modp
分析:如果3k+7为质数,根据定理,那么
(
3
k
+
6
)
!
+
1
3
k
+
7
\dfrac{(3k+6)!+1}{3k+7}
3k+7(3k+6)!+1和
⌊
(
3
k
+
6
)
!
3
k
+
7
⌋
\lfloor\dfrac{(3k+6)!}{3k+7}\rfloor
⌊3k+7(3k+6)!⌋相减就为1,否则这两个数差为0。
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int N = 1e6;
int isprime[N];
int s[N];
void init(){
memset(isprime,0,sizeof(isprime));
for(int i = 2; i < 3*N; i++)
{
if(!isprime[i])
for(int j = 2*i; j < 3*N; j += i)
{
isprime[j] = 1;
}
}
s[0] = 0,s[1] = 0;
for(int i = 2; i <= N; i++)
{
int tmp = 3 * i + 7;
s[i] = s[i-1] + (1 - isprime[tmp]);
}
}
int main(){
int t;
init();
scanf("%d",&t);
while(t--){
int n;
scanf("%d",&n);
printf("%d\n",s[n]);
}
return 0;
}