素数筛
1.朴素算法
1.1试除法
使用试除法判断n是否为素数,则在[2,n)范围内逐一枚举数字,尝试用n去除。
- 若能被范围内任意数整除,则说明枚举的数字时n的因子。素数n非素数。
- 若范围内没有能被n整除的数字,则n为素数
代码:
int isprime(int n)
{
int i;
if(n<2)//最小的素数是2,因此小于2都不是素数
return 0;
for(i=2;i<n;i++)//枚举[2,n)内有无n的因子
if(n%2==0)
return 0;
return 1;
}
1.2优化
上述方法试除范围过大,一旦n的规模变大就需要耗费很长时间,所以可以从试除范围入手进行优化。
数字的因子有个特点就是成对出现(两个成对的因子相乘等于原来的数),且成对的因子会分布在sqrt(n)的两侧。假设k=sqrt(n),那么k*k=n。如果两个因子不分布在k的两侧,那么两个因子相乘要么大于n要么小于n,既然相乘结果不等于n那么就不是n的因子,与开始的假设相违背。
综上,可将试除范围缩小到[2,sqrt(n)]。(如果[2,sqrt(n)]内没有因子,那么[sqrt(n),n-1]也不会有因子,因为因子是成对出现的)
代码如下:
int isprime(int n)
{
int i,x=sart(n);
if(n<2)//最小的素数是2,因此小于2都不是素数
return 0;
for(i=2;i<=x;i++)//枚举[2,sqrt(n)]内有无n的因子
if(n%2==0)
break;
if(i>x)
return 1;
else
return 0;
}
2.埃氏筛
反向构造
我们事先不知道一个数是否素数,但可以构造合数。那么当我们把某个范围内的合数标记出来了,那剩下的数就是素数。
算法
1.定义两个数组,一个数组包括范围内所有数字,用来标记该数字是否被访问。另外一个数组用来存储已经筛选出来的素数。
2.将数组的标记设为未访问,将0,1设为已访问。
3.从2开始循环直到范围上界,判断每一个数是否被访问过,访问过的是合数,未访问是素数。
4.对该数字进行倍增,从两倍开始,直至若干倍到达上界为止。其中得到的结果都标记为已访问。
代码如下:
#include<stdio.h>
#define N 1000001
int primes[N]={0}, pri[N] = {0};
int x = 1000000,pp=0;
//埃氏筛法
void Eth()
{
pri[0]=1,pri[1]=1;
for (int i = 2; i <= x; i++)
{
if (!pri[i])//未被访问且最小
primes[++pp] = i;//将素数存入primes数组中
for (int j = 2; i*j <= x; j++)//标记合数
pri[j* i] = 1;//将倍数标记为合数
}return;
}
void main()
{
int n,k;
Eth();
scanf("%d", &n);
for (k = 1; primes[k] <= n&&primes[k]!=0&&k<=pp; k++)
printf("%d\n", primes[k]);
}
3.欧拉筛
埃氏筛存在一个问题:一个数字可能会被重复筛除。比如12可能在2的时候被筛,也可能在6的时候被筛。
欧拉筛可以很好的解决这个问题,保证每个合数只被筛一次,从而提高效率。
根据唯一分解定理,每个合数都有一个最小质因子,保证每个合数都被其最小质因子筛去就🆗了。
代码如下:
#include<stdio.h>
#define N 1000001
int primes[N]={0}, pri[N] = {0};
int x = 1000000,pp=0;
//欧拉筛法
void ola()
{
pri[0]=1,pri[1]=1;
for (int i = 2; i <= x; i++)
{
if (!pri[i])//为0未被标记则为素数
primes[++pp] = i;//将素数存入primes数组中
for (int j = 1; primes[j] * i <= x; j++)//标记合数 条件判断防止数组越界
{
pri[primes[j] * i] = 1;//将i的质因数倍标记为合数
if (i % primes[j] == 0)//判断primes[j]是否为i的最小质因数,如果是则标记合数到此为止
break;
}
}return;
}
void main()
{
int n,k;
ola();
scanf("%d", &n);
for (k = 1; primes[k] <= n&&primes[k]!=0&&k<=pp; k++)
printf("%d\n", primes[k]);
}
埃氏筛是,只要是素数的倍数就倍增。而欧拉筛是用当前遍历到的数字i,去乘以已经在素数表中的素数。 欧拉筛和埃氏筛的框架大致相同,区别在于第二层循环的倍增过程的操作。
Problem D of 1725 :Jack的字符串问题
题目描述
阿操最讨厌写字符串的题目,看到眼前密密麻麻的字母就烦。这件事的起因就是很久很久以前的一次找重复字符的位置。 现在有一个字符串,我们要找出其中的重复的字符并输出这些字符和字符的位置,如:aabcaabc22 输出 a,0;a,1;a,4;a,5,b,2;b,6,c,3;c,7;2,8;2,9
输入
输入一行字符串(字符串中只含数字和字母)。其长度不超过100。包括多组输入。
输出
根据样例的格式将重复出现的字符位置输出
样例输入
aabcaabc22
样例输出
a:0,a:1,a:4,a:5 b:2,b:6 c:3,c:7 2:8,2:9
思路:
根据题目描述,我们可以使用双指针来遍历数组,使用指针p指向当前需要判断的字符,指针q遍历p之后的字符串判断是否存在与其相同的字符,若相同则输出字符与其对应位置。不过这样会存在一个问题:以上述样例为例,当p指向第二个a时,会输出与第一轮判断重复的结果。因此,我们需要一个标记来解决这个问题:输入的串中只包含数字和字母,所以可用‘*’来标记已判断过的重复字符。在遍历到‘*’跳过即可。
代码实现:
#include<stdio.h>
#define N 200
void main()
{
int i,flag;
char s[N],c;
while (scanf("%s", s) != EOF)
{
for (i = 0; s[i] != 0; i++)
{
flag = 0;
char* p = s+i;
char* q = p+1;
while (*p!='*'&& * q != 0)//遍历查找p之后的串中是否有与其相同的字符
{
if (*p == *q)//找到相同字符
{
if (flag == 0)//flag=0,当前字符为首个重复字符
{
printf("%c:%d", *p, p - s);//此时需将p输出
printf(",%c:%d", *p, q - s);
flag = 1;
*q = '*';//将重复字符标记为'*'
}
else if (flag == 1)flag=1,当前字符非首个重复字符
{
printf(",%c:%d", *p, q - s);//只需输出q
*q = '*';
}
}
else
q++;
}
if (flag == 1)//存在与p相同字符,将p标记
{
*p = '*';
printf("\n");
}
}
}
}
Problem H of 1725 :密码锁问题
题目描述
一个微调密码锁是这样的一种锁,这种锁你仅能转动密码盘。这是一种常见的密码盘,通过仅在允许的组中改变这些密码盘以微调某个值。
设想一行有D个编号的密码盘,每个密码盘顺序有0到9共九个数字。这类似于密码箱的组合锁。
下面是一系列B按钮,每个按钮标记有D位数字。例如,D可能是4标记就是1000 1200 1002 0111.按标记为1000的按钮,则仅转动第一个转盘一次,而其他转盘不动,而按按钮1002则转第一个转盘一次,转第四个转盘两次,剩下的不动。每个盘按循环的方式转动,即如果转到9,再转一次,就又转回0.
你的任务是模仿这样一个上锁的微调密码锁,给出最终的各密码盘的读数。输入
输入的每个测试数据的第一行包含有D个数字(至多10个),表示密码盘的起始位置。接下来的每一行有一排有标记的按钮,表示下一次会按的按钮。
输出
对每个测试用例用一行输出最终各密码盘的读数。
样例输入
0001 1003 0206 0034 1111 1003
样例输出
3348
思路:
首先,我们可以使用两个一维数组pssw和turn来分别存储密码的起始位置以和转动次数标记,然后定义一个函数实现转动密码盘功能。每次读入标记就调用该函数实现密码微调。再将数组turn清空,重复上一步操作。
有了以上思路就可以试着实现函数的功能,在此可以定义两个形参psw和t来接收主函数传递过来的起始密码和标记。接下来在函数体内循环遍历每一个密码盘,对每个盘完成对应的操作。这里我们用循环加1的方法来模拟转动密码盘,转动一次则使*psw+1。当密码盘数字为9时再转动一次回到起点,此时需要将*psw重置为0。
函数实现:
void fun(char* psw, char* t)
{
while (*t != 0)
{
int n = *t - 48;//记下当前密码盘需要转动的次数
for (int i = 0; i < n; i++)
{
if (*psw-48 + 1 >= 10)//转回起点
*psw = '0';
else
*psw += 1;//加1模拟转动一次
}
t++; psw++;//定位到下一个密码盘
}
}
完整代码如下:
#include<stdio.h>
#include<string.h>
#define N 15
void fun(char* psw, char* t)
{
while (*t != 0)
{
int n = *t - 48;//记下当前密码盘需要转动的次数
for (int i = 0; i < n; i++)
{
if (*psw-48 + 1 >= 10)//转回起点
*psw = '0';
else
*psw += 1;//加1模拟转动一次
}
t++; psw++;//定位到下一个密码盘
}
}
void main()
{
char pssw[N],turn[N];
scanf("%s", pssw);
while (scanf("%s", turn) != EOF)//读入标记
{
fun(pssw, turn);//转动密码盘
}
printf("%s\n", pssw);
}
Problem J of 1725:字符串操作
题目描述
输入一长度为n的字符串,若其n为偶数,则将字符串从中间反转,若为奇数,则将前后各(n-1)/2个字符反转,中间字符不动。
输入
输入有多组,每行有一个字符串
输出
对应每行输出反转后的字符串
样例输入
asdfghjkl qwerty
样例输出
fdsaglkjh ewqytr
思路:
我们首先忽略奇偶问题,对于反转字符,可以使用双指针解决。
定义指针p和q。让p指向字符串首,q指向字符串末尾,交换p与q所指字符。完成一次操作后p向后移动q向前移动,继续交换,直至指针p与指针q碰撞即可完成字符反转。
接下来用函数实现上述功能:
void rvs(char a[N],int n){
char* p, * q, t;
p = a;
q = a + n - 1;
while (p < q)//若指针p与指针q碰撞则结束循环,完成字符反转
{
t = *p;
*p = *q;
*q = t;
p++;//向后移动
q--;//向前移动
}
}
有了函数实现反转,接下来考虑长度奇偶问题(题目将字符串分为两个部分,分别反转左半边和右半边)。要解决这个问题很简单,首先需要找出奇偶之间对应操作的不同。
当n为偶数时,右半边起点为n/2。当n为奇数时,右半边起点为n/2+1。
完整代码如下:
#include<stdio.h>
#include<string.h>
#define N 200
void rvs(char a[N],int n){
char* p, * q, t;
p = a;
q = a + n;
while (p < q)//若指针p与指针q碰撞则结束循环,完成字符反转
{
t = *p;
*p = *q;
*q = t;
p++;//向后移动
q--;//向前移动
}
}
void main()
{
char s[N];
while (scanf("%s",s) != EOF)
{
int n = strlen(s);//求串长
if (n % 2 == 0)//偶数情况
{
rvs(s,n/2-1);//左半边奇偶一致
rvs(s+n/2,n-1);
}
else//奇数情况
{
rvs(s,n/2-1);
rvs(s+n/2+1,n-1);
}
printf("%s\n", s);
}
}
okk。