等比数列二分求和取模

最新推荐文章于 2021-02-26 01:31:21 发布

BRCOCOLI

最新推荐文章于 2021-02-26 01:31:21 发布

阅读量936

点赞数 1

分类专栏： POJ+百练数论

本文链接：https://blog.csdn.net/qq_34446253/article/details/52126329

版权

POJ+百练同时被 2 个专栏收录

120 篇文章 2 订阅

订阅专栏

数论

9 篇文章 0 订阅

订阅专栏

Sn= a+a2+...+an

要求 Sn mod p
如果用公式算，可能溢出，因此用二分法求
1) 若 n是偶数
Sn= a+...+ a^(n/2) + a^(n/2+1) + a^(n/2+2) +...+ a^(n/2+n/2)

=(a+...+ a^(n/2)) + a^(n/2)*(a+...+ a^(n/2))
=Sn/2+ a^(n/2)Sn/2
=(1+ a^(n/2))Sn/2

2) 若n是奇数，那么n-1是偶数，只要用上面推出来的公式S(n-1)再加上a^n就好

Sn=( 1+a^((n-1)/2) )*S(n-1)/2+a^n

#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long LL;
int PowMod(int a,int n,int m)  //a^n 结果对m取模 
{
	if(n==0) return 1;
	int x = PowMod(a,n/2,m); 
	long long ans = (long long) x*x%m;
	if(n%2==1) ans=ans*a%m; //但次数为奇数时还要多乘一次 
	return (int)ans;
}
int PowSumMod(int a,int n,int m)
{
	if(n==1) return a%m;
	if(n%2==0) return (1+PowMod(a,n/2,m))*PowSumMod(a,n/2,m)%m;
	else return ((1+PowMod(a,(n-1)/2,m))*PowSumMod(a,(n-1)/2,m)%m + PowMod(a,n,m))%m;
}