2018 Multi-University Training Contest 4

太尼玛难了，告辞。
A
x是k-magic数：小于x的正整数有超过K个因子多于x。
重要结论：一个数如果是k-magic数，那么它的倍数也是k-magic数。
发现非k-magic数很少——枚举k最大时的所有非k-magic数。
一开始知道1是非k-magic数。从小到大枚举质数。
对于每个质数p，将已知的数乘以p的幂加入数列直到超出上限，然后筛掉全部k-magic数（按数字从小到大维护前k多因数数量）。
直到连续几个p没有改变数列停止。此时p约为293。
在k是233的数列中以同样的方法筛出k是0到233的234个数列。
二分。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const ll maxn=1.1e18;
ll t,k;
ll n;
bool isprime[300];
ll prime[300],p;
ll ji[250];

struct number{
    ll x,num;
}a[1000000],tmpa[1000000],nums[235][50000];

ll pa,tmppa,ps[250];

bool cmp(number a, number b){
    return a.x<b.x;
}

void init(){
    p=0;
    memset(isprime,1,sizeof(isprime));
    for(ll i=2; i<300; i++){
        if(isprime[i]){
            prime[p++]=i;
            for(ll j=i*i; j<300; j+=i)
                isprime[j]=0;
        }
    }
    pa=1;
    a[0]=(number){1,1};
    for(int i=0; i<p; i++){
        tmppa=0;
        for(int j=0; j<pa; j++){
            ll tmp=a[j].x;
            tmpa[tmppa++]=a[j];
            for(int c=1; tmp<maxn/prime[i]; c++){
                tmp*=prime[i];
                tmpa[tmppa++]=(number){tmp,a[j].num*(c+1)};
            }
        }
        sort(tmpa,tmpa+tmppa,cmp);
        memset(ji,0,sizeof(ji));
        ji[0]=maxn;
        pa=0;
        for(ll j=0; j<tmppa; j++){
            if(tmpa[j].num>=ji[234]){
                a[pa++]=tmpa[j];
                for(ll k=234; k>0; k--){
                    ji[k]=ji[k-1];
                    if(tmpa[j].num<ji[k]){
                        ji[k]=tmpa[j].num;
                        break;
                    }
                }
            }
        }
    }
    memset(ji,0,sizeof(ji));
    ji[0]=maxn;
    memset(ps,0,sizeof(ps));
    for(ll i=0; i<pa; i++){
        ll j=234;
        while(a[i].num>=ji[j-1]) j--;
        //cout<<j<<endl;
        for(ll k=j-1; k<=233; k++)
            nums[k][ps[k]++]=a[i];
        for(ll k=234; k>j; k--) ji[k]=ji[k-1];
        ji[j]=a[i].num;
    }

    //cout<<pa<<":"<<endl;
    //for(int i=0; i<100; i++) cout<<a[i].x<<" "<<a[i].num<<endl;
    //cout<<a[840].x<<" "<<a[841].x<<endl;
    //cout<<nums[233][840].x<<" "<<nums[233][841].x<<endl;
    //cout<<nums[100][286].x<<" "<<nums[100][287].x<<endl;
    return;
}

int main(){
    init();
    scanf("%lld",&t);
    while(t--){
        scanf("%lld%lld",&n,&k);
        ll l=0,r=ps[k]-1,mid;
        while(l<r-1){
            //cout<<l<<" "<<r<<endl;
            mid=(l+r)/2;
            if(n+mid+1<=nums[k][mid].x)
                r=mid;
            else l=mid;
        }
        ll ans=n+l+1;
        printf("%lld\n",ans);
    }
    return 0;
}

B
发现可以$O(1)$转移，然后莫队一下。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
ll power(ll n,ll p){
    ll ans=1;
    ll base=n;
    while(p){
        if(p&1){
            ans=ans*base%mod;
        }
        base=base*base%mod;
        p>>=1;
    }
    return ans;
}
ll getinv(ll n){
    return power(n,mod-2);
}
int T;
int bsiz=sqrt(100000);
struct node{
    int x,y,id;
}rec[100003];
bool cmp(node a,node b){
    if(a.x/bsiz==b.x/bsiz){
        if((a.x/bsiz)%2)return a.y<b.y;
        return a.y>b.y;
    }
    return a.x<b.x;
}
int tot=0;
ll fac[100003];
ll inv[100003];
inline ll C(ll n,ll m){
    return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
ll tmp=1;
ll tx,ty;
ll ans[100003];
int main(){
    fac[0]=inv[0]=1;
    for(int i=1;i<=100000;i++){
        fac[i]=fac[i-1]*i%mod;
        inv[i]=getinv(fac[i]);
    }
    scanf("%d",&T);
    while(T--){
        ++tot;
        scanf("%d%d",&rec[tot].x,&rec[tot].y);
        rec[tot].id=tot;
    }
    sort(rec+1,rec+1+tot,cmp);
    tx=1,ty=0;
    for(int i=1;i<=tot;i++){
        while(tx<rec[i].x){
            tmp=(2*tmp-C(tx,ty)+mod)%mod;
            tx++;
        }
        while(ty>rec[i].y){
            tmp=(tmp-C(tx,ty)+mod)%mod;
            ty--;
        }
        while(tx>rec[i].x){
            tmp=(tmp+C(tx-1,ty))*inv[2]%mod;
            tx--;
        }
        while(ty<rec[i].y){
            tmp=(tmp+C(tx,ty+1))%mod;
            ty++;
        }
        ans[rec[i].id]=tmp;
    }
    for(int i=1;i<=tot;i++)printf("%lld\n",ans[i]);
}

C
并查集维护边都是1的联通块和大小，边是1和2的联通块和大小，边是1的联通块相邻的用3连接的边是1和2的联通块的总大小。最后一个可以把每个联通块加到树的father上，查询的时候只要查询儿子的总和以及father的贡献，于是变成了维护两种联通块的边是3的直接father。细节很多。

#include <bits/stdc++.h>
using namespace std;
map<pair<int,int> ,int>mp;
struct edge{
    int to,next,v;
}e[200003];
int head[100003];
int cnt;
void init(int n){
    for(int i=0;i<=n;i++)head[i]=-1;
    cnt=0;
}
void add(int u,int v,int w){
    e[cnt].to=v;
    e[cnt].next=head[u];
    e[cnt].v=w;
    head[u]=cnt++;
}
int n,m;
int f1[100003],f2[100003],up1[100003],up2[100003];
int s1[100003],s2[100003],s3[100003];
int getfa1(int u){
    if(u==f1[u])return u;
    return f1[u]=getfa1(f1[u]);
}
int getfa2(int u){
    if(u==f2[u])return u;
    return f2[u]=getfa2(f2[u]);
}
void Union1(int u,int v){
    u=getfa1(u);
    v=getfa1(v);
    f1[u]=v;
    s1[v]+=s1[u];
    s3[v]+=s3[u];
    up1[v]=max(up1[u],up1[v]);
}
void Union2(int u,int v){
    u=getfa2(u);
    v=getfa2(v);
    f2[u]=v;
    if(up2[u]>up2[v]){
        int fa=getfa1(up2[u]);
        s3[fa]+=s2[v];
    }
    else if(up2[u]<up2[v]){
        int fa=getfa1(up2[v]);
        s3[fa]+=s2[u];
    }
    s2[v]+=s2[u];
    up2[v]=max(up2[u],up2[v]);
}
void update(int u,int v){
    if(u>v)swap(u,v);
    pair<int,int> pr=make_pair(u,v);
    int typ=mp[pr];
    if(mp[pr]==1)return;
    else if(mp[pr]==2){
        Union1(u,v);
        mp[pr]--;
    }
    else {
        int u1=getfa1(u);
        int u2=getfa2(u);
        int v1=getfa1(v);
        int v2=getfa2(v);
        if(up1[u1]==v){
            s3[v1]-=s2[u2];
            up1[u1]=up2[u2]=0;
        }
        else {
            s3[u1]-=s2[v2];
            up1[v1]=up2[v2]=0;
        }
        Union2(u,v);
        mp[pr]--;
    }
}
void query1(int u,int v){
    int u2=getfa2(u);
    int v1=getfa1(v);
    if(getfa1(up2[u2])==v1||getfa2(up1[v1])==u2||getfa2(u)==getfa2(v)){
        printf("1 ");
    }
    else {
        printf("0 ");
    }
}
void query2(int u){
    int u1=getfa1(u);
    int u2=getfa2(u);
    int ans=s3[u1]+s2[u2];
    int fa=up1[u1];
    if(fa!=0){
        int fa2=getfa2(fa);
        ans+=s2[fa2];
    }
    printf("%d\n",ans);
}
void dfs(int u,int fa){
    for(int i=head[u];~i;i=e[i].next){
        int v=e[i].to;
        int typ=e[i].v;
        if(v==fa)continue;
        if(typ==3){
            int v1=getfa1(v),v2=getfa2(v);
            up1[v1]=u;
            up2[v2]=u;
            s3[getfa1(u)]+=s2[v2];
        }
        dfs(v,u);
    }
}
void debug(){
    cout<<"f1:  ";
    for(int i=1;i<=n;i++)printf("%d ",f1[i]);
    cout<<endl<<"s1:  ";
    for(int i=1;i<=n;i++)printf("%d ",s1[i]);
    cout<<endl<<"up1: ";
    for(int i=1;i<=n;i++)printf("%d ",up1[i]);
    cout<<endl<<"f2:  ";
    for(int i=1;i<=n;i++)printf("%d ",f2[i]);
    cout<<endl<<"s2:  ";
    for(int i=1;i<=n;i++)printf("%d ",s2[i]);
    cout<<endl<<"up2: ";
    for(int i=1;i<=n;i++)printf("%d ",up2[i]);
    cout<<endl<<"s3:  ";
    for(int i=1;i<=n;i++)printf("%d ",s3[i]);
    cout<<endl<<endl;
}
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        mp.clear();
        scanf("%d%d",&n,&m);
        init(n);
        for(int i=1;i<=n;i++){
            f1[i]=f2[i]=i;
            s1[i]=s2[i]=1;
            s3[i]=0;
            up1[i]=up2[i]=0;
        }
        for(int i=1;i<n;i++){
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            if(u>v)swap(u,v);
            mp[make_pair(u,v)]=w;
            add(u,v,w);
            add(v,u,w);
            if(w==1){
                Union1(u,v);
                Union2(u,v);
            }
            else if(w==2){
                Union2(u,v);
            }
        }
        dfs(1,0);
        //debug();
        for(int i=1;i<=m;i++){
            int u,v,a,b;
            scanf("%d%d%d%d",&u,&v,&a,&b);
            update(u,v);
        //  debug();
            query1(b,a);
            query2(a);
        }
    }
}

D
不知道怎么回事就猜了个和出题人一样的错误结论就过了。

#include <bits/stdc++.h>
using namespace std;
int a[103],b[103];
bool vis[103];
int n,m;
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++){
            scanf("%d%d",&a[i],&b[i]);
        }
        memset(vis,false,sizeof(vis));
        while(1){
            int id=-1,mx=0;
            for(int i=1;i<=n;i++){
                if(vis[i])continue;
                int res=m/(b[i]+1);
                if(res>mx){
                    mx=res;
                    id=i;
                }
            }
            if(id==-1)break;
            vis[id]=true;
            m=mx;
        }
        int ans=0;
        for(int i=1;i<=n;i++)if(vis[i])ans++;
        printf("%d\n",ans);
    }
}

E
好像找找循环节就行了。

#include<bits/stdc++.h>
using namespace std;
#define ll long long

int t,l,q,xu,xd,yl,yr,A[20];
ll sum=0;

ll M(int x, int y){
    return A[(1ll*(x+y)*(x+y+1)/2+(x))%l];
}

ll getans(int x, int y, int len){
    int lef=len%(2*l);
    ll ans=0, tmp=0;
    for(int i=0; i<lef; i++){
        for(int j=0; j<=i; j++){
            ans+=M(x+j,y+i-j);
        }
    }
    for(int i=0; i<2*l; i++){
        for(int j=0; j<=i+lef; j++){
            tmp+=M(x+j,y+i+lef-j);
        }
    }
    int num=(len/(2*l));
    ans+=tmp*num;
    ans+=4ll*l*sum*(1ll*num*(num-1)/2);
    return ans;
}

int main(){
    scanf("%d",&t);
    while(t--){
        scanf("%d",&l);
        sum=0;
        for(int i=0; i<l; i++){
            scanf("%d",&A[i]);
            sum+=A[i];
        }
        scanf("%d",&q);
        /*
        while(q--){
            int x,y;
            cin>>x>>y;
            cout<<M(x,y)<<endl;
        }
        */

        while(q--){
            scanf("%d%d%d%d",&xu,&yl,&xd,&yr);
            ll ans=getans(xu,yl,(xd-xu)+(yr-yl)+1);
            ans-=getans(xu,yr+1,(xd-xu));
            ans-=getans(xd+1,yl,(yr-yl));
            printf("%lld\n",ans);
        }

    }
    return 0;
}

G
dfs直到前缀不能完全匹配，每次计算贡献。每个节点需要维护当前儿子的左边的有多少个儿子没有dfs过，这个可以用vector版本的bit维护下。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
ll power(ll n,ll p){
    ll ans=1;
    ll base=n;
    while(p){
        if(p&1){
            ans=ans*base%mod;
        }
        base=base*base%mod;
        p>>=1;
    }
    return ans;
}
ll getinv(ll n){
    return power(n,mod-2);
}
int T,n;
vector<int>g[1000003];
vector<int>c[1000003];
int p[1000003];
ll fac[1000003];
ll inv[1000003];
ll ans,tmp;
bool flag;
int clk;
int lowbit(int x){
    return x&(-x);
}
void update(int u,int x){
    int siz=g[u].size();
    while(x<=siz){
        c[u][x]++;
        x+=lowbit(x);
    }
}
int query(int u,int x){
    int siz=g[u].size();
    int ret=0;
    while(x){
        ret+=c[u][x];
        x-=lowbit(x);
    }
    return ret;
}
int bisearch(int u,int x){
    int siz=g[u].size();
    if(g[u][siz-1]<x)return siz;
    int l=0,r=siz-1;
    while(1){
        if(r-l<=1){
            if(g[u][l]>=x)return l;
            return r;
        }
        int mid=(l+r)/2;
        if(g[u][mid]>=x)r=mid;
        else l=mid;
    }
}
void dfs(int u,int fa){
    clk++;
    int siz=g[u].size();
    for(int i=0;i<siz;i++){
        if(i==siz-1&&fa!=0)return;
        int bd=p[clk+1];
        int id=bisearch(u,bd);
        int cnt=id-query(u,id);
        //cout<<u<<' '<<id<<' '<<g[u][id]<<' '<<cnt<<' '<<tmp<<endl;
        if(id==siz||g[u][id]>fa){
            cnt--;
        }
        if(fa==0)cnt++;
        if(fa==0)tmp=tmp*inv[siz-i]%mod;
        else tmp=tmp*inv[siz-i-1]%mod;
        //cout<<siz<<' '<<i<<endl;
        //cout<<tmp<<' '<<cnt<<endl;
        ans=(ans+tmp*cnt)%mod;
        if(id==siz||g[u][id]!=p[clk+1]){
            flag=false;
            return;
        }
        update(u,id+1);
        dfs(g[u][id],u);
        if(!flag)return ;
    }
}
int main(){
    fac[0]=1;
    for(int i=1;i<=1000000;i++)fac[i]=fac[i-1]*i%mod;
    for(int i=1;i<=1000000;i++)inv[i]=getinv(i);
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(int i=1;i<=n;i++){
            scanf("%d",&p[i]);
            g[i].clear();
            c[i].clear();
            c[i].push_back(0);
        }
        ans=0,tmp=1;
        for(int i=1;i<n;i++){
            int u,v;
            scanf("%d%d",&u,&v);
            g[u].push_back(v);
            g[v].push_back(u);
            c[u].push_back(0);
            c[v].push_back(0);
        }
        for(int i=1;i<=n;i++){
            sort(g[i].begin(),g[i].end());
            tmp=tmp*fac[g[i].size()-1]%mod;
        }
        for(int i=1;i<p[1];i++){
            ans=(ans+g[i].size()*tmp)%mod;
        }
        flag=true;
        clk=0;
        tmp=tmp*g[p[1]].size()%mod;
        dfs(p[1],0);
        printf("%lld\n",ans);
    }
}

J
爆搜就完了啊，比赛的时候怎么写挂了啊？

#include <bits/stdc++.h>
using namespace std;
char s[17][17];
int r[4][4][4];
int c[4][4][4];
int ans;
int getnum(char c){
    if(c>='0'&&c<='9')return c-'0';
    return c-'A'+10;
}
void print(int x){
    for(int i=0;i<16;i++)printf("%d",(1<<i)&x?1:0);
    puts("");
}
void rotate(int x,int y){
    int tmp[4];
    for(int i=0;i<4;i++){
        tmp[i]=r[x][y][3-i];
    }
    for(int i=0;i<4;i++){
        r[x][y][i]=c[x][y][i];
    }
    for(int i=0;i<4;i++){
        c[x][y][i]=tmp[i];
    }
}
bool check(int x,int y){
    for(int i=0;i<4;i++){
        int tmp=0;
        for(int j=0;j<x;j++){
            tmp|=c[j][y][i];
            if(tmp&c[j+1][y][i])return false;
        }
    }
    for(int i=0;i<4;i++){
        int tmp=0;
        for(int j=0;j<y;j++){
            tmp|=r[x][j][i];
            if(tmp&r[x][j+1][i])return false;
        }
    }
    return true;
}
void dfs(int flor,int tim){
    if(flor==16){
        ans=min(ans,tim);
    }
    int x=flor/4,y=flor%4;
    for(int i=0;i<4;i++){
        if(check(x,y)){
            dfs(flor+1,tim+i);
        }
        rotate(x,y);
    }
}
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        for(int i=0;i<16;i++){
            scanf("%s",s[i]);
        }
        ans=100;
        memset(r,0,sizeof(r));
        memset(c,0,sizeof(c));
        for(int i=0;i<4;i++){
            for(int j=0;j<4;j++){
                for(int k=0;k<4;k++){
                    for(int l=0;l<4;l++){
                        int x=i*4+k;
                        int y=j*4+l;
                        c[i][j][l]|=(1<<getnum(s[x][y]));
                        r[i][j][k]|=(1<<getnum(s[x][y]));
                    }
                }
            }
        }
        dfs(0,0);
        printf("%d\n",ans);
    }
}

K
分情况讨论下就行了。

#include <bits/stdc++.h>
using namespace std;
char s[1003];
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%s",s);
        int n=strlen(s);
        for(int i=0;i<n;i++){
            if(s[i]=='?'){
                if((i>1&&s[i-1]=='0'&&(s[i-2]<'0'||s[i-1]>'9'))||(i==1&&s[i-1]=='0'))s[i]='+';
                else s[i]='1';
            }
        }
        bool ok=true;
        if(s[0]<'0'||s[0]>'9')ok=false;
        if(s[n-1]<'0'||s[n-1]>'9')ok=false;
        for(int i=0;i<n;i++){
            if(s[i]<'0'||s[i]>'9'){
                if(i==0)continue;
                if(s[i-1]<'0'||s[i-1]>'9')ok=false;
            }
            else {
                if(i==0||s[i-1]<'0'||s[i-1]>'9'){
                    if(s[i]=='0'){
                        if(s[i+1]>='0'&&s[i+1]<='9')ok=false;
                    }
                }
            }
        }
        if(!ok)printf("IMPOSSIBLE\n");
        else printf("%s\n",s);
    }
}

L
签到。

#include <bits/stdc++.h>
using namespace std;
int T,n;
int w[100003];
int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(int i=1;i<=n;i++)scanf("%d",&w[i]);
        int s=w[1];
        int t=w[n];
        printf("%d\n",(int)sqrt(abs(t-s)));
    }
}