逆元总结

逆元是个好东西，但太久不写就忘了。

一般来说，逆元可以用费马小定理来求解，但扩展GCD也可以。

费马小定理可以转化a/p%mod为a*(p^(mod-2))%mod，可以用快速幂。

long long ny(long long x,long long mod)
{
	long long cs=mod-2,anss=1;
	while(cs!=0)
	{
		if(cs%2!=0)
		{
			cs--;
			anss=((anss*x)%mod);
		}
		cs/=2;
		anss=((anss*x)%mod);
		anss=((anss*x)%mod);
	}
	return anss;
}

3777. 最短路（shortest） (Standard IO)
      Time Limits: 3000 ms Memory Limits: 262144 KB Detailed Limits
      Goto ProblemSet

Description
小Y最近学得了最短路算法，一直想找个机会好好练习一下。话虽这么说，OJ上最短路的题目都被他刷光了。正巧他的好朋友小A正在研究一类奇怪的图，他也想凑上去求下它的最短路。

   小A研究的图可以这么看：在一个二维平面上有任意点(x,y)（0<=x<=N,0<=y<=M,且x,y均为整数），且(x,y)向(x-1,y)（必须满足1<=x）和(x,y-1)（必须满足1<=y）连一条边权为0的双向边。

   每个点都有一个非负点权，不妨设(x,y)的权值为F[x][y]，则有：

   1.x=0或y=0：F[x][y]=1；2.其他情况：F[x][y]=F[x-1][y]+F[x][y-1]。

   现在，小Y想知道(0,0)到(N,M)的最短路，即使得经过的点的权值之和最小。为了炫耀自己学过最短路算法，他决定和你进行一场比赛，看谁的程序跑得快。然则小Y没有学过高精度算法，所以他希望输出答案时只输出答案模1000000007后的值。

Input
一行两个正整数N,M，表示图的大小。

Output
一行一个整数Ans，表示答案模1000000007后的值。

Sample Input
1 2

Sample Output
6

Data Constraint
见左图。

Hint
10%的数据满足N,M<=20；

30%的数据满足N,M<=100；

60%的数据满足min(N,M)<=100；

100%的数据满足N*M<=10^12。

#include<bits/stdc++.h>
using namespace std;
long long n,m,i,j,ans,f[2222][2222],dp[2222][2222],mod=1000000007,x=1,y=1;
long long ny(long long xx,long long yy)
{
    long long ba=xx,nw=1;
    while(yy!=0)
    {
        if(yy&1)
        {
        	nw=((nw*ba)%mod);
		}
        ba=((ba*ba)%mod);
        yy=yy/2;
    }
    return nw;
}
int main()
{
	scanf("%lld%lld",&n,&m);
	if(n>m) swap(n,m);
    for(i=m+2;i<=m+n+1;i++)
	{
		x=((x*i)%mod);
	}
    for(i=1;i<=n;i++)
	{
		y=((y*i)%mod); 
	}
    ans=((x*ny(y,mod-2)%mod+m)%mod);
	printf("%lld\n",ans%mod);
}