题目描述

给出一个有理数c=\frac{a}{b},求c\ \bmod 19260817的值。

输入输出格式

输入格式:

一共两行。

第一行,一个整数a。
第二行,一个整数b。

输出格式:

一个整数,代表求余后的结果。如果无解,输出Angry!

输入输出样例

输入样例#1: 

233
666

输出样例#1: 

18595654

说明

对于所有数据,0\leq a,b \leq 10^{10001}


 个人思路:

  • 一开始还以为直接a/b就行了,后来发现范围太大了(而且\frac{a\boldsymbol{mod}P}{b\boldsymbol{mod}P}\neq \frac{a}{b}\boldsymbol{mod}P
  • 看了一下题解,思路是用费马小定理(a^{P-1}\equiv 1(\boldsymbol{mod }P))求逆元即可(因为p刚好是质数)
  • 也就是说,因为a*b=1,所以a*b=a^{P-1} \rightarrow a^{P-2}=b
  • 在获取数字的过程中,为了避开使用高精度,一位一位地读,每读完一位就%一次P即可保证数字不超范围(取模的基本性质)
  • 因为我们利用了取模的基本性质,所以最终答案当然不能直接a/b,所以我们需要用快速幂求b,也就是计算a^{P-2},也就计算出了b的值
  • 之后,a与b相乘并模P即可

#include<cstdio>
#include<iostream>
#include<cstring>
#define ll long long
using namespace std;
int p=19260817;
char c[10010];
int read(){
	scanf("%s",c);
	ll length=strlen(c);
	ll num=0;
	for(int i=0;i<length;i++){
		num*=10;
		num+=(c[i]-'0');
		num=num%p;
	}
	return num;
}
long poww(ll a,ll b){
	long ans=1;
	while(b){
		if(b&1){
			ans=(ans*a)%p;
		}
		a=(a*a)%p;
		b=b>>1;
	}
	return ans;
}
int main(){
	ll a,b;
	a=read(),b=read();
	if(b==0){
		cout<<"Angry!"<<endl;
		return 0;
	} 
	cout<<a*poww(b,p-2)%p<<endl;
	return 0;
}