bzoj 1089

很久没有写博客了呢。

今天来补一补。

ps：

要摆脱题解依赖了！

∑(っ °Д °;)っ (๑•̀ㅂ•́)و✧

先说这道dp方程十分简单，却要写高精的题。

这个dp【i】表示的是深度小于等于i的数目。所以最后结果需要剪掉它

找了一个重载运算符的代码。感觉真是异常的高深。学到了很多东西。

然后代码比别人写的短哈哈哈

再补充一句：

自己的一些小理解：

高精就是把一个数的每一位当成一个单独的字符进行运算，他们之间的唯一联系就在于进位。

而此题由于数很多。我们把三个字符砍成一截进行存储。所以就有了%1000这种东西。

其实本nc宝宝 一开始就有点翻不过来，尤其是到输出还有倒着输出。

但是后来你想想就明白了在单个的一节中 我们读肯定是 百十个、

输出顺序也是百十个。但是你的n节从左往右开始递增意思是右边最大。

先从右开始，在每个节中输出又是从做开始的所以正好啊。

！还有一点要注意的是每个字节都要沾满3个字符，不够就先补0

/**************************************************************
Problem: 1089
User: mars_ch
Language: C++
Result: Accepted
Time:36 ms
Memory:1476 kb
****************************************************************/

#include<cstdio>
#include<algorithm>
#include<cstring>
#define rad 1000
//by mars_ch
using namespace std;
int n,d;
struct data
{
int l;
int v[5005];
}dp[31];
data operator*(data a,data b)     //重新定义乘法的算法
{
data ans;
for(int i=1;i<=a.l+b.l;i++) ans.v[i]=0;  //!!!!!!    没进行一会新的乘法，不能够累加上上回的结果
for(int i=1;i<=a.l;i++)
{
for(int j=1;j<=b.l;j++)
{
ans.v[i+j-1]+=a.v[i]*b.v[j];    //先按乘法原理乘一遍，就跟普通的一样
}
}
ans.l=a.l+b.l;    //最多吗
for(int i=1;i<=ans.l;i++)
{
if(ans.v[i]>=rad)     //开始按3位处理了
{
if(i == ans.l) ans.l++;
ans.v[i+1]+=ans.v[i]/rad;
ans.v[i]%=rad;
}
}
while(ans.l>0 && ans.v[ans.l] == 0) ans.l--;     如果前导是0删掉。
return ans;
}
data operator+(data a,int p)
{
a.v[1]+=p;
int id=1;
while(a.v[id]>=rad)
{
a.l=max(id+1,a.l);
a.v[id+1]+=a.v[id]/rad;
a.v[id]%=rad;
id++;
a.l=max(id,a.l);
}
return a;

}
data operator-(data a,data b)
{
for(int i=1;i<=a.l;i++)
{
a.v[i]-=b.v[i];
if(a.v[i]<0)
{
a.v[i]+=rad;
a.v[i+1]--;
}
}
while(a.l>0 && a.v[a.l]==0) a.l--;
return a;
}
data operator^(data a,int p)
{
data ans;
ans.l=1;
ans.v[1]=1;
for(int i=p;i>=1;i--)
{
ans=ans*a;
}
return ans;
}
void print(data a)
{
printf("%d",a.v[a.l]);
for(int i=a.l-1;i>=1;i--)
{
printf("%03d",a.v[i]);
}
printf("\n");
}
int main()
{
scanf("%d%d",&n,&d);
if(d == 0)
{
puts("1");
return 0;
}
dp[0].l=1;
dp[0].v[1]=1;
for(int i=1;i<=d;i++)
{
dp[i]=(dp[i-1]^n)+1;
}
print(dp[d]-dp[d-1]);
return 0;
}