CF Gym100917 C

要找到和为给定值的所有的等比数列.

1肯定是要特判一下.

我的想法是先找到所有等比为1的,等比为1就是将这个数分为相同的一些数,总共就是这个数的所有约数个数减一(有一个约数为1,题目要求至少分成两个数).

对于其他的等比不为1 的,用等比数列的求和公式暴力一发就行了.

#include<cstdio>
#include<cstring>
#include<cmath>

using namespace std;

typedef long long ll;
const int MAXN=1e5+5;
int prime[MAXN];
int check[MAXN];
int tot=0;
ll n;
ll ans;


void creat_prime()
{
	for(int i=2;i<MAXN;i++)
	{
		if(!check[i]) prime[tot++]=i;
		for(int j=0;j<tot && prime[j]*i<MAXN;j++)
		{
			check[i*prime[j]]=true;
			if(i%prime[j]==0) break;
		}
	}
}

void init()
{
	ans=1;
	ll m=n;
	for(int i=0;i<tot;i++)
	{
		int tmp=0;
		while(m%prime[i]==0)
		{
			tmp++;
			m/=prime[i];
		} 
		ans*=(tmp+1);
		if(m==1) break;
	}
	ans--;
}

ll quick_pow(ll a,ll b)
{
	ll ret=1;
	while(b)
	{
		if(b&1) ret*=a;
		b>>=1; a*=a;
	}
	return ret;
}

int main()
{
	//printf("%lld",quick_pow(2,10));
	creat_prime();
	while(~scanf("%I64d",&n))
	{
		if(n==1)
		{
			ans=0;
			printf("%I64d\n",ans);
			continue;
		}
		init();
		for(int i=2;i<n;i++)
		{
			for(int j=2;j<n;j++)
			{
				ll tmp=(quick_pow(i,j)-1)/(i-1);
				if(tmp>n) break;
				if(n%tmp==0)
				{
					ans++;
				}
			}
		}
		printf("%I64d\n",ans);
	} 
	
}