自然数幂之和

简单的说，问题就是：

不同的问题随着k的大小变化解法也要求相应的进化。

zoj A very easy task

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1864

相关数据： (0 <= n <= 10⁵⁰, 1 <= k <= 100)

二项式展开：
进一步得到：
把n变成n-1,n-2,''',1：




上式累加：

换元：

k不断减小，直到递归到k=1，那么S(n,1)=(n+1)n/2; 时间复杂度：O(n^2)
import java.util.*;
import java.math.BigInteger;
public class Main {
	static final int N=110;
    static BigInteger []ans=new BigInteger [N];
    static BigInteger [][] C=new BigInteger[N][N];
    static void init(){
    	for(int i=1;i<N;i++){
    		C[i][0]=C[i][i]=BigInteger.ONE;
    		for(int j=1;j<i;j++){
    			C[i][j]=C[i-1][j-1].add(C[i-1][j]);
    		}
    	}
    }
    static BigInteger solve(BigInteger n,int k){
    	if(ans[k].compareTo(BigInteger.valueOf(-1))!=0){
    		return ans[k];
    	}
    	if(k==1){
    		ans[k]=(n.add(BigInteger.ONE)).multiply(n).divide(BigInteger.valueOf(2));
    	    return ans[k];
    	}
    	BigInteger t=BigInteger.ONE,temp=BigInteger.ZERO;
    	for(int i=0;i<k+1;i++){
    		t=t.multiply(n.add(BigInteger.ONE));
    	}
    	t=t.subtract(n.add(BigInteger.ONE));
    	for(int i=1;i<k;i++){
    		temp=temp.add(C[k+1][k+1-i].multiply(solve(n,i)));
    	}
    	t=t.subtract(temp);
    	ans[k]=t.divide(BigInteger.valueOf(k+1));
    	return ans[k];
    }
	public static void main(String[] args) {
		init();
		BigInteger n;
		int k;
		Scanner cin=new Scanner (System.in);
		while(cin.hasNext()){
			n=cin.nextBigInteger();
			k=cin.nextInt();
			for(int i=0;i<N;i++){
				ans[i]=BigInteger.valueOf(-1);
			}
			System.out.println(solve(n,k));
		}

	}

}

51nod 1228 序列求和

http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1228

相关数据：(1 <= N <= 10^18, 1 <= K <= 2000)

T(n) = n^k，S(n) = T(1) + T(2) + ...... T(n)。给出n和k，求S(n)。
 
 

  
  例如k = 2，n = 5，S(n) = 1^2 + 2^2 + 3^2 + 4^2 + 5^2 = 55。
 
 
 
 

  
  由于结果很大，输出S(n) Mod 1000000007的结果即可。
 
 
 
 

  
  分析：有5000case，自然数幂之和的普通解法O(k^2)，故另辟蹊径。借助ACdreamer大神的思路，该问题和伯努利数联系。有这样的等式：

 
 
 
 

  
  顺便回忆了一下逆元的相关知识，万能的求解，这里需要K=1-->2000，k+1 mod 1000000007的逆，可以使用递推逆元的方法预处理：

 
 
 
 

  
  相关证明：http://blog.csdn.net/acdreamers/article/details/8220787
 
 
 
 

  
  
  
  

   
   组合数C(n,m)同样预处理。
  
  
  
  

   
   最后的伯努利：
  
  
  
  

 
 
 
 

  
  mod操作，小心又小心，不要少加%mod。

 
 
 
 

  
  
  
  
import java.util.*;
public class Main {
 static final int N=2005,mod=1000000007;
    static long [][]C=new long [N][N];
 static long []inv=new long [N],B=new long [N],nadd1=new long [N];
    static void initC(){
  C[0][0]=1;
  for(int i=1;i<N;i++){
   C[i][0]=C[i][i]=1;
   for(int j=1;j<i;j++)
       C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
  }
 }
 static void initinv(){
  inv[1]=1;
  for(int i=2;i<N;i++){
   inv[i]=(mod-mod/i)*inv[mod%i]%mod;
  }
 }
 static void initB(){
  B[0]=1;
  for(int i=1;i<N;i++){
   if(i==N-1) break;
   for(int j=0;j<i;j++){
    B[i]=(B[i]+C[i+1][j]*B[j]%mod)%mod;
   }
   B[i]=B[i]*(-inv[i+1]);
   B[i]=(B[i]%mod+mod)%mod;
  }
 }
    public static void main(String[] args) {
  int t,k;
  long n;
  Scanner cin=new Scanner (System.in);
  t=cin.nextInt();
  initC();  
  initinv();
  initB();  
  while(t>0){
   t--;
   n=cin.nextLong();
   n=n%mod; //!!
   k=cin.nextInt();
   nadd1[0]=1;
   for(int i=1;i<=k+1;i++){
    nadd1[i]=nadd1[i-1]*(n+1)%mod;
   }
   long ans=0;
   for(int i=1;i<=k+1;i++){
    ans=(ans+C[k+1][i]*B[k+1-i]%mod*nadd1[i]%mod)%mod;
   }
   ans=(ans*inv[k+1]%mod+mod)%mod;
   System.out.println(ans);
  }
 }

}

  
  当k达到万级时，又是一个工程，我还得继续学习。。。||-_-