HDU2973 YAPTCHA(威尔逊定理)

计算Sn=k=1n[(3k+6)!+13k+7(3k+6)!3k+7]S_n=\sum^{n}_{k=1}[\frac{(3k+6)!+1}{3k+7}-\lfloor\frac{(3k+6)!}{3k+7}\rfloor]
其中x\lfloor x\rfloor表示不大于xx的最大整数。

威尔逊定理:当且仅当p为素数时,(p1)! 1mod  p( p − 1)! \ ≡−1\mod p
分析:如果3k+7为质数,根据定理,那么(3k+6)!+13k+7\dfrac{(3k+6)!+1}{3k+7}(3k+6)!3k+7\lfloor\dfrac{(3k+6)!}{3k+7}\rfloor相减就为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;
}
posted @ 2019-02-25 22:58  Mr.doublerun  阅读(47)  评论(0)    收藏  举报