bzoj 1223: [HNOI2002]Kathy函数 数位DP 高精度
1223: [HNOI2002]Kathy函数
Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 207 Solved: 90
[Submit][Status]
Description
Input
仅有一行，为正整数m
Output
输出仅有一个正整数，表示所有的满足f(n)=n,(n<=m) 的自然数的个数。
Sample Input
Sample Output
这道题的高精度模板相对又有完善，但还存在bugs，使用时尽量讲位数多开几倍。
这道题由于递归方程转移有2、4，所以应该往二进制那边去想，可以通过归纳法证明f(n)=n当且仅当n二进制为回文串（相当神奇），然后就是数位DP，最开始我将数据范围想成了m<=2^100，还说刚好long long可以存，但是10^100就必须写高精了。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<unistd.h>
using namespace std;
typedef long long qword;
#define MAXN 20100
class number//四位
{
public:
const static int max_val=;
const static int max_len=;
number()
{
clear();
}
number(int x)
{
*this=x;
}
bool is_odd()const
{
return numb[]%==;
}
bool is_even()const
{
return numb[]%==;
}
int size()const
{
return topn;
}
int length()const
{
int x=numb[topn];
int ret=;
while (x)
{
ret++;
x/=;
}
int y=;
x=number::max_val;
while (x)
{
y++;
x/=;
}
y--;
ret+=topn*y;
return ret;
}
number pow(int x)const
{
number ret,a;
a=*this;
ret=;
while (x)
{
if (x&)ret=ret*a;
a=a*a;
x>>=;
}
return ret;
}
void operator =(int x)//{{{
{
int now=;
clear();
numb[now]=x;
while (numb[now]>=number::max_val)
{
numb[now+]+=numb[now]/number::max_val;
numb[now]%=number::max_val;
now++;
if (now>topn)topn=now;
}
}//}}}
void operator =(number num)
{
topn=num.topn;
memcpy((this->numb),num.numb,sizeof(num.numb[])*(topn+));
}
void operator +=(number num)//{{{
{
int i;
topn=max(topn,num.topn);
for (i=;i<=topn;i++)
{
numb[i]+=num.numb[i];;
if (numb[i]>=number::max_val)
{
numb[i+]+=numb[i]/number::max_val;
numb[i]%=number::max_val;
}
}
while (numb[topn+])
{
topn++;
numb[topn+]+=numb[topn]/number::max_val;
numb[topn]%=number::max_val;
}
}//}}}
void operator +=(int x)//{{{
{
int now=;
if (topn==-)topn=;
numb[now]+=x;
while (numb[now]>=number::max_val)
{
numb[now+]+=numb[now]/number::max_val;
numb[now]%=number::max_val;
now++;
if (now>topn)topn=now;
}
}//}}}
void operator *=(int x)//{{{
{
int i;
for (i=;i<=topn;i++)
{
numb[i]*=x;
}
for (i=;i<=topn;i++)
{
if (numb[i]>=number::max_val)
{
numb[i+]+=numb[i]/number::max_val;
numb[i]%=number::max_val;
}
}
while (numb[topn+])
{
topn++;
numb[topn+]+=numb[topn]/number::max_val;
numb[topn]%=number::max_val;
}
}//}}}
void operator -=(number &num)//{{{
{
if (*this<num)throw "Error!\n->void operator -=(number &num)\n";
int i;
for (i=;i<=topn;i++)
{
numb[i]-=num.numb[i];
}
for (i=;i<=topn;i++)
{
while (numb[i]<)
{
numb[i]+=number::max_val;
numb[i+]--;
}
}
while (topn&&!numb[topn])topn--;
}//}}}
void operator --(int)//{{{
{
if (topn==&&numb[]==)throw "Error!\n->void operator --(int)\n";
int now=;
numb[now]--;
while (numb[now]<)
{
numb[now+]--;
numb[now]+=number::max_val;
}
while (topn&&!numb[topn])topn--;
}//}}}
private:
int numb[max_len];
int topn;
void clear()
{
topn=;
memset(numb,,sizeof(numb)); }
friend bool operator <(number num1,number num2);
friend bool operator <=(number num1,number num2);
friend bool operator ==(number num1,number num2);
friend ostream& operator <<(ostream &out,number num);
friend istream& operator >>(istream &in,number &num);
friend number operator *(number num1,number num2);
friend number operator *(number num,int x);
friend number operator +(number num1,number num2);
//a=a+b远没有a+=b快
};
bool operator <(number num1,number num2)//{{{
{
if (num1.topn!=num2.topn)
{
return num1.topn<num2.topn;
}
int i;
for (i=num1.topn;i>=;i--)
{
if (num1.numb[i]!=num2.numb[i])
{
return num1.numb[i]<num2.numb[i];
}
}
return false;
}//}}}
bool operator <=(number num1,number num2)//{{{
{
if (num1.topn!=num2.topn)
{
return num1.topn<num2.topn;
}
int i;
for (i=num1.topn;i>=;i--)
{
if (num1.numb[i]!=num2.numb[i])
{
return num1.numb[i]<num2.numb[i];
}
}
return true;
}//}}}
bool operator ==(number num1,number num2)//{{{
{
if (num1.topn!=num2.topn)return false;
for (int i=;i<=num1.topn;i++)
{
if (num1.numb[i]!=num2.numb[i])return false;
}
return true;
}//}}}
ostream& operator <<(ostream &out,number num)//{{{
{
int i;
out<<num.numb[num.topn];
for (i=num.topn-;i>=;i--)
{
//压六位时
//if (num.numb[i]<100000)out<<"0";
//if (num.numb[i]<10000)out<<"0";
if (num.numb[i]<)out<<"";
if (num.numb[i]<)out<<"";
if (num.numb[i]<)out<<"";
out<<num.numb[i];
}
return out;
}//}}}
istream& operator >>(istream &in,number &num)//{{{
{
string str;
in>>str;
int i;
num.clear();
for (i=(int)str.length()-,num.topn=;i>=;i-=,num.topn++)
{
if (i-<(int)str.length())
{
num.numb[num.topn]=(str[i]-'')+*(str[i-]-'')+*(str[i-]-'')+*(str[i-]-'');
}else
{
if (i-<(int)str.length())num.numb[num.topn]+=*(str[i-]-'');
if (i-<(int)str.length())num.numb[num.topn]+=*(str[i-]-'');
if (i<(int)str.length())num.numb[num.topn]+=(str[i]-'');
}
}
num.topn--;
return in;
}//}}}
number operator *(number num,int x)//{{{
{
number ret;
ret=num;
ret*=x;
return ret;
}//}}}
number operator *(number num1,number num2)
{
number ret;
ret.topn=num1.topn+num2.topn;
for (int i=;i<=num1.topn;i++)
for (int j=;j<=num2.topn;j++)
{
ret.numb[i+j]+=num1.numb[i]*num2.numb[j];
ret.numb[i+j+]+=ret.numb[i+j]/number::max_val;
ret.numb[i+j]%=number::max_val;
}
for (int i=;i<=ret.topn;i++)
{
ret.numb[i+]+=ret.numb[i]/number::max_val;
ret.numb[i]%=number::max_val;
}
while (ret.numb[ret.topn+])
{
ret.topn++;
ret.numb[ret.topn+]+=ret.numb[ret.topn]/number::max_val;
ret.numb[ret.topn]%=number::max_val;
}
while (ret.topn && !ret.numb[ret.topn])
{
ret.topn--;
}
return ret;
}
number operator +(number num1,number num2)//{{{
{
number ret;
ret=num1;
ret+=num2;
return ret;
}//}}} char str[MAXN];
int num0[MAXN];
bool num[MAXN];
number rec[MAXN][][];
int n,m;
number search_m(int now,int fl,int bo)
{
int x,y;
x=now-;y=n-now;
if (x>y)
{
number ret;
ret=(int)(fl || (!fl && !bo));
return ret;
}
number &res=rec[now][fl][bo];
number aa;
if (aa<rec[now][fl][bo])return res;
res=;
if (!fl)
{
if (num[y])
{
//
if (num[x]==)
res+=search_m(now+,false,bo);//
else
res+=search_m(now+,false,true);// if (now>)//首位非0
{
if (num[x]==)
res+=search_m(now+,true,false);//
else
res+=search_m(now+,true,bo);//
}
}else if (!num[y])
{
if (now>)//特殊处理m==0情况
{
if (num[x]==)
res+=search_m(now+,false,false);//
else
res+=search_m(now+,false,bo);//
}
}
}else
{
res+=search_m(now+,true,true);//
res+=search_m(now+,true,true);//
}
return res;
}
int main()
{
//freopen("input.txt","r",stdin);
scanf("%s\n",str);
int i,j,x,y;
n=strlen(str);
x=y=-;
for (i=n-;i>=;i-=)
{
x++;
for (j=max(,i-);j<=i;j++)
num0[x]=num0[x]*+str[j]-'';
}
n=x+;
x=;
while (~n)
{
num[x++]=num0[]%;
for (i=n-;i>=;i--)
{
if (i)num0[i-]+=(num0[i]&)*;
num0[i]/=;
}
while (n>= && !num0[n-])n--;
}
n=x;
/*for (i=n-1;i>=0;i--)
{
printf("%d",num[i]);
}
printf("\n");*/
/*for (i=0;i<=n;i++)
for (j=0;j<2;j++)
for (k=0;k<2;k++)
rec[i][j][k]=-1;*/
number ans;
ans=search_m(,false,false);
number xx;
for (i=;i<n;i++)
{
xx=;
xx=xx.pow((i-)/);
ans+=xx;
}
cout<<ans<<endl;
return ;
}
bzoj 1223: [HNOI2002]Kathy函数 数位DP 高精度的更多相关文章
- 【SPOJ 2319】 BIGSEQ - Sequence （数位DP+高精度）
BIGSEQ - Sequence You are given the sequence of all K-digit binary numbers: 0, 1,..., 2K-1. You need ...
- BZOJ.4513.[SDOI2016]储能表(数位DP)
BZOJ 洛谷 切了一道简单的数位DP,终于有些没白做题的感觉了...(然而mjt更强没做过这类的题也切了orz) 看部分分,如果\(k=0\),就是求\(\sum_{i=0}^n\sum_{j=0} ...
- [bzoj 1026]windy数（数位DP）
题目:http://www.lydsy.com/JudgeOnline/problem.php?id=1026 分析: 简单的数位DP啦 f[i][j]表示数字有i位,最高位的数值为j的windy数总 ...
- bzoj 1026 [SCOI2009]windy数 数位dp
1026: [SCOI2009]windy数 Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeOnline ...
- bzoj 3131 [Sdoi2013]淘金（数位DP+优先队列）
Description 小Z在玩一个叫做<淘金者>的游戏.游戏的世界是一个二维坐标.X轴.Y轴坐标范围均为1..N.初始的时候,所有的整数坐标点上均有一块金子,共N*N块. 一阵风吹 ...
- BZOJ 4513: [Sdoi2016]储能表 [数位DP ！]
4513: [Sdoi2016]储能表 题意:求\[ \sum_{i=0}^{n-1}\sum_{j=0}^{m-1} max((i\oplus j)-k,0) \] 写出来好开心啊...虽然思路不完 ...
- BZOJ 4521 CQOI 2016 手机号码 数位DP
4521: [Cqoi2016]手机号码 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 539 Solved: 325[Submit][Status ...
- BZOJ 2425 [HAOI2010]计数：数位dp + 组合数
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2425 题意: 给你一个数字n,长度不超过50. 你可以将这个数字: (1)去掉若干个0 ( ...
- bzoj 3131 [Sdoi2013]淘金（数位dp）
题目描述 小Z在玩一个叫做<淘金者>的游戏.游戏的世界是一个二维坐标.X轴.Y轴坐标范围均为1..N.初始的时候,所有的整数坐标点上均有一块金子,共N*N块. 一阵风吹过,金子的位置发生了 ...
随机推荐
- Java通过代理server上网
完整代码 package com.proj.net; //导入编码的jar文件 import it.sauronsoftware.base64.Base64; import java.io.Buffe ...
- CCEditBox用法
1.以下是CCEditBox的相关函数和类型说明: /* 编辑框的一些函数 setText("字符串"); //设置文本 setFontColor(color); //设置文本颜色 ...
- PureMVC(JS版）源码解析（一）：观察者模式解析
假设一种情景,在程序开发中,我们需要在某些数据变化时,其他的类做出相应,例如在游戏中,升级一件装备,会触发玩家金币数量改变,背包数据改变和冷却队列数据改变等等.我们不可能设置setInte ...
- 前后端分离--mock
fekit mock 数据 > fekit server -m mock.config 配置mock.config 支持正则 module.exports = { /queryProductDe ...
- CentOS7上Nginx的使用
Nginx 的启动 指定配置文件的方式启动nginx # nginx -c /etc/nginx/nginx.conf 对于yum安装的nginx,使用systemctl命令启动 # systemct ...
- C#将DataTable导出Execl、Word、Xml
/// <summary> /// 将DT转换为Execl的方法 /// </summary> /// <param name=" ...
- C#语法糖之第三篇: 匿名类 & 匿名方法
今天时间有点早,所以上来在写一篇文章吧,继续上一篇的文章,在我们平时编程过程中有没有遇到过这样的一个情景,你定义的类只是用来封装一些相关的数据,但并不需要相关联的方法.事件和其他自定义的功能.同时,这 ...
- [Mime] MediaTypes--电子邮件类型类 （转载）
点击下载 MediaTypes.rar 这个类是关于 电子邮件类型类的操作,在发送电子邮件是规定以什么样的格式发送,Xml,HTML,文本等方式1.电子邮件类型帮助类,Xml格式,HTML格式等看下面 ...
- Chart图形 [GDI+] OWCChart统计图的封装类 （转载）
点击下载 OWCChart.zip 利用OWC11进行作统计图的封装类. /// <summary> /// 类说明:进行作统计图的封装类 /// 联系方式:361983679 /// 更 ...
- js EasyUI前台 全选的实现
aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAAAWcAAAEQCAIAAADj/SKjAAAgAElEQVR4nO1dz48ry1W+/5N3swaFEP ...