赚钱买房(java大数+逆元+等比数列求和+快速幂)

这一题的思想就是等比数列求和,求和过程中需要用到快速幂方法!


求完和之后对这个数求模!求模用到了逆元,什么是逆元呢?


今天我们来探讨逆元在ACM-ICPC竞赛中的应用,逆元是一个很重要的概念,必须学会使用它。

 

对于正整数,如果有,那么把这个同余方程中的最小正整数解叫做的逆元。

 

逆元一般用扩展欧几里得算法来求得,如果为素数,那么还可以根据费马小定理得到逆元为

 

推导过程如下

                            

 

求现在来看一个逆元最常见问题,求如下表达式的值(已知

 

           

 

当然这个经典的问题有很多方法,最常见的就是扩展欧几里得,如果是素数,还可以用费马小定理。

 

但是你会发现费马小定理和扩展欧几里得算法求逆元是有局限性的,它们都会要求互素。实际上我们还有一

种通用的求逆元方法,适合所有情况。公式如下

 

         

 

现在我们来证明它,已知,证明步骤如下

 

         

 

这个就是逆元!


这个公式就是我们要求的数!

逆元就想上面哪个公式求!

#include<cstdio>
#include<cmath>
#include<iostream>
#define ll long long
using namespace std;
ll pow_mod(ll a,ll n,ll m)
{
    ll ans = 1;
    a = a % m;
    while(n>0)
    {
        if(n%2 ==1)
            ans=(ans*a)%m;
        n=n/2;
        a=(a*a) % m;
    }
    return ans;
}
int main()
{
    ll a,n,b,sum1,sum2,i,sum3;
    while(scanf("%lld%lld%lld",&a,&n,&b)!=EOF)
    {
       sum1=9*b;
       cout<<"模式:"<<pow_mod(10,n,sum1)<<endl;
       sum3=(pow_mod(10,n,sum1)-1%sum1)%sum1;
       cout<<sum3<<endl;
       sum2=((a%sum1)*(sum3%sum1))%sum1;
       sum2/=9;
       cout<<b-sum2<<endl;
    }
}
/*
0 9999999 2147483646
5 9999999  456465465
2012456840 1000 3
2054879653 1000 546813
2147483646 0 54
4568 9999999 5687
68456213 9999999 2147483646
9999999 9999999 9999999

2147483646
109893150
1
230594
54
2219
326455227
9999999
*/
C++在运算过程中超了ll的值,所以我们需要重新想办法!

这里我用了java大数!正好连第二题也给过了!

import java.math.BigInteger;
import java.util.Scanner;
import java.math.*;
public class Main {
	public static void main(String[] args)
	{
		Scanner cin=new Scanner(System.in);
		BigInteger sum2,a,b,sum1,sum3;
		int n;
		while (cin.hasNext()) 
		{
			a=cin.nextBigInteger();
			n=cin.nextInt();
			b=cin.nextBigInteger();
			sum1=BigInteger.valueOf(9).multiply(b);
			sum3=(pow_mod(BigInteger.valueOf(10),n,sum1).subtract((BigInteger.valueOf(1).remainder(sum1)))).remainder(sum1);
			sum2=((a.remainder(sum1)).multiply(sum3.remainder(sum1))).remainder(sum1);
			sum2=sum2.divide(BigInteger.valueOf(9));
			System.out.println(b.subtract(sum2));
		}
		
	}
	private static  BigInteger pow_mod(BigInteger a, int n,BigInteger m) 
	{
		BigInteger ans = BigInteger.valueOf(1);
	    a = a.remainder(m);
	    while(n>0)
	    {
	        if(n%2 ==1)
	            ans=(ans.multiply(a)).remainder(m);
	        n=n/2;
	        a=a.multiply(a).remainder(m);
	    }
	    return ans;
		
	}
}

发一下第二题的题目!

ok!呵呵!







  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值