P5956-[POI2017]Podzielno【数学】
正題
題目鏈接:https://www.luogu.com.cn/problem/P5956
題目大意
BBB進(jìn)制下,給出序列aaa,aia_iai?表示數(shù)字iii有多少個(gè)。求一個(gè)最大的XXX在BBB進(jìn)制下,由給出的數(shù)字組成(不一定要用完),且其是B?1B-1B?1的倍數(shù)。
qqq次詢問XXX的第kkk位是幾。
2≤B≤106,1≤q≤105,1≤ai≤106,0≤k≤10182\leq B\leq 10^6,1\leq q\leq 10^5,1\leq a_i\leq 10^6,0\leq k\leq 10^{18}2≤B≤106,1≤q≤105,1≤ai?≤106,0≤k≤1018
解題思路
設(shè)xix_ixi?表示第iii位的話就是
(∑i=0xi×Bi)%(B?1)=0?∑i=0(xi×Bi%(B?1))=0\left(\sum_{i=0}x_i\times B^i\right)\%(B-1)=0\Rightarrow \sum_{i=0}\left(x_i\times B^i\%(B-1)\right)=0(i=0∑?xi?×Bi)%(B?1)=0?i=0∑?(xi?×Bi%(B?1))=0
拆開單獨(dú)的一個(gè)來看
xi×Bi%(B?1)=(xi%(B?1))×(Bi%(B?1))x_i\times B^i\%(B-1)=(\ x_i\%(B-1)\ )\times (\ B^i\%(B-1)\ )xi?×Bi%(B?1)=(?xi?%(B?1)?)×(?Bi%(B?1)?)
=xi%(B?1)×1=x_i\%(B-1)\times 1=xi?%(B?1)×1
所以其實(shí)就是各位數(shù)字的和為B?1B-1B?1的倍數(shù)就好了。
然后再回頭看題目發(fā)現(xiàn)有限制ai≥1a_i\geq 1ai?≥1。這樣如果用上所有數(shù)字的和對(duì)B?1B-1B?1取模為ttt的話,若ttt不為000,我們就讓ata_tat?減去一個(gè)111就好了。
然后對(duì)于詢問求一個(gè)前綴和然后二分
時(shí)間復(fù)雜度O(B+qlog?B)O(B+q\log B)O(B+qlogB)
code
#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=1e6+10; ll B,q,a[N]; signed main() {scanf("%lld%lld",&B,&q);ll t=0;for(ll i=0;i<B;i++){scanf("%lld",&a[i]);(t+=a[i]*i)%=B-1;}if(t)a[t]--;for(ll i=0;i<B;i++)a[i]+=a[i-1];while(q--){ll x;scanf("%lld",&x);x++;if(x>a[B-1])puts("-1");else printf("%lld\n",lower_bound(a,a+B,x)-a);}return 0; }總結(jié)
以上是生活随笔為你收集整理的P5956-[POI2017]Podzielno【数学】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: bzoj4025-二分图【线段树分治,并
- 下一篇: 网站怎么icp备案(网址icp如何备案)