P1630 求和(快速幂&规律)

求和

题目描述

1 b + 2 b + ⋯ + a b 1^b+2^b+\cdots + a^b 1b+2b++ab 的和除以 1 0 4 10^4 104 的余数。

输入格式

本题有多组数据。

第一行一个整数 N N N,表示共有 N N N 组测试数据。

对于每组数据,一行两个整数 a , b a,b a,b

输出格式

对于每组数据,一行一个整数,表示答案。

样例 #1

样例输入 #1

1
2 3

样例输出 #1

9

提示

对于 30 % 30\% 30% 的数据, N ≤ 10 N \le 10 N10 a , b ≤ 1 0 3 a,b \le 10^3 a,b103

对于 100 % 100\% 100% 的数据, 1 ≤ N ≤ 100 1 \le N \le 100 1N100 1 ≤ a , b ≤ 1 0 9 1 \le a,b \le 10^9 1a,b109

解析

可以观察到,将a分为多段,其中包含多段1-10000和一段10000以内的段。
前缀和,快速幂求出10000以内的前缀和,然后计算即可
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int mod=1e4;
const int N=1e4+5;
int a,b,sum[N];
int qmi(int m,int k,int p){		
	int res=1%p,t=m;
	while(k){
		if(k&1) res=res*t%p;
		t=t*t%p;
		k>>=1;
	}
	return res%mod;
}
void solve(){
	scanf("%lld%lld",&a,&b);
	for(int i=1;i<=1e4;i++){
		sum[i]=sum[i-1]+qmi(i,b,mod);
		sum[i]%=mod;
	}
	printf("%lld\n",(a/10000*sum[10000]+sum[a%10000])%mod);
}
signed main(){
	int t=1;
	scanf("%lld",&t);
	while(t--) solve();
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值