HDU- 2973 YAPTCHA(威尔逊定理)

题目链接:点击查看

题目大意:求解:

S_{n}=\sum_{k=1}^{n}\left[ \frac{(3k+6)!+1}{3k+7} - \left[ \frac{(3k+6)!}{3k+7}\right] \right]

题目分析:我们需要提前知道威尔逊定理:如果p是素数,则(p-1)!\equiv -1(mod\ p)

所以我们可以分情况讨论:

  1. 当3k+7是素数时,\frac{(3k+6)!+1}{3k+7}=\frac{(3k+6)!}{3k+7}+\frac{1}{3k+7},因为威尔逊定理,(3k+6)!\equiv -1(mod\ (3k+7)),所以\frac{(3k+6)!+1}{3k+7}一定是一个整数,我们设为x,这样一来后面的\frac{(3k+6)!}{3k+7}=x-\frac{1}{3k+7},因为\frac{1}{3k+7}一定小于等于1,所以\frac{(3k+6)!}{3k+7}取整后一定是x-1,这样一来f(k)=1
  2. 当3k+7不是素数时,(3k+6)!一定能被(3k+7)整除(显然成立),然后上述公式就转换成了\left[\frac{(3k+6)!}{3k+7} +\frac{1}{3k+7}- \left[ \frac{(3k+6)!}{3k+7}\right]\right]=\left[\frac{1}{3k+7}\right]=0,所以f(k)=0

综上所述,先用欧拉线性筛打一个3e6的素数表,因为3k+7最大能到3e6+7,然后预处理一个前缀和,就可以O(1)回答每一个查询了

代码:

#include<iostream>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<climits>
#include<cmath>
#include<cctype>
#include<stack>
#include<queue>
#include<list>
#include<vector>
#include<set>
#include<map>
#include<sstream> 
using namespace std;

typedef long long LL;

const int inf=0x3f3f3f3f;

const int N=1e7+100;

int cnt=0;

bool vis[N];

int pri[N];

int sum[N];

void P()
{
	memset(vis,false,sizeof(vis));
	for(int i=2;i<N;i++)
	{
		if(!vis[i])
			pri[cnt++]=i;
		for(int j=0;j<cnt&&pri[j]*i<N;j++)
		{
			vis[pri[j]*i]=true;
			if(i%pri[j]==0)
				break;
		}
	}
	memset(vis,false,sizeof(vis));
	for(int i=0;i<cnt;i++)
		vis[pri[i]]=true;
}

void init()
{
	P();
	sum[0]=0;
	for(int i=1;i<1e6+10;i++)
	{
		sum[i]=sum[i-1];
		if(vis[3*i+7])
			sum[i]++;
	}
}

int main()
{
//	freopen("input.txt","r",stdin);
	init();
	int w;
	cin>>w;
	while(w--)
	{
		int n;
		scanf("%d",&n);
		printf("%d\n",sum[n]);
	}
	
	
	
	
	
	
	
	
	
	
	

	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Frozen_Guardian

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值