codeforces359C 2000分数学

题目传送门

题意:

n个数,第i个数是 \dpi{150}a_i 。一个素数x。

表示n个分数之和,即 \sum_{i=1}^n \frac{1}{x^{a_i}} 。

把n个分数之和表示为 \frac{t}{s} 的形式,其中 s = x^{\sum_{i = 1}^n a_i} 。

计算 gcd(s , t)\; \% \;mod 。

题解:

\frac{t}{s} = \frac{\sum_{i=1}^{n} x ^ {sum-a_i}}{x^{sum}} , 其中 sum = \sum_{i=1}^{n}a_i 。

你可以整理成 \frac{x^y * m}{x^{sum}} , 其中m不是x的倍数。

你可以把分子看成是一个多项式,按照幂次从低到高去整理系数,如果系数不能整除x,那么你找到了y。

用map去整理系数方便一些,一直处理map的首元素,map的第一关键字表示幂次,第二关键字表示系数。

答案是 x^{min(y,sum)} \;\%\;mod ,快速幂计算。

感受:

写好后,直接编译成功,直接AC。太少见了。

代码:

#include<bits/stdc++.h>
using namespace std ;
typedef long long ll ;
typedef pair<int , int> pii ;
const int maxn = 1e5 + 5 ;
const ll mod = 1e9 + 7 ;
int n ;
ll x , a[maxn] , s = 0 ;
map<ll , ll> num ;
ll qpow(ll a , ll b)
{
    ll ans = 1 ; 
    while(b)
    {
        if(b & 1)
          ans = (ans * a) % mod ;
        b >>= 1 ;
        a = (a * a) % mod ;
    }
    return ans % mod ;
}
void solve()
{
	ll ans = 0 ;
	while(1)
	{	
		ll y = (*num.begin()).first ;
		ll z = num[y] ;
		if(z % x != 0){ans = y ; break ;}
		num.erase(y) ;
		while(z % x == 0)
		{
			z /= x ;
			y ++ ;
		}
		num[y] += z ;
	}
	ans = min(ans , s) ;
	printf("%lld\n" , qpow(x , ans)) ;
}
int main()
{
	scanf("%d%lld" , &n , &x) ;
	for(int i = 1 ; i <= n ; i ++)  scanf("%lld" , &a[i]) , s += a[i] ;
	for(int i = 1 ; i <= n ; i ++)  num[s - a[i]] ++ ;  
	solve() ;
	return 0 ;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值