[bzoj]1041: [HAOI2008]圆上的整点

原题链接：圆上的整点

看到下面提示里的视频就点进去看了。。直接套结论（大雾

提示里的视频不完整，完整视频链接：【官方双语】隐藏在素数规律中的π

引入一个x(kai)函数（周期函数）1,0,-1,0,1,0,-1,0……

总之x(2n)=0，x(4n+1)=1,x(4n+3)=-1。

将r^2分解质因数，然后圆上的整点个数就是

***

简单来说就是把r分解质因数O(logn)，然后对于每一个质因数，把1到p（p为该质因数的个数）的x值加起来O(logn)，最后累乘再成4得到。。代码实现就很简单了。。

#include <iostream>
#include <cstdio>
using namespace std;
struct kk{
    int num;int cal;
};
int r,qwq=0,ans=1,tmp;
kk a[1000000];
int ksm(int k,int p);
int kai(int qaq);
int main()
{
    int j=1;
    scanf("%d",&r);
    while(r>1){
        j++;
        if(r%j==0)a[++qwq].num=j;
        while(r%j==0){
            r/=j;
            a[qwq].cal++;
        }
    }
    /*for(int i=1;i<=qwq;i++)
        cout<<a[i].num<<"  "<<a[i].cal<<endl;*/
    for(int i=1;i<=qwq;i++){
        a[i].cal*=2;
        tmp=0;
        for(int po=0;po<=a[i].cal;po++){
            tmp+=kai(ksm(a[i].num,po));
            //cout<<tmp<<' '<<po<<endl;;
        }
        ans*=tmp;
    //cout<<endl;
    }
    ans*=4;
    cout<<ans<<endl;
    return 0;
}
int kai(int qaq){
    if(!(qaq&1))return 0;
    if(!((qaq-1)%4))return 1;
    if(!((qaq-3)%4))return -1;
}
int ksm(int k,int p){
    int ans=1;
    while(p){
        if(p&1)ans*=k;
        k*=k;
        p>>=1;
    }
    return ans;
}

总复杂度大概就是O(logn)（瞎算，逃）