题解
A - Prime Gap
题目:读入一个n,是素数就输出0,不是就输出它左右两边最近的两个素数形成的区间的长度。
- 素数线性筛
- 思路:埃式筛法打表(素数表),用结合比较快的素数判定函数。
- 注意点:上限是1299709,可以设成
const int maxn = 1299709+100;
- 判定素数函数部分
bool isprime(int n)
{
if(n < 2)return 0;
if(n ==2 || n == 3)return 1;
if(n % 6 != 1 && n % 6 != 5)return 0;
for(int i = 5; i * i <= n; i += 6)
{
if(n % i == 0 || n % (i + 2) == 0)
return 0;
}
return 1;
}
- 打表部分(埃式)
//构造素数表
memset(bl, 0, sizeof bl);
for(int i=2; i<maxn; ++i)
for(int k=i*2; k<maxn; k+=i)
bl[k]=1;
int num=0;
for(int i=2; i<maxn; ++i)//素数队列
if(bl[i]==0)
prime[num++]=i;
- 逻辑部分
int n;
while(cin>>n)
{
if(n==0)
break;
if(isprime(n))
{
cout<<0<<endl;
continue;
}
for(int i=0; i<num; ++i)
{
if(prime[i]<n&&prime[i+1]>n)
cout<<(prime[i+1]-prime[i])<<endl;
}
}
B - Relatives
题目:求单个欧拉函数
- 求单个欧拉函数
- 样题
#include <iostream>
using namespace std;
int eular(int n) //求单个数的欧拉函数
{
int ans = n;
for(int i = 2; i*i <= n; i++)
{
if(n % i == 0)
{
ans = ans/i*(i-1); //ans=ans*(1-1/i)=ans*(i-1)/i
// =ans/i*(i-1) //先除防止溢出
while(n % i == 0) //消除i因子
n /= i;
}
}
if(n > 1)ans = ans/n*(n-1); //最后可能还剩下一个质因数没有除
return ans;
}
int main()
{
int n;
while(cin >> n && n)
{
cout << eular(n) << endl;
}
return 0;
}
C - Sum of Consecutive Prime Numbers
题目:读入一个正整数n,求有多少种方式可以把它表示成几个连续的素数的和。
- 素数线性筛(一个正整数转换多个连续的素数的和)
- 打表(埃式筛法)+暴力枚举
- 注意点:几个
连续的
素数的和 - 还用A题的打表模板,下面是逻辑部分
int n;
while(cin>>n&&n)
{
int cnt(0);
for(int i=0; prime[i]<=n; i++)
{
int sum(0);
for(int j=i; ; j++)
{
if(sum == n)
{ cnt++; break; }
else if(sum > n) break;
else sum += prime[j];
}
}
cout << cnt << endl;
}
D - Farey Sequence
题目:读入正整数n,求1到n里面有几对互质数。(如果 g c d ( a , b ) = 1 gcd(a,b)=1 gcd(a,b)=1,则称a,b为互质数)
- 欧拉函数
- 思路:
- 分析:找到所有的互质数对,从2开始每一个数找它前面的存在的互质数,恰好符合欧拉函数的特点。
- 做法:构造欧拉表(打表)
- 构造欧拉表(埃式;模板)。直接函数读入
maxn
的值就可以打表了。
void euler(int n)
{
for (int i=1; i<=n; i++) phi[i]=i;
for (int i=2; i<=n; i++)
{
if (phi[i]==i)//这代表i是质数
{
for (int j=i; j<=n; j+=i)
{
phi[j]=phi[j]/i*(i-1); //把i的倍数更新掉
}
}
}
}
- 另一种方法
memset(a,0,sizeof a);
//构造欧拉表
for(int i=2; i<=mx; i++)
{
if(a[i]==0)
{
for(int j=i; j<=mx; j+=i)
{
if(a[j]==0)
a[j]=j;
a[j]=a[j]/i*(i-1);
}
}
}
E - Visible Lattice Points
题目;方格网中取 ( n + 1 ) ∗ ( n + 1 ) (n+1)*(n+1) (n+1)∗(n+1)个坐标点,左下角为(0,0),求从(0,0)点可以看到几个点。(从原点看去,在同一条线上的点只能看见最前面的点)
- 多个数的欧拉函数
- 思路:
- 分析:可以知道,可见的点的总数,等于按照对角线分开的其中一个三角形的可见点的二倍加上斜率为1的那条线上的一个点。对一个三角形分析可知,n增加的时候,点数增加的数对应着新的斜率。如图。
(图片上有点错误,n=4时,不是两个1/3,而是一个1/3一个2/3) - 做法:欧拉表加输出调整;按照欧拉函数求那些斜率大于0小于1的线对应的点数的二倍,最后加3(一个斜率为0,一个斜率为1,一个没有斜率)
- 对应代码
//n*n的方格,从(0,0)点可以看到几个点
#include <iostream>
using namespace std;
const int mx = 1000;//开大了也可能出现run的问题问题
int phi[mx+10];
void euler(int n)
{
for (int i=1; i<=n; i++) phi[i]=i;
for (int i=2; i<=n; i++)
{
if (phi[i]==i)//这代表i是质数
{
for (int j=i; j<=n; j+=i)
{
phi[j]=phi[j]/i*(i-1); //把i的倍数更新掉
}
}
}
}
int main()
{
int n, t;
euler(mx);
int cnt=1;
cin >> t;
for(; cnt<=t; cnt++)
{
cin >> n;
int sum(0);
for(int i=2; i<=n; i++)
sum+=phi[i];
//就是改了个输出,分析可知n的增加,对应着增加了1/n,2/n,3/n````(n-1)/n,而且
//其中的假分数代表被挡到,看不见了,所以还是欧拉数列,每个n是前面那个数对应的欧拉值的和
//再加上斜率为0和斜率为1的三个点
cout << cnt << " " << n << " " <<sum*2+1+2<< endl;
}
return 0;
}
F - Almost Prime Numbers
题目:定义
almost prime
就是除了1和本身仅可以被一个素数整除的非素数,求一个给定区间里的almost prime
数量。
如4,8就是符合条件的数
- 素数线性筛
- 思路:打表+暴力枚举
- 对应代码
- 注意点:这里的
long long
很重要
//自己的按照a题弄的,错在了没用ll,一直没发现
#include<iostream>
#include<string.h>
using namespace std;
#define ll long long
const int maxn = 1000005;
ll prime[maxn],vis[maxn];
int main()
{
//构造素数表
memset(vis, 0, sizeof vis);
for(ll i=2; i<=maxn; i++)
for(ll k=i*2; k<=maxn; k+=i)
vis[k]=1;
ll len=0;
for(ll i=2; i<=maxn;i++)//素数队列
if(vis[i]==0)
prime[len++]=i;
ll t;
cin>>t;
while (t--)
{
ll cnt = 0;
ll l,r;
cin>>l>>r;
for (ll i = 0; i < len; i++)
for ( ll j = prime[i]*prime[i]; j <= r; j*=prime[i] )
if ( j >= l )
cnt++;
cout<<cnt<<endl;
}
return 0;
}