HDU 4704 Sum 【隔板原理+费马小定理+快速幂】

Sum

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 2745    Accepted Submission(s): 1143

Problem Description

 

Sample Input

  
  

   
   2
  
  

 

Sample Output

  
  

   
   2

   
   

    
    

     
     Hint
    
    
1.  For N = 2, S(1) = S(2) = 1.

2.  The input file consists of multiple test cases. 

   
   
   
    
  
  

 

Source

2013 Multi-University Training Contest 10

 

1.隔板原理
1~N有N个元素，每个元素代表一个1.分成K个数，即在(N-1)个空挡里放置（K-1）块隔板
即求组合数C(N-1,0)+C(N-1,1)+...+C(N-1，N-1)。
组合数求和公式：
C(n,0)+C(n,1)+C(n,2)+.+C(n,n)=2^n
证明如下：
利用二项式定理（a+b)^n=C(n,0)a^n+C(n,1)a^(n-1)b+C(n,2)a^(n-2)b^2 +.+C(n,n)b^n
令a=b=1左边就是2^n
所以题目即求2^（n-1）%（1e9+7）
设MOD为1e9+7
 
2.费马小定理（降幂）
a^(p-1)%p=1（a,p互质）
 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<algorithm>
using namespace std;
#define ll long long
#define ms(a,b)  memset(a,b,sizeof(a))
const int M=1e7+10;
const int inf=0x3f3f3f3f;
const int mod=1e9+7;
int i,j,k,n,m;
char s[M];

ll quickpow(ll a,ll b){

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

int main()
{
    while(~scanf("%s",s)){
     int len=strlen(s);
     ll num=0;
     for(int i=0;i<len;i++){
        num=(num*10+(s[i]-'0'))%(mod-1);
     }
      printf("%I64d\n",quickpow(2,num-1));
    }
    return 0;
}