B

題解

$f[i][(gcd(prime[j]*prime[k]\%P,P))]=\sum\limits_{k=1}^{k<=num} f[i-1][k]*phi(\frac{P}{prime[j]})$

關于$phi(\frac{P}{prime[j]})$理解

$phi(\frac{P}{prime[j]})$是求$prime[j]$代表的數的個數

$P=k_0*prime[j]$

$x_1=k_1*prime[j]$

$x_2=k_2*prime[j]$

.......

?

要求代表$prime[j]$數個數就是求$k_1$,$k_2$個數$(k_0,k_1,k_2等互質)$(不互質$gcd就是別的數了$)

移項顯然與$k_0$互質數個數就是$phi(\frac{P}{prime[j]})$

代碼

?

#include<bits/stdc++.h> using namespace std; #define ll long long #define A 55 const ll mod=1e9+7; ll f[A][23333],phi[A*A],prm[A*A],to[A*A][A*A]; ll n,m,P; ll meng(ll x,ll k){ll ans=1;for(;k;k>>=1,x=x*x%mod)if(k&1)ans=ans*x%mod;return ans; } ll gcd(ll x,ll y){if(y==0) return x;return gcd(y,x%y); } ll p(ll x){ll ans=x;for(ll i=2;i*i<=x;i++){if(x%i==0){ans=ans/i*(i-1);while(x%i==0) x/=i;}}if(x>1) ans=ans/x*(x-1);return ans; } void fen(ll x){for(ll i=1;i*i<=x;i++){if(x%i==0){prm[++prm[0]]=i;if(i*i!=x) prm[++prm[0]]=x/i;}}sort(prm+1,prm+prm[0]+1); } void pre_work(){fen(P);for(ll i=1;i<=prm[0];i++)phi[i]=p(P/prm[i]);for(ll j=1;j<=prm[0];j++){for(ll pre=1;pre<=prm[0];pre++){ll g=gcd(prm[j]*prm[pre]%P,P);to[j][pre]=lower_bound(prm+1,prm+prm[0]+1,g)-prm;}} } void work(){for(ll i=2;i<=n;i++)for(ll j=1;j<=prm[0];j++)for(ll pre=1;pre<=prm[0];pre++)(f[i][to[j][pre]]+=f[i-1][pre]*phi[j]%mod)%=mod; } void sub_task(){pre_work();for(ll j=1;j<=prm[0];j++)f[1][j]=phi[j];work();for(ll i=1,a;i<=m;i++){scanf("%lld",&a); a=lower_bound(prm+1,prm+prm[0]+1,gcd(a,P))-prm;printf("%lld ",f[n][a]*meng(phi[a],mod-2)%mod)%mod;}printf("\n"); } int main(){scanf("%lld%lld%lld",&n,&m,&P);sub_task();return 0; } View Code

?

C

題解

三分,對于怎么看出來三分,這可能是個套路,你覺得這個題你用貪心做不了(但非常像貪心),二分答案會被hack,然后你$dp$也難以做,你三分就可以了

三分$check$貪心做,很水,瞎jb差分一下,我會說貪心我考試時就寫對了嗎?

注意細節,細節很多,不要死于細節

代碼

/* n*log^2 */ #include<bits/stdc++.h> using namespace std; #define ll long long #define A 222222 struct node{ll l,r;friend bool operator < (const node & a,const node &b){return a.l==b.l?a.r>b.r:a.l<b.l;} }wat[A]; ll n,m,t,ans=0x7fffffffffffff,maxx=0,cnt=0; ll p[A],lef[A],now[A],c[A]; ll check(ll x){ll sum=0;for(ll i=1;i<=n;i++){now[i]=max(p[i]-x,0ll);c[i]=0;}for(ll i=1;i<=n;i++){ll nowid=lef[i];c[i]+=c[i-1];if(lef[i]==0) continue ;ll cha=c[i]; // printf("i=%lld c[i]=%lld c[i-1]=%lld\n",i,c[i],c[i-1]); // printf("i=%lld now+cha=%lld\n",i,now[i]+cha);if(now[i]+cha>0){sum+=now[i]+cha;c[i]-=(now[i]+cha);c[wat[nowid].r+1]+=now[i]+cha;now[i]=0;}}for(ll i=1;i<=n;i++){ll cha=c[i]; // printf("now=%lld x=%lld cha=%lld\n",now[i],x,cha);if(now[i]+cha>0) return 0x7fffffffff; }return sum+x*t; } int main(){ // freopen("da.in","r",stdin); // freopen("ans.bf","w",stdout);scanf("%lld%lld%lld",&n,&m,&t);for(ll i=1;i<=n;i++){scanf("%lld",&p[i]);maxx=max(maxx,p[i]);}for(ll i=1;i<=m;i++){scanf("%lld%lld",&wat[i].l,&wat[i].r);}sort(wat+1,wat+m+1);for(ll i=1;i<=m;i++){if(!lef[wat[i].l])lef[wat[i].l]=i;}for(ll i=1;i<=n;i++)if(wat[lef[i-1]].r>=i){if(wat[lef[i-1]].r>wat[lef[i]].r)lef[i]=lef[i-1]; }ll l=0,r=maxx;while(l<r){ll len=(r-l);ll lmid=l+len/3,rmid=r-len/3;ll lnow=check(lmid),rnow=check(rmid); // printf("l=%lld r=%lld\n",l,r);if(lnow>=rnow) l=lmid+1;else r=rmid-1; ans=min(ans,lnow);ans=min(ans,rnow);} // printf("%lld\n",check(5));printf("%lld\n",ans); } View Code

?

轉載于:https://www.cnblogs.com/znsbc-13/p/11524078.html

總結

以上是生活随笔為你收集整理的csp-c模拟测试43「A·B·C」的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。