【题解】sdoj2597[2018.8.10集训]A.细胞分裂 费马小定理降幂+快速幂取模

题目链接
问题描述

小 A 养了一大坨细胞。
最初小 A 只有 1 个细胞。 每秒， 小 A 的每个细胞都会分裂成 2 个细胞。
已知： 现在离“最初”已经过去了 x 秒， 那么现在的细胞数当然是可以计算的。
小 A 想知道的当然不是当前的细胞数。 小 A 知道他养的细胞的习性： 每 y 个细胞
会聚成一团。 经常会有剩下的细胞， 那么我们称这些细胞是孤独的。
小 A 想知道的就是孤独的细胞个数

输入

输入文件共一行， 为两个整数 x y， 以空格隔开

输出

输出文件共一行， 为一个整数， 即孤独的细胞个数。

输入样例

3 3

输出样例

2

数据规模和约定

对于 10%的数据， x<2^6。
对于 20%的数据， x<2^17。
对于 40%的数据， x<2^64。
对于 70%的数据， x<2^2333。
对于 100%的数据， 0≤x<2^233333， y 是 3 到 1000 之间（含两端） 的质数

---

题解

比赛时看到这题毫无头绪（我太弱了），首先打了个线性筛把素数表找到，然后开始打表找规律，发现对于每个素数，答案是循环的。于是把每个素数的循环长度和每一位的值存下来，再对x模上每一个素数的循环长度（我不知道大概有多少位所以都只读取每一位的数），本来想着能拿点分就好，结果AC了。正解是费马小定理降幂，对于指数x模上(p-1)，可以把x降得较小，然后快速幂取模。
正解

#include<cstdio>
#include<cstring>
const int N=1e5+10;
char num[N];
int qpow(int a,int x,int p)
{
	int ret=1;
	while(x)
	{
		if(x&1)ret=(ret*a)%p;
		a=(a*a)%p;
		x>>=1;
	}
	return ret%p;
}
int main()
{
	int p;
	scanf("%s%d",num,&p);
	int nowx=0;
	for(int i=0;i<strlen(num);i++)
	nowx=(nowx*10+num[i]-'0')%(p-1);
	printf("%d\n",qpow(2,nowx,p));
	return 0;
}

比赛AC代码

#include<cstdio>
#include<vector>
#include<cstring>
#define PB(v) push_back(v)
using namespace std;
typedef long long ll;
int prime[1010],p;
bool iscomp[1010];
vector<int>v[1010];
char x;
int xh[1010];
void primetable()
{
	for(int i=2;i<=1000;i++)
	{
		if(!iscomp[i])prime[p++]=i;
		for(int j=0;j<p&&i*prime[j]<=1000;j++)
		{
			iscomp[i*prime[j]]=1;
			if(i%prime[j]==0)break;			
		}
	}
}
int main()
{
	//freopen("cell.in","r",stdin);
	//freopen("cell.out","w",stdout);
	memset(xh,0,sizeof(xh));
	memset(iscomp,0,sizeof(iscomp));
	primetable();
	int y;
	for(int i=1;i<p;i++)
	{
		int base=1;
		for(int j=1;j<=1000;j++)
		{
			base*=2;
			base%=prime[i];
			if(v[prime[i]].size()>=1&&base==v[prime[i]][0])break;
			v[prime[i]].PB(base);
		}
	}
	while(scanf("%c",&x)&&x!=' ')
	{
		for(int i=1;i<p;i++)
		{
			xh[prime[i]]=xh[prime[i]]*10+x-'0';
			xh[prime[i]]%=v[prime[i]].size();
		}
	}
	scanf("%d",&y);
	if(xh[y]==0)printf("%d\n",v[y][v[y].size()-1]);
	else printf("%d\n",v[y][xh[y]-1]);
	return 0;
}