雅礼学习10.4

Posted 2021-01-10 kuaileyongheng

雅礼学习10.4

上午考试

各题状况

T1

莫名其妙20分了。

考场上一眼秒掉了这个题目：这不就是个并查集捆绑+快速幂么

然后开开心心这个点如果有这个质因子的话(fa)就指向这个质因子，而每个数字有多个质因子。。。

多个质因子在相互指(fa)的时候指乱套了。。。。

对拍的时候看出来的，然后用(1)个多小时来调这份代码，最后自己都不知道这东西在干嘛了，就凉了。

T2

写了个暴力枚举，期望(20)实际(20)

T3

看到成绩之后：这题怎么会爆(long long)的？？？

然后讲题的时候。。。出题人说第二个部分分里面(a_i)的取值范围给错了，少了三个零。。。

凉凉，(30)挂没了

题目及考场代码

T1

#include <cstdio>
#include <cstring>

inline int read()
{
    int n=0,w=1;register char c=getchar();
    while(c<'0' || c>'9'){if(c=='-')w=-1;c=getchar();}
    while(c>='0' && c<='9')n=n*10+c-'0',c=getchar();
    return n*w;
}
const int N=1e5+1,mod=1e9+7;
int n,tot,a[N],prime[N],head[N*10];
bool notprime[N],vis[N];
struct Edge{
    int v,nxt;
}edge[N<<1];
inline void add(int u,int v)
{edge[++tot]=(Edge){v,head[u]};head[u]=tot;}

inline void Get_prime()
{
    for(int i=2;i<N;++i)
    {
        if(!notprime[i])
            prime[++tot]=i;
        for(int j=1;i*prime[j]<=N && j<=tot;++j)
        {
            notprime[i*prime[j]]=true;
            if(i%prime[j] == 0)
                break;
        }
    }
}

inline void work(int now)
{
    int x=1e5,id=now;
    now=a[now];
    while(now^1)
    {
        ++x;
        if(now%prime[(int)(x-1e5)]==0)
        {
            while(now%prime[x]==0)
                now/=x;
            add(id,x);
        }
    }
}

void dfs(int now)
{
    vis[now]=true;
    for(int i=head[now];i;i=edge[i].nxt)
        dfs(edge[i].v);
}

inline int ksm(long long x,int b)
{
    long long res=1;
    for(;b;b>>=1,x=x*x%mod)
        if(b&1)
            res=res*x%mod;
    return res;
}

int main()
{
    freopen("x.in","r",stdin);
    freopen("x.out","w",stdout);

    int t=read();
    Get_prime();
    tot=0;
    while(t--)
    {
        memset(vis,false,sizeof vis);
        memset(head,0,sizeof head);
        tot=0;
        n=read();
        for(int i=1;i<=n;++i)
        {
            a[i]=read();
            work(i);
        }
        tot=0;
        for(int i=1;i<=n;++i)
            if(!vis[i])
            {
                ++tot;
                dfs(i);
            }
        printf("%d
",(ksm(2,tot)-2+mod)%mod);
    }

    fclose(stdin),fclose(stdout);
    return 0;
}

T2

/*
 * 这啥题啊。。。
 * ao原来d是长度啊。。。。
 *
 * 那不是。。。暴力就很好写了
 * string+map水一发
 */
#include <string>
#include <cstdio>
#include <map>

inline int read()
{
    int n=0,w=1;register char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
    while(c>='0'&&c<='9')n=n*10+c-'0',c=getchar();
    return n*w;
}

const int N=91;
int n,m,d,tot,ans,head[N];
std::string s;
std::map<std::string,int> mp;
struct Edge{
    int v,w,nxt;
}edge[N*(N-1)*2];
inline void add(int u,int v,int w)
{edge[++tot]=(Edge){v,w,head[u]};head[u]=tot;}

void dfs(int now,int step)
{
    if(step==d)
    {
        if(!mp[s])
            mp[s]=++ans;
        return ;
    }
    for(int i=head[now];i;i=edge[i].nxt)
    {
        s.push_back(edge[i].w+'0');
        dfs(edge[i].v,step+1);
        s.erase(s.end()-1);
    }
}

int main()
{
    freopen("y.in","r",stdin);
    freopen("y.out","w",stdout);

    n=read(),m=read(),d=read();
    for(int u,v,w,i=1;i<=m;++i)
    {
        u=read(),v=read(),w=read();
        add(u,v,w),add(v,u,w);
    }
    dfs(1,0);
    printf("%d",ans);

    fclose(stdin);fclose(stdout);
    return 0;
}

T3

/*
 * 我旁边的dalao用一小时切了这个题，然后到处说：哇T3好水我贪心能100啊
 * 然后有人问他：你复杂度多少的
 * 他：100啊
 * 那人：？？？100是个什么复杂度啊，我问的是复杂度
 * 他：aoao。。。我看看啊。。。 O(n*q)的
 * 那人：。。。你再去看看数据范围
 * 他：？？？
 *
 * 笑死我了hhhhh
 * 写了个模拟
 */
#include <cstdio>

inline int read()
{
    int n=0,w=1;register char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
    while(c>='0'&&c<='9')n=n*10+c-'0',c=getchar();
    return n*w;
}

const int N=1e5+1;
int n,q,x[N];

int main()
{
    freopen("z.in","r",stdin);
    freopen("z.out","w",stdout);
    n=read(),q=read();
    for(int i=1;i<=n;++i)
        x[i]=read();
    int len,l,r,ans;
    while(q--)
    {
        len=read();
        ans=l=0,r=len;
        for(int i=1;i<=n;++i)
        {
            if(x[i]>r)
            {
                ans+=x[i]-r;
                l+=x[i]-r;
                r=x[i];
            }
            else
                if(x[i]<l)
                {
                    ans+=l-x[i];
                    r-=l-x[i];
                    l=x[i];
                }
        }
        printf("%d
",ans);
    }

    fclose(stdin);fclose(stdout);
    return 0;
}

正解及代码

T1

对所有(gcd e 1)的两个数字连边，最后统计联通块数量（假设为(cnt)），那么答案为(2^{cnt}-2)

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

const int maxn=1e5+10,maxa=1e6+10,mod=1e9+7;
int t,n,last[maxa],ans;
bool vis[maxn];
vector<int> g[maxn];
int pcnt,prime[maxa],minp[maxa];
bool prm[maxa];

inline void init(){
    for(int i=2;i<maxa;++i){
        if(!prm[i]){
            prime[++pcnt]=i;
            minp[i]=i;
        }
        for(int j=1;j<=pcnt&&i*prime[j]<maxa;++j){
            prm[i*prime[j]]=true;
            minp[i*prime[j]]=prime[j];
            if(i%prime[j]==0)
                break;
        }
    }
}
void dfs(int pos){
    vis[pos]=true;
    for(int i=0;i<g[pos].size();++i)
        if(!vis[g[pos][i]])
            dfs(g[pos][i]);
}

int main(){
    freopen("x.in","r",stdin);
    freopen("x.out","w",stdout);
    init();
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        for(int i=1;i<=pcnt;++i)
            last[prime[i]]=0;
        for(int i=1,x;i<=n;++i){
            vis[i]=false;
            g[i].clear();
            scanf("%d",&x);
            while(x>1){
                int fac=minp[x];
                while(x%fac==0)
                    x/=fac;
                if(last[fac]){
                    g[i].push_back(last[fac]);
                    g[last[fac]].push_back(i);
                }
                last[fac]=i;
            }
        }
        ans=1;
        for(int i=1;i<=n;++i)
            if(!vis[i])
                ans=ans*2%mod,dfs(i);
        printf("%d
",(ans+mod-2)%mod);
    }
    return 0;
}

T2

(f[i][j][mask])表示从(i)出发，到(j)结束，是否存在一条表示为(mask)的路径

但是显然这个复杂度不够优

考虑meet in the middle ，对于每种可能的路径，枚举中间的位置进行判断

(O(2^{frac{d}{2}} imes n imes (n+m))+2^d imes n)

可以用(bitset)进行优化，能让整体除以(64)，由此可以解决(N)更大的问题

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

const int maxn=90+10,maxmask=1<<20/2+1;
int n,m,d,d1,d2,ans;
bitset<maxn> g0[maxn],g1[maxn],dp[maxmask],f[maxmask];

int main(){
    freopen("y.in","r",stdin);
    freopen("y.out","w",stdout);
    scanf("%d%d%d",&n,&m,&d);
    for(int i=1,u,v,c;i<=m;++i){
        scanf("%d%d%d",&u,&v,&c);
        if(c)
            g1[u][v]=g1[v][u]=true;
        else
            g0[u][v]=g0[v][u]=true;
    }
    d2=d/2;d1=d-d2;
    for(int u=n;u;--u){
        for(int i=0;i<maxmask;++i)
            dp[i].reset();
        dp[1][u]=true;
        for(int x=1;x<1<<d1;++x)
            for(int v=1;v<=n;++v)
                if(dp[x][v]){
                    dp[x<<1]|=g0[v];
                    dp[x<<1|1]|=g1[v];
                }
        for(int x=0;x<1<<d1;++x)
            f[x][u]=dp[1<<d1|x].any();
    }
    for(int i=0;i<1<<d1;++i)
        for(int j=0;j<1<<d2;++j)
            if((dp[1<<d2|j]&f[i]).any())
                ++ans;
    printf("%d
",ans);
    return 0;
}

T3

如果存在一个(x_i)使得(x_{i-1}=x_i)或者(x_{i-1}lt x_ilt x_{i+1})，那么可以删掉它（当完成上一个(x)的时候同时完成了这个任务）

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

typedef long long ll;
const int maxn=1e5+10;
int n,m;
ll tot,ans[maxn];
vector<int> x;
vector<pair<int,int> > a;
map<int,int> mp;

inline ll calc(ll k){
    if(!mp.empty()&&mp.begin()->second<0)
        return tot-(mp.size()-1)*k;
    else
        return tot-mp.size()*k;
}
inline void solve(){
    priority_queue<pair<int,int>,vector<pair<int,int> >,greater<pair<int,int> > > q;
    int t=0;
    for(int i=0;i<x.size();++i){
        tot+=abs(x[i]);
        mp[i]=x[i];
        q.push(make_pair(abs(x[i]),i));
    }
    while(!q.empty()){
        int id=q.top().second,tmp=q.top().first;q.pop();
        map<int,int>::iterator p=mp.lower_bound(id);
        if(p==mp.end()||p->first!=id||abs(p->second)!=tmp)
            continue;
        while(t<a.size()&&abs(p->second)>a[t].first)
            ans[a[t].second]=calc(a[t].first),++t;
        if(p!=mp.begin())
            if(p!=prev(mp.end())){
                tmp=p->second,tot-=abs(p->second);
                tmp+=prev(p)->second,tot-=abs(prev(p)->second);
                tmp+=next(p)->second,tot-=abs(next(p)->second);
                mp.erase(prev(p));
                mp.erase(next(p));
                p->second=tmp,tot+=abs(tmp);
                q.push(make_pair(abs(tmp),id));
            }
            else{
                tot-=abs(p->second);
                mp.erase(p);
            }
        else if(p->second>0)
            if(p!=prev(mp.end())){
                tmp=p->second,tot-=abs(p->second);
                tmp+=next(p)->second,tot-=abs(next(p)->second);
                mp.erase(next(p));
                if(tmp){
                    p->second=tmp,tot+=abs(tmp);
                    q.push(make_pair(abs(tmp),id));
                }
                else
                    mp.erase(p);
            }
            else{
                tot-=abs(p->second);
                mp.erase(p);
            }
    }
    while(t<a.size())
        ans[a[t].second]=calc(a[t].first),++t;
}

int main(){
    freopen("z.in","r",stdin);
    freopen("z.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=0,p,last=0;i<n;++i){
        scanf("%d",&p);
        if(p==last)
            continue;
        if(!x.empty()&&(x.back()<0&&p<last||x.back()>0&&p>last))
            x.back()+=p-last;
        else
            x.push_back(p-last);
        last=p;
    }
    for(int i=0,l;i<m;++i){
        scanf("%d",&l);
        a.push_back(make_pair(l,i));
    }
    sort(a.begin(),a.end());
    solve();
    for(int i=0;i<m;++i)
        printf("%lld
",ans[i]);
    return 0;
}

下午讲课：图题

例一

一个(n)个点(m)条边的无向连通图从(1)号点开始(bfs)，可能得到的(bfs)序有很多，取决于出边的访问顺序。现在给出一个(1)到(n)的排列，判断是否可能是一个(bfs)序

(n,mle 2 imes 10^5)

解：

先假设给出的这个(bfs)序列是合法的，令每个节点的权值为在给定序列中的位置

然后从(1)号点开始进行(bfs)，出边的访问顺序按照权值从小到大访问

最后将这样跑(bfs)得到的(bfs)序跟给定的序列进行比较，如果完全一致，那么这个序列是合法的

例二

一个(n)个点(m)条边的无向连通图中每个点都有一个权值，现在要求给每条边定一个权值，满足每个点的权值等于所有相连的边权之和，权值可负。

解：

如果图是一棵树，那么方案就是唯一的，直接判就可以了。
否则随便搞一棵生成树，先不管其他边，跑一遍。
这时根节点可能还不满足条件。
这时考虑其他的边，一条非树边会形成一个环。
如果是一个偶环，那么无论这条非树边怎么变，都不会对根节点产生影响。
而如果是奇环，那么若给这条非树边增加或减少权值，根节点会发生2 的权值变
化。
那么就可以了。

例三

给定一个(n)个点(m)条边的有向带权图，对于一条边权为(w)的边，经过时将获得(w)的收益，之后(w=lfloor frac{w}{2} floor)

问从(1)号点出发最多能获得多少收益

解：

一个强连通分量内部所有的边肯定可以被走到底。
所以缩点后(dp)即可

例四

有(m)个人，(n)张椅子，第(i)个人只能坐在第(u_i)或第(v_i)张椅子上。求有多少种方案满足没有人坐在同一张椅子上。

解：

把椅子当做点，人当做边，变成一个图。
每个连通块可以分开考虑。
假设某个连通块中有(v)个点，(e)条边，由于连通，有(v?1le e)，并且若(egt v)则无解，所以(e)只有(v?1)和(v)两种取值。
假如(e=v?1)，那么该连通块有(v)种方案：考虑枚举每个点不放的情况，其他的点都可以唯一确定。
假如(e=v)且环长(gt 1)，那么该连通块有(2)种方案：考虑环上的一条边，这条边的放法确定后其他的都可以唯一确定。

例五

给定一个(v)个点(e)条边的带权无向图，在图上有(n)个人，第(i)个人位于点(x_i)，一个人通过一条边需要花费这条边的边权的时间。现在每个人可以自由地走。求最短多少时间后满足结束后有人的节点数(ge m)

(n,vle 500)

解：

先(floyd)预处理出两两之间的距离。
然后可以二分答案。二分答案之后，每个人向能走到的点连边。
可以发现合法的条件就是最大匹配数(ge m)。
跑二分图匹配就可以了。