牛客网暑期ACM多校训练营day4

A  根据题意找出递推式 发现上一项的ans作为下一项的幂次 所以欧拉降幂 (写错递归 背锅背锅

#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
#include<map>
#include<set>
#include<cstring>
#include<cstdio>
#define N 3000006
#define ll long long
using namespace std;
const int MAXN=2e5+7;
int cnt;
char s[MAXN];
int mul[N];
bool notprime[N];
int prime[N];
ll c1,c2,c3,c4,c5,c6,c0,c7;
long long qpow(long long a,long long b,long long c){
     long long t=1;long long e=a;
     while(b!=0){
        if((b&1)==1) t=t*e%c;
        e=e*e%c;
        b=b>>1;
     }
     return t;
}
  
void Prime(){
    notprime[1]=1;
    cnt=0;
    mul[1]=1;
    for(int i=2;i<N;i++){
        if(notprime[i]==0){
            prime[++cnt]=i;
            mul[i]=i-1;
        }
        for(int j=1;j<=cnt&&prime[j]*i<N;j++){
            if(notprime[i*prime[j]]==0)  notprime[i*prime[j]]=1;
            if(i%prime[j]==0){
                mul[i*prime[j]]=mul[i]*prime[j];break;
            }
            else{
                mul[i*prime[j]]=mul[i]*(prime[j]-1);
            }
        }
    }
}
ll dfs(int x, ll mod) {
    if(mod==1||x==0)return 0;
    if(s[x]=='2') {
        ll p;
        if(mod==c0)p=c1;
        else if(mod==c1)p=c2;
        else if(mod==c2)p=c3;
        else if(mod==c3)p=c4;
        else if(mod==c4)p=c5;
        else if(mod==c5)p=c6;
        else if(mod==c6)p=c7;
        else p=mul[mod];
        ll tmp=(qpow(2,(dfs(x-1,p)+p)%p,mod)*6LL%mod-3LL+mod)%mod;
        return tmp;
    }
    if(s[x]=='0') {
        return (dfs(x-1,mod)+1LL+mod)%mod;
    }
    if(s[x] =='1'){
        return (2LL*dfs(x-1,mod)%mod+2LL+mod)%mod;
    }
    return 0;
}
int main(){
    int m;
    Prime();
    scanf("%d",&m);
    c0=1000000007;c1=1000000006;c2=500000002;c3=243900800;c4=79872000;c5=19660800;c6=5242880;
    c7=2097152;
    int n;
    while(m--){
        scanf("%s",s+1);int len=strlen(s+1);
        printf("%lld
",dfs(len,1e9+7));
    }
    return 0;
}

 B DP

队友博客https://blog.csdn.net/qkoqhh/article/details/81267007

C 数位dp

队友博客https://blog.csdn.net/qkoqhh/article/details/81274954

D 构造 发现偶数项转移规律

#include<bits/stdc++.h>
using namespace std;
int a[205][205];
int n;
void ycl()
{
    a[1][1] = 0;a[1][2] = -1;
    a[2][1] = 1;a[2][2] = 1;
   for(int nn = 4;nn<=200;nn+=2)
   {
    for(int j = 1;j<=nn;j++) a[nn-1][j] = a[j][nn] = -1,a[nn][j] = a[j][nn-1] = 1;
    a[nn-1][nn-1] = 0;
   }
   return ;
}
int main()
{
    ycl();
    int _;cin>>_;
   while(_--)
   {
      cin>>n;
    if(n%2 == 1) cout<<"impossible"<<endl;
    else
    {
        cout<<"possible"<<endl;
        for(int i = 1;i<=n;i++)
        {
        for(int  j = 1;j<=n;j++)
        {
            cout<<a[i][j]<<" ";
        }  
        cout<<endl;
        }
    }      
   }
}

 E 题意题  ....看懂题以后 发现是扫描线线段树板子题  不想写离散化的辣鸡写动态开点调了一天...然后是没有取膜也没开ll 做法还是很明显:动态开点+扫描线线段树

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int MAXN=1e5+10;
const ll mod=1e9+7;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
ll ksm(ll a,ll b,ll c){
    ll ans=1;
    while(b){
    if(b&1)ans=ans*a%c;
    a=a*a%c;b=b>>1;
    }
    return ans;
}
typedef struct node{
    int l,r;ll sum,flag;
}node;
node d[MAXN*25];int cnt;ll ans=0;
void newnode(int &x,int len){x=++cnt;;d[x].l=d[x].r=0;d[x].sum=len;d[x].flag=1;}
void up(int x,int l,int r){
    int mid=(l+r)>>1;
    if(!d[x].l)d[x].sum=mid-l+1;else d[x].sum=d[d[x].l].sum*d[d[x].l].flag%mod;
    if(!d[x].r)d[x].sum+=r-mid;else d[x].sum+=d[d[x].r].sum*d[d[x].r].flag%mod;
    d[x].sum%=mod;
}
ll tag,tag1,vul1;
void update(int &x,int l,int r,int ql,int qr,ll vul){
    if(!x)newnode(x,r-l+1);
    if(ql<=l&&r<=qr){
    d[x].flag=d[x].flag*tag1%mod;
    vul1=vul*d[x].flag%mod;
    ans=(ans+tag*(r-l+1)%mod)%mod;
    ans=(ans+((-1ll*((vul1*(tag*(d[x].sum%mod)%mod))%mod))%mod+mod)%mod)%mod;
    return ;
 }
    int mid=(l+r)>>1;
    if(ql<=mid)update(d[x].l,l,mid,ql,qr,vul*d[x].flag%mod);
    if(qr>mid) update(d[x].r,mid+1,r,ql,qr,vul*d[x].flag%mod);
    up(x,l,r);
}
int n;
typedef struct Node{
    ll x,y;ll vul;
    friend bool operator<(Node aa,Node bb){
    if(aa.y==bb.y)return aa.x>bb.x;
    return aa.y>bb.y;
    }
}Node;
Node que[MAXN];
vector<int>vec;
void slove(){
    n=read();ll t1,t2;ll Xmax=-1,Ymax=-1;vec.clear();cnt=0;ans=0;
    for(int i=1;i<=n;i++)que[i].x=read(),que[i].y=read(),vec.push_back(que[i].y),t1=read(),t2=read(),Xmax=max(Xmax,que[i].x),Ymax=max(Ymax,que[i].y),que[i].vul=1ll*(t2-t1)*ksm(t2,mod-2,mod)%mod;
    sort(que+1,que+n+1);ll last=Xmax;tag=0;vec.push_back(0);
    sort(vec.begin(),vec.end());
    int sz=unique(vec.begin(),vec.end())-vec.begin();
    t1=sz-1;t2=sz-2;
    que[n+1].x=que[n+1].y=0;
    Ymax++;
    int tot=1;tag1=1ll;int rt=0;
    while(tot<=n){
    tag=vec[t1]-vec[t2];
    if(Xmax!=que[tot].x)update(rt,1,last,que[tot].x+1,Xmax,1ll);//cout<<que[tot].x+1<<"====="<<Xmax<<" "<<ans<<endl;
    Xmax=que[tot].x;tag1=tag1*que[tot].vul%mod;
    if(que[tot+1].y!=que[tot].y&&Xmax>=1){update(rt,1,last,1,Xmax,1ll);Xmax=0;}
    if(que[tot+1].y!=que[tot].y&&!Xmax){Xmax=last;tag1=1ll;t1=t2;t2--;}
    tot++;
    }
    printf("%lld
",ans);
}
int main(){
    int _;_=read();
    while(_--){
    slove();
    }
    return 0;
}

 F 签到

#include<bits/stdc++.h>
using namespace std;
char maps[2005][2005];
int n,m;
int main()
{
    ios::sync_with_stdio(false);
    int T;
    cin>>T;
    while(T--)
    {
    cin>>n>>m;
    for(int i = 1;i<=n;i++)
    {
        for(int j = 1;j<=m;j++) cin>>maps[i][j];
    }
    int heng= 0,shu = 0;
    for(int i = 1,j = n;i<j-1;i++,j--)
    {
    int flag = 0;
       for(int k = 1;k<=m/2;k++)
       {
       if(maps[i][k] != maps[i][m-k+1] || maps[j][k] !=maps[j][m-k+1]) {
           flag = 1;
           break;
       }
       }      
       if(flag == 0)
       {
       heng++;
       }
       else break;
    }
    for(int i = 1,j=m;i<j-1;i++,j--)
    {
    int flag = 0;
       for(int k = 1;k<=n/2;k++)
       {
       if(maps[k][i] != maps[n-k+1][i] || maps[k][j] !=maps[n-k+1][j]) {
           flag = 1;
           break;
       }
       }      
       if(flag == 0)
       {
       shu ++;
       }
       else break;
    }
       cout<<heng*shu<<endl;
    }
}

 G 贪心 二分前缀和即可

#include <bits/stdc++.h>
const int MAXN=1e5+10;
#define ll long long
using namespace std;
int d[MAXN],b[MAXN];
vector<int>vec;
ll num[MAXN];int a[MAXN];
int main(){
    int _;scanf("%d",&_);
    while(_--){
        int n,m;scanf("%d%d",&n,&m);
        int t;
        for(int i=1;i<=n;i++)scanf("%d",&a[i]),vec.push_back(a[i]);
        sort(vec.begin(),vec.end());
        int sz=unique(vec.begin(),vec.end())-vec.begin();
        for(int i=1;i<=n;i++)a[i]=lower_bound(vec.begin(),vec.begin()+sz,a[i])-vec.begin()+1,d[a[i]]++;
        for(int i=1;i<=sz;i++)b[i]=d[i];
        sort(d+1,d+sz+1);
        bool flag=0;
        for(int i=1;i<=sz;i++)num[i]=num[i-1]+d[i];
        for(int i=sz-1;i>=0;i--){
            int t=b[i+1];
            int pos=lower_bound(d+1,d+sz+1,t)-d;
            ll num1=num[sz]-num[pos-1]-t;
            ll t1=num1-1ll*(t-1)*(sz-pos);
            if(t1<=m){flag=1;printf("%d
",vec[i]);break;}
        }
        if(!flag)puts("-1");
        vec.clear();memset(num,0,sizeof(num));memset(d,0,sizeof(d));
    }
}

 J

一个比较神奇的题 首先明白每个不在对应位置的vul 他将收到一段区间的控制(具体是什么区间 自行理解) 如果这段区间最小值是-1 则是不存在的序列 然后想办法和这区间的点建图 但是明显 这个图建图是n^2的 考虑到分块的思想 结合线段树中的分块 我们按照线段树建树优化建图  通过拓扑解决问题(有环 则这个序列不成立)

#include <bits/stdc++.h>
#define ll long long
const int MAXN=2e5+10;
const int inf=1e9+7;
using namespace std;
int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int a[MAXN];
typedef struct node{
    int id,vul;
    friend bool operator<(node aa,node bb){
    return aa.vul>bb.vul;
    }
}node;
priority_queue<node>que;
int dp[MAXN][21];
int n;int ma[MAXN];
int key[MAXN<<2];
vector<int>vec[MAXN<<2];
int pos[MAXN];
int du[MAXN<<2];
void built(int rt,int l,int r){
    key[rt]=-1;du[rt]=0;
    if(l==r){pos[l]=rt;key[rt]=a[l];return ;}
    vec[rt<<1].clear();vec[rt<<1|1].clear();
    vec[rt<<1].push_back(rt);
    vec[rt<<1|1].push_back(rt);
    du[rt]+=2;
    int mid=(l+r)>>1;
    built(rt<<1,l,mid);
    built(rt<<1|1,mid+1,r);
}
void update(int rt,int l,int r,int ql,int qr,int vis){
    if(ql<=l&&r<=qr){vec[rt].push_back(vis);du[vis]++;return ;}
    int mid=(l+r)>>1;
    if(ql<=mid)update(rt<<1,l,mid,ql,qr,vis);
    if(qr>mid)update(rt<<1|1,mid+1,r,ql,qr,vis);
}
void st(){
    for(int i=1;i<=n;i++)dp[i][0]=a[i];
    for(int j=1;(1<<(j-1))<=n;j++){
    for(int i=1;i+(1<<j)<=n+1;i++){
        dp[i][j]=min(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
    }
    }
}
int rmq(int l,int r){
    if(l>r)return inf;
    int k=ma[r-l+1];int k1=(1<<k);
    return min(dp[l][k],dp[r-k1+1][k]);
}
vector<int>ans;
void bfs(){
    for(int i=1;i<=n;i++){
    if(!du[pos[i]]&&a[i]!=-1)que.push((node){pos[i],a[i]});
    }
    if(que.empty())return ;
    while(!que.empty()){
    node t=que.top();que.pop();
    if(t.vul!=-1)ans.push_back(t.vul);
    for(int i=0;i<vec[t.id].size();i++){
        du[vec[t.id][i]]--;
        if(!du[vec[t.id][i]])que.push((node){vec[t.id][i],key[vec[t.id][i]]});
    }
    }
    for(int i=1;i<=n;i++){
    if(du[pos[i]]){
        ans.clear();break;
    }
    }
}
int main(){
    ma[0]=-1;
    for(int i=1;i<MAXN;i++)if((i&(i-1))==0)ma[i]=ma[i-1]+1;else ma[i]=ma[i-1];
    int _;_=read();
    while(_--){
    n=read();ans.clear();
    while(!que.empty())que.pop();
    //memset(num,0,sizeof(num));
    int cnt=0;
    for(int i=1;i<=n;i++)a[i]=read();
    built(1,1,n);
    st();bool flag=0;
    for(int i=1;i<=n;i++){
        if(a[i]==-1)continue;
        cnt++;
        if(a[i]%n==i-1){continue;}
        int t1=a[i]%n+1;
        if(i>t1){
        int t2=rmq(t1,i-1);
        if(t2==-1){flag=1;break;}
        update(1,1,n,t1,i-1,pos[i]);
        }
        else{
        int t2=min(rmq(t1,n),rmq(1,i-1));
        if(t2==-1){flag=1;break;}
        update(1,1,n,t1,n,pos[i]);
        if(i!=1)update(1,1,n,1,i-1,pos[i]);
        }
    }
    if(!cnt){puts("");continue;}
    if(flag){puts("-1");continue;}
    bfs();
    if(ans.size()<cnt){puts("-1");continue;}
    for(int i=0;i<ans.size();i++){
        if(i==0)printf("%d",ans[i]);
        else printf(" %d",ans[i]);
    }
    printf("
");
    }
    return 0;
}