2023年第7周学习总结

素数筛

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。

  • 3
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值