hdu 3524 Perfect Squares【打表、除法取余、快速幂】

先打表，找规律

#include<iostream> #include<set> using namespace std; int main() { set<long long>myset; for(long long i=1;i<20;i++,cout<<endl) { myset.clear(); for(long long j=0;j<1000000;j++) myset.insert(j*j%(1<<i)); cout<<i<<' '<<myset.size(); } return 0; }

f(i)=4*f(i-2)-5 通项 f(i)=(2^(n-1)+5)/3 i为奇数
f(i)=4*f(i-2)-4 f(i)=(2^(n-1)+4)/3 i为偶数
即 f(i)=(2^(n-1)+4+n%1)/3
若b%c=1，则a/b%c=a%c=d，证明如下：
b=c*k1+1 a=b*(c*k2+d)=(c*k1+1)*(c*k2+d)=c*(k1*k2+k1*d+k2)+d=c*k3+d
因此a%c=d=a/b%c
求解f(i)%mod，a/b%c=(a*k)/(b*k)%c=a*k%c，（b*k%c=1）
即求任一k使b＊k%c=1, 即b*k=c*k'+1，设x=k，y=-k’，则x*b+y*c=1，extgcd或者打表求一个x即可

#include<iostream> using namespace std; int extgcd(int &x,int &y,int a,int b) { if(b==0){x=1,y=0;return a;} int d=extgcd(x,y,b,a%b); int t=x;x=y;y=t-a/b*y; return d; } int main() { int x,y; extgcd(x,y,3,10007); cout<<x<<endl; }

#include<iostream> using namespace std; int main() { for(int i=1;i<1000000;i++) if(i*3%10007==1) { cout<<i<<endl;break; } }

所以f(i)=(2^(n-1)+4+n%1)*3336%10007.

#include<iostream> using namespace std; typedef long long ll; const ll mod=10007; ll pow(ll n) { ll ans=1,base=2; while(n) { if(n&1) ans=ans*base%mod; base=base*base%mod,n>>=1; } return ans; } int main() { ll t,n,cas=1; cin>>t; while(t--) { cin>>n; cout<<"Case #"<<cas++<<": "<<3336*(pow(n-1)+4+(n&1))%mod<<endl; } return 0; }