C语言快速幂取余算法

#include<stdio.h>
#include<math.h>
//a为底数，b为次数，p为余数

long long fun(long long a, long long b,long long p)
{
	long long ans = 1;
 
	while (a&&b)
	{
		if (b & 1) ans = ((ans*a) % p) % p;
 
		a = (a*a) % p;
 
		b >>= 1; //移位运算，右移一位
	}
 
	return ans;
}

int main(){
	long long int n,p=1;
	scanf("%lld",&n);
	p=fun(2,n,20180520);
	printf("%lld",p);
}