[bzoj1223] [HNOI2002]Kathy函数

　　首先由题解可得TAT，f(i)=i当且仅当i在二进制下为回文串。

　　那么问题就变成了1~n中有多少个二进制下的回文串。

　　把m转成2进制后就是正常的统计了= =。

　　f[i]表示二进制下，有多少个i位的回文串（包括前导0）。f[1]=f[2]=2,f[i]=f[i-2]*2.....也就是f[i]=2^((i+1)/2)

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #include<algorithm>
 #define ri register int
 using namespace std;
 const int modd=;
 char s[];
 int two[][];
 int a[],l[],len,len2,mp[],ten[];
 int i,j,k,n,m;
 int ans[],ansl;

 inline bool notsmall(int x){
     if(l[x]!=len)return l[x]>len;
     for(ri i=len;i;i--)if(two[x][i]!=a[i])return two[x][i]>a[i];
     return ;
 }
 inline bool notbig(int x){
     if(l[x]!=len)return l[x]<len;
     for(ri i=len;i;i--)if(two[x][i]!=a[i])return two[x][i]<a[i];
     return ;
 }
 inline void add(int x){
     l[x]=l[x-];int y=x-;
     for(ri i=;i<=l[x];i++){
         two[x][i]+=two[y][i]<<;
         if(two[x][i]>=modd)two[x][i]-=modd,two[x][i+]++;
     }
     if(two[x][l[x]+])l[x]++;
 }
 inline void addtoans(int x){
     for(ri i=;i<=l[x];i++){
         ans[i]+=two[x][i];
         if(ans[i]>=modd)ans[i+]++,ans[i]-=modd;
     }
     while(ans[ansl+])ansl++;
 }
 inline void dec(int x){
     for(ri i=;i<=l[x];i++){
         a[i]-=two[x][i];
         if(a[i]<)a[i]+=modd,a[i+]--;
     }
     while(!a[len]&&len>)len--;
 }
 inline void run(){
     ri i,j,k,len;
     two[][l[]=]=;
     for(i=;!notsmall(i);)i++,add(i);
     for(;i;i--)
         if(notbig(i))dec(i),mp[i]=,len2=max(len2,i);
 //    for(i=len2;i;i--)printf(" %d",mp[i]);puts("");
     for(i=;i<len2;i++)addtoans((i+)>>);
     int flag=;
     for(i=len2>>;i;i--)swap(mp[i],mp[len2-i+]);
     for(i=;i<=len2&&flag;i++){
         if(mp[i]){
             if(i<=(len2>>))addtoans((len2-(i<<)+)>>);
             else{
                 if(mp[len2-i+]==||len2-i+==i)addtoans();
             }
         }
         if(i>(len2>>))flag&=(mp[len2-i+]==mp[i]);
     }
     if(flag==||len2==)addtoans();
 }
 int main(){
     scanf("%s",s);j=strlen(s);
     for(i=ten[]=;i<=;i++)ten[i]=ten[i-]*;
     for(i=j-;i>=;i--){
         if(!k)len++;
         a[len]+=(s[i]-'')*ten[k];
         k++;if(k>)k=;
     }
     run();
     for(printf("%d",ans[ansl]),i=ansl-;i>;i--){
         for(j=;j<modd;j*=)if(ans[i]<j)putchar('');
         printf("%d",ans[i]);
     }puts("");
     return ;
 }

代码又丑又慢>_<