NOIP模拟题——细胞分裂

【问题描述】
小 A 养了一大坨细胞。 最初小 A 只有 1 个细胞。每秒，小 A 的每个细胞都会分裂成 2 个细胞。

已知：现在离“最初”已经过去了x秒，那么现在的细胞数当然是可以计算的。 小 A 想知道的当然不是当前的细胞数。小 A 知道他养的细胞的习性：每 y 个细胞会聚成一团。经常会有剩下的细胞，那么我们称这些细胞是孤独的。 小 A 想知道的就是孤独的细胞个数。
【输入文件】
输入文件为 cell.in。 输入文件共一行，为两个整数 xy，以空格隔开。
【输出文件】
输出文件为 cell.out。 输出文件共一行，为一个整数，即孤独的细胞个数。
【输入样例】

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 之间（含两端）的质数。

法1(70分)：高精度+快速幂

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cstring>
 4 using namespace std;
 5 int x[300000],y[300000];
 6 int qwe;
 7 long long solve(int *x)
 8 {
 9     if((x[0]==0||x[0]==1)&&x[1]==1)return 2;
10     int pre=0;int temp=0;
11     bool pd=true;
12     if(x[1]/2!=0)temp++;
13     if(x[x[0]]%2==0)pd=false;//偶数
14     for(int i=1;i<=x[0];i++)
15     {
16         int k=pre*10+x[i];
17         y[temp++]=k/2;
18         pre=k%2;
19     }
20     if(x[1]/2!=0)y[0]=x[0];else y[0]=x[0]-1;
21     for(int i=1;i<=y[0];i++)
22     x[i]=y[i];
23     x[0]=y[0];
24     if(pd==true)//奇数
25     {
26         long long b=solve(x)%qwe;
27         b=b*b*2%qwe;
28         return b;
29     }
30     else
31     {
32         long long b=solve(x)%qwe;
33         b=b*b%qwe;
34         return b;
35     }
36 }
37 int main()
38 {
39     freopen("cell.in","r",stdin);
40     freopen("cell.out","w",stdout);
41     int temp=0;
42     char c=getchar();
43     while(c<‘0‘||c>‘9‘)c=getchar();
44     while(c>=‘0‘&&c<=‘9‘)
45     {
46         x[++temp]=c-‘0‘;
47         c=getchar();
48     }
49     x[0]=temp;
50     scanf("%d",&qwe);
51     long long ans=solve(x);
52     printf("%I64d",ans%qwe);
53     return 0;
54 }

法2：费马小定理：

假如p是质数，且gcd(a,p)=1，那么 a^(p-1)≡1（mod p）。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;
char s[100005];int x,y;
void solve()
{
    char c;
    int len=strlen(s);
    for(int i=0;i<len;i++)
    {
        c=s[i];
        x=((c-‘0‘)+x*10)%(y-1);
    }
}

int main()
{
    freopen("cell.in","r",stdin);
    freopen("cell.out","w",stdout);
    scanf("%s %d",s,&y);
    solve();
    int ans=1;
    for(int i=1;i<=x;i++)
    ans=(ans*2)%y;
    printf("%d",ans);
}