HDU 5525:Product 欧拉定理

原创 2015年11月17日 20:34:12

Product

 
 Accepts: 21
 
 Submissions: 171
 Time Limit: 6000/3000 MS (Java/Others)
 
 Memory Limit: 131072/131072 K (Java/Others)
问题描述
给n个数{A}_{1},{A}_{2}....{A}_{n}A1,A2....An,表示N=\prod_{i=1}^{n}{i}^{{A}_{i}}N=i=1niAi。求N所有约数之积。
输入描述
输入有多组数据.
每组数据第一行包含一个整数n.(1\leq n\leq {10}^{5})(1n105)
第二行n个整数{A}_{1},{A}_{2}....{A}_{n}A1,A2....An,保证不全为0.(0\leq {A}_{i}\leq {10}^{5})(0Ai105).
数据保证 \sum n\leq 500000n500000.
输出描述
对于每组数据输出一行为答案对{10}^{9}+7109+7取模的值.
输入样例
4
0 1 1 0
5
1 2 3 4 5
输出样例
36
473272463

官方题解:


做了一道div1的第三题,很好的一道题

参考了别人的代码,其中一处是(p^index)%mod,如果index很大的话,根据欧拉定理,就可以变成(p^(index%(phi(mod))))%mod。记录一下。

代码:

#pragma warning(disable:4996)  
#include <iostream>  
#include <algorithm>  
#include <cmath>  
#include <vector>  
#include <string>  
#include <cstring>  
using namespace std;

typedef long long ll;

ll N;
const int maxn = 100010;
const ll mod = 1e9 + 7;
const ll mod2 = 2 * (mod - 1);
ll index[maxn], L[maxn], R[maxn];
int num, pri[maxn], vis[maxn];
vector<int>have[maxn];

void init()
{
	num = 0;
	int i, j, n;
	for (i = 2; i < maxn; i++)
	{
		if (vis[i])
			continue;
		pri[++num] = i;
		for (j = i + i; j < maxn; j = j + i)
		{
			vis[j] = 1;
		}
	}
	for (i = 1; i < maxn; i++)
	{
		n = i;
		for (j = 1; j <= num&&pri[j] <= n; j++)
		{
			while (n%pri[j] == 0)
			{
				have[i].push_back(j);//记录所含有的质数,用质数的下标记录
				n /= pri[j];
			}
		}
	}
}

ll getresult(ll A, ll n, ll k)
{
	ll b = 1;
	while (n > 0)
	{
		if (n & 1)
		{
			b = (b*A) % k;
		}
		n = n >> 1;
		A = (A*A) % k;
	}
	return b;
}

void solve(int cishu, int n)
{
	int temp;
	int si = have[cishu].size();
	for (int i = 0; i < si; i++)
	{
		temp = have[cishu][i];
		index[temp] = (index[temp] + n) % mod2;
	}
}

int main()
{
	//freopen("i.txt", "r", stdin);
	//freopen("o.txt", "w", stdout);

	int x;
	ll k, n, ans;
	init();
	while (scanf("%d", &N) == 1)
	{
		memset(index, 0, sizeof(index));
		for (int i = 1; i <= N; i++)
		{
			scanf("%d", &x);
			solve(i, x);
		}
		int p = 1;
		while (pri[p] < N)
			p++;
		N = p;
		L[0] = R[N + 1] = 1;
		for (int i = 1; i <= N; i++)
			L[i] = L[i - 1] * (index[i] + 1) % mod2;
		for (int i = N; i >= 1; i--)
			R[i] = R[i + 1] * (index[i] + 1) % mod2;
		ans = 1;
		for (int i = 1; i <= N; i++)
		{
			k = L[i - 1] * R[i + 1] % mod2;
			n = index[i] * (index[i] + 1) / 2 % mod2;
			ans = ans*getresult(pri[i], n*k%mod2,mod) % mod;
		}
		printf("%lld\n", ans);
	}
	//system("pause");
	return 0;
}




版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

hdu5525 Product 费马小定理

Product    Accepts: 21    Submissions: 171  Time Limit: 6000/3000 MS (Java/Others)    Memory Limit:...

HDU 5525(Product-费马小定理)

给你一个数NN,求其所有约数的积。令N=pa11∗pa22∗...∗pannN=p_1^{a_1}*p_2^{a_2}*...*p_n^{a_n}于是答案=∏ni=1pi(ai+1)∗ai2∗∏nj=...

hdu 5525 Product(数论)

题目链接:hdu 5525 Product解题思路首先要把每个数拆分成质因子考虑,对于质因子p,假设出现的次数为c,那么对于包含p0,p1,p2...pcp^0,p^1,p^2...p^c的因子个数是...

HDU 5525 Product

Problem Description Given a number sequence A1,A2....An,indicating N=∏ni=1iAi.What is the prod...

hdu5525Product+约数的乘积

Product Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others) Tota...

欧拉定理——PowMod ( HDU 5728 )(2016 Multi-University Training Contest 1 1006)

题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5728 分析: 题中有 k = ∑(m,i=1) φ(i∗n) mod 1000000007,即...
  • FeBr2
  • FeBr2
  • 2016-07-22 20:01
  • 175

HDU 1395(欧拉定理)

欧拉φ函数的值  通式:φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数,x是不为0的整数。φ(1)=1(唯...

hdu 4549 M斐波那契数列 【矩阵+快速幂+欧拉定理】

M斐波那契数列Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others) Total Su...

hdu 3221 欧拉定理 +快速幂+矩阵快速幂

首先这道题考的是递归调用的计数问题 最终返回时分别是1时执行a次和2时执行b次,根据计数原理: 根据递归的调用机制和排列组合的乘法原则,得出f[n]=f[n-1]*f[n-2]; 而递归的增长随...

Hdu 1395 2^x mod n = 1 (欧拉定理 分解素因数)

先吐槽一下:考完试后第一天实验室就断网,简直不可思议。 题意:给出n,找出最小的x使其满足2^x mod n = 1 思路:欧拉定理(费马小定理的欧拉推广)可得如果n为1或2的倍数则无解,n为大于...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)