HDU - 4704 Sum （费马小定理 + 快速幂）

最新推荐文章于 2020-07-08 16:12:02 发布

Joyyiwei

最新推荐文章于 2020-07-08 16:12:02 发布

阅读量810

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/u011345136/article/details/38389693

版权

数论专栏收录该内容

181 篇文章 0 订阅

订阅专栏

Description

Sample Input

Sample Output

Hint

1. For N = 2, S(1) = S(2) = 1. 2. The input file consists of multiple test cases.

题意：求1-n中，组成n的不同种数

思路：首先我们要先推出这个公式：2^(n-1)是结果，结果的推论是：我们想象n个1排一排，那么我们能用隔板法得到的结果是C(n-1, 0), C(n-1, 1) ..... 所以结果就是2^(n-1)了

但是n那么大，显然不是我们能处理来的，所以这里用到了

费马小定理：：假如a是一个整数，p是一个质数，如果a不是p的倍数，这个定理也可以写成 $a^{p-1} \equiv 1 \pmod{p}$

所以我们将n拆成多个 a*(p-1) + k 的话，那么2^n = 2^[a*(p-1) + k ] % M = 2^k % M,然后就是再减一处理的快速幂取模了

#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <algorithm>
typedef long long ll;
using namespace std;
const ll M = 1e9+7;

char str[1000010];

ll cal(ll m) {
	int len = strlen(str);
	ll ans = 0;
	for (int i = 0; i < len; i++) 
		ans = (ans * 10 + str[i] - '0') % m;
	return ans;
}

ll pow_mod(ll a, ll b) {
	ll ans = 1;
	while (b) {
		if (b & 1)
			ans = ans * a % M;
		a = a * a % M;
		b >>= 1;
	}
	return ans;
}

int main() {
	while (scanf("%s", str) != EOF) {
		ll k = cal(M-1);
		cout << pow_mod(2, k-1) << endl;
	}
	return 0;
}