Hdu-1215 七夕节【算术基本定理应用】

news/2024/7/10 5:23:31 标签: 优化

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1215

解题思路:

这道题要求的就是一个整数的因子之和。

但是这道题在杭电数据比较弱,各种方法都能过。暴力都能过,但是在我们OJ,只能使用优化的代码。因为数据是1-500000,暴力的话,估计得跑个几十秒。

下面说下这道题的优化解法:

1.素数打表,这个就不说了,基本功。

2.接下来就是这个问题的关键:素分解。不会的可以移步这篇文章:http://wenku.baidu.com/view/e55ca209ba1aa8114431d98a.html


代码如下:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<cmath>
#include<algorithm>
using namespace std;

#define CLR(arr, what) memset(arr, what, sizeof(arr))
const int N = 500002;

int prim[N];
bool visit[N];

void fast_prim() //素数打表
{
	memset(visit, true, sizeof(visit));
	int limit = (int)sqrt(N * 1.0);
	for(int i = 2; i < limit; ++i)
		if(visit[i])
			for(int j = i * i; j < N; j += i)
				visit[j] = false;
	for(int i = 2, j = 0; i < 10000; ++i)
		if(visit[i])
			prim[j++] = i;
}

int prim_reduce(int n) //整数素分解
{
	int limit, num, sum, total = 1, temp = n;
	limit = sqrt(N * 1.0);
	for(int i = 0; prim[i] * prim[i] <= n; ++i)
	{
		num = sum = 1;
		if(n == 1)
			break;
		while(n % prim[i] == 0)
		{
			num *= prim[i];
			n /= prim[i];
			sum += num;
		}
			total *= sum;
	}
	if(n != 1)
		total *= (n + 1);
	return total - temp;
}

int main()
{
	fast_prim();
	int ncase;
	int num;
	scanf("%d", &ncase);
	while(ncase--)
	{
		scanf("%d", &num);
		if(num == 1)
		{
			printf("0\n");
			continue;
		}
		else if(visit[num])
		{
			printf("1\n");
			continue;
		}
		else
			printf("%d\n", prim_reduce(num));
	}
	return 0;
}

第二种方法:

在于打表的过程不同,第一种方法是解决这种问题的通法,但这道题我们可以用另一种打表方面,a[i]代表的就是这个因子之和,这样,我们就可以一次打表,O(1)输出。

代码如下:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<cmath>
#include<algorithm>
using namespace std;

#define CLR(arr, what) memset(arr, what, sizeof(arr))
const int N = 500002;
const int M = 250000;
int a[N] = {0, 0};

void fash()
{
	for(int i = 2; i < N; ++i) //因子初始化为1
		a[i] = 1;
	for(int i = 2; i <= M; ++i) //一半即可
		for(int j = 2 * i; j < N; j += i)
			a[j] += i;
}

int main()
{
	fash();
	int ncase;
	scanf("%d", &ncase);
	while(ncase--)
	{
		int num;
		scanf("%d", &num);
		printf("%d\n", a[num]);
	}
	return 0;
}



http://www.niftyadmin.cn/n/1516679.html

相关文章

NYOJ-453 小珂的烦恼【模拟||找规律】

题目链接&#xff1a;http://acm.nyist.net/JudgeOnline/problem.php?pid453 解题思路&#xff1a; 简单的可以模拟&#xff0c;用一个数组来维护数组中未使用过的最小数字&#xff0c;若使用&#xff0c;就把它对应的另一个数也标记为使用过。模拟一下即可。 代码如下&…

len是python的内置函数吗_Python3.6中是否可以模仿内置len()函数?

在Python3.6中&#xff0c;可以模拟内置的len()函数吗&#xff1f;在我有一个类&#xff0c;它定义了一个依赖于len()函数的简单方法&#xff0c;如下所示&#xff1a;class MyLenFunc(object):def is_longer_than_three_characters(self, some_string):return len(some_string…

C++this指针的理解

先要理解class的意思。class应该理解为一种类型&#xff0c;象int,char一样&#xff0c;是用户自定义的类型。&#xff08;虽然比int char这样build-in类型复杂的多&#xff0c;但首先要理解它们一样是类型)。用这个类型可以来声明一个变量&#xff0c;比如int x, myclass my等…

NYOJ-421 打扑克【模拟】

题目链接&#xff1a;http://acm.nyist.net/JudgeOnline/problem.php?pid421 解题思路&#xff1a; 这道题的原型大家貌似都知道&#xff0c;各地都有不同的叫法。&#xff08;汗&#xff0c;我竟然不知道&#xff5e;&#xff09;题意很好理解&#xff0c;本来想的是暴力的…

如何理解python面向对象_如何理解python面向对象编程

在本篇文章里小编给大家整理的是关于如何理解python面向对象编程的相关知识点&#xff0c;需要的朋友们可以学习下。类是面向对象程序设计的一部分。面向对象程序设计或者简称为 OOP 致力于创建可重用代码块称之为类。当你想在你的程序中使用类时&#xff0c;你会从类中创建一个…

NYOJ-485 A*B Problem【九余数定理】

题目链接&#xff1a;http://acm.nyist.net/JudgeOnline/problem.php?pid485 解题思路&#xff1a; 经典的九余数定理&#xff0c;可以参考这篇文章&#xff1a;http://blog.csdn.net/niushuai666/article/details/6856109 代码如下&#xff1a; #include<iostream> #…

WinCE BSP工程文件pbcxml分析

一.pbcxml分析 每一个BSP都有一个工程文件&#xff0c;里面描述了BSP的信息。 下面就来介绍一下BSP的pbcxml文件。 文件的大致格式应该是这样的&#xff1a; <CatalogFile …> <FileInformation ….> …. </FileInformation> <BSP …> … <…

NYOJ-452 ShippingCubes【数学】

题目链接&#xff1a;http://acm.nyist.net/JudgeOnline/problem.php?pid452 解题思路&#xff1a; 这道题是topcoder的一次比赛题目&#xff0c;但是数据太水&#xff0c;暴力可过。3层循环无压力。 但是如果是正规做法&#xff0c;应该是先将n开3次方&#xff0c;因为要3…