题目描述
我们知道,整数做除法时,有时得到有限小数,有时得到无限循环小数。
如果我们把有限小数的末尾加上无限多个0,它们就有了统一的形式。
本题的任务是:在上面的约定下,求整数除法小数点后的第n位开始的3位数。
输入
一行三个整数:a b n,用空格分开。a是被除数,b是除数,n是所求的小数后位置(0<a,b,n<1000000000)
输出
一行3位数字,表示:a除以b,小数后第n位开始的3位数字。
样例输入
1 8 1
样例输出
125
思路:a/b的值,理论上都可以表示成无限循环小数。为什么呢?根据鸽巢原理,在a/b的次数大于等于b次的时候,就肯定会出现重复。那么也肯定就是可以表示成不循环部分+循环部分了。那么再根据n的关系去找出对应位置的数。
代码如下:
#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll a,b,n;
int main()
{
scanf("%lld%lld%lld",&a,&b,&n);
string s="";
a%=b;
int pos=-1;
map<ll,int> mp;
int cnt=0;
while(1)
{
a*=10ll;
char c=(char)(a/b+'0');
s+=c;
a%=b;
if(mp[a])
{
pos=mp[a];
break;
}
else mp[a]=++cnt;
}
int len=s.length()-pos;
if(n>pos)
{
n-=(pos);
n%=len;
n=pos+(n?n-1:len-1);
}
else n--;
for(int i=1;i<=3;i++)
{
cout<<s[n];
n++;
if(n>=s.length()) n=pos;
}
return 0;
}
说句实话,这个题目,蓝桥的数据水了,b没有到1e9的,因此才不会超时,否则这个算法,也会存在超时的风险的。
自我感觉的正解:
因为求第n位开始的三位小数,那么我们只求到n+2就可以了。
那么也就是求(a/b)*10^(n+2)%1000就可以了。
1000并为素数所以不能用费马定理和扩展欧几里得。需要用下面的方法来实现:
(a/b)*10^(n+2)%1000=(a *10 ^ (n+2))%(b * 1000)/b;
那么代码的主要内容就是求一个快速幂即可;
代码如下:
#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll a,b,n;
inline ll pow(ll x,ll y,ll mod)
{
ll ans=1ll;
while(y)
{
if(y&1) ans=(ans*x)%mod;
x=(x*x)%mod;
y>>=1;
}
return ans;
}
int main()
{
scanf("%lld%lld%lld",&a,&b,&n);
ll mod=1000*b;
ll temp=pow(10ll,n+2,mod);
ll c=(a*temp)%mod/b;
printf("%03d\n",c);
return 0;
}
努力加油a啊,(o)/~