纪中2018暑假培训day3提高a组改题记录（混有部分b组）

day3

模拟赛，看了看a组题，发现是博弈论，非常开心（因为好玩），于是做的a组。结果差点爆零，死命纠结t1的sg函数，但其实只是一个dp，不用扯到sg函数的那种。

t1：

Description

被污染的灰灰草原上有羊和狼。有N只动物围成一圈，每只动物是羊或狼。
该游戏从其中的一只动物开始，报出[1,K]区间的整数，若上一只动物报出的数是x，下一只动物可以报[x+1,x+K]区间的整数，游戏按顺时针方向进行。每只动物报的数字都不能超过M。若一只动物报了M这个数，它所在的种族就输了。问以第i只动物为游戏的开始，最后哪种动物会赢？

Input

第一行输入三个正整数N，M，K。
接下来一行N个正整数，分别表示N只动物的种类，以顺时针的方向给出。0代表羊，1代表狼。

Output

一行输出N个整数，表示若从第i只动物开始，赢的动物的种类。同上，0代表羊，1代表狼。

就是dp,f[i][j]表示第i个动物报数最大为j是否必胜。然后加一个后缀和（应该能叫后缀和吧）（或者说区间和？）优化。

中间还wa了一次，因为在处理s数组时减去的那个f数组越界了，需要特判一下。

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

using namespace std;

int n,m,kk;

int f[10010][5010],a[10010],s[10010][5010];

int main()

{

    freopen("vode.in","r",stdin);

    freopen("vode.out","w",stdout);

    scanf("%d%d%d",&n,&m,&kk);

    for(int i=1;i<=n;i++)

        scanf("%d",&a[i]);

    int p=1;

    for(int i=n+1;i<=n+m;i++)

    {

        if(p>n)

            p=1;

        a[i]=a[p];

        p++;

    }

    for(int i=n+m-1;i>=1;i--)

        for(int j=m-1;j>=1;j--)

        {

                if(a[i]==a[i+1])

                {

                    s[i][j]=s[i][j+1];

                    if(j+kk<m)

                        s[i][j]-=f[i][j+kk];

                    if(s[i+1][j+1]>0)

                    {

                        f[i][j]=1;

                        s[i][j]++;

                    }

                }

                else

                {

                    s[i][j]=s[i][j+1];

                    if(j+kk<m)

                        s[i][j]-=f[i][j+kk];

                    if(s[i+1][j+1]==0)

                    {

                        f[i][j]=1;

                        s[i][j]++;

                    }

                }

        }

    for(int i=1;i<=n;i++)

        for(int j=1;j<=kk;j++)

            if(f[i][j])

            {

                printf("%d ",a[i]);

                break;

            }

            else if(j==kk)

            {

                printf("%d ",a[i]^1);

            }

    fclose(stdin);

    fclose(stdout);

    return 0;

}

不想改a组了，先写了一个b组t3玩。

Description

  话说，小X是个数学大佬，他喜欢做数学题。有一天，小X想考一考小Y。他问了小Y一道数学题。题目如下：
      对于一个正整数N，存在一个正整数T（0<T<N），使得
aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAAAD4AAABECAYAAADZeIbjAAADRUlEQVR4nO2aMXKjMBSGf3ZyFG8KT04AJ7DduHKbTpROs11Kd2lwmXRpt6IJnMA+gceF4S5ayQYvBqEgGYF24ZuRCzPw+JH03pP0HMrAAPnR9wv0xSh8aNgnPPbhOA782KwZi4Sn2HoOnPlHJ9YsEj7BekdBkwBuB9YsEt4to/ChMQofGqPwoVEVnm7hscyJZ0+Ot2VpRfmyd7mWNdMZVj0x/MJ7iJuHbVlADq0jIny5SklUc80NaFJ7sx4RwdnmtclslN4hCVx2D6F/XzeiRHJ/rfCIuNR1xca5EeEH6ZCIFEUmNKi8a0SDoL5raoSzBwVR1usuvb2fGykatQHWu2yEuBKhZcTOLf3CET+B2RIEe7y8FSdywq5N+VV7iEPwpc3T46TxLWLhCTBd8IfM8IuNIXyEuEpnRg6rBZqbME96OrBfguWs+T1C4XF4RP7xJutX9sgPbDL3mJ6A1cIq2fj6vWczUm0UCoSnON0M5RmWbALtX95Yr3MjQN2Iiv3vwosDrza+6MKmHtetOgorsz4JKCk7Cfafew5tzIm07M4BNG5CZGFXQrXHk2O1SycLrPhUn89xmLbr1uglsjRqInTmd264gCRUCUNb31zCGDRG4VX4JfPJh5VIIPsopP1sTY8sYSlNBRX9Dv+5Y6T+s4yrMzvItpiv4c+HqcWfRcK56GfgM/fi0TlxmpsSb8j7qBMFVYea5Q8qi4+m2NPjszXW5dRr8ognQ+bsEV6La2RtYHc44weI4RL0XTUt+56H1p/YGszZbYBo175ojrXCY597+B3MyLZUON/JDZc7vBtc9lvn3LjoZ3yiOK3TrV+/TaxL6wHyDm4XSg23mTWxRnhlT73QTCQwdoczg1g3x7vi7NX5Suh/pzywH0R/DoHBDvVR+NAYhTei82qJO6seZGilPV1VS9xZ9SBDa6jH4QEuP1LaCHr9dAB5XbdyjByHfD2eP0t8KjrVPbJW/1Z9VUuoVz3IUO/xvqolNKoeZKgL76laQvtUtAZl4f1US+hVPchQFN5XtYRm1YMMJY/QcbXEFc2qBxlqPd5xtURO2/P7TPNv1Fe1hH7Vg4xGwvuplri/6kHGuOc2NEbhQ2MUPjRG4UPjD4ped/B8yB0SAAAAAElFTkSuQmCC" alt=" " />
的值是正整数。
      小X给出N，让小Y给出所有可能的T。如果小Y不回答这个神奇的大佬的简单数学题,他学神的形象就会支离破碎。所以小Y求你帮他回答小X的问题。

Input

一个整数N。

Output

第一个数M，表示对于正整数N，存在M个不同的正整数T，使得
aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAAAD4AAABECAYAAADZeIbjAAADRUlEQVR4nO2aMXKjMBSGf3ZyFG8KT04AJ7DduHKbTpROs11Kd2lwmXRpt6IJnMA+gceF4S5ayQYvBqEgGYF24ZuRCzPw+JH03pP0HMrAAPnR9wv0xSh8aNgnPPbhOA782KwZi4Sn2HoOnPlHJ9YsEj7BekdBkwBuB9YsEt4to/ChMQofGqPwoVEVnm7hscyJZ0+Ot2VpRfmyd7mWNdMZVj0x/MJ7iJuHbVlADq0jIny5SklUc80NaFJ7sx4RwdnmtclslN4hCVx2D6F/XzeiRHJ/rfCIuNR1xca5EeEH6ZCIFEUmNKi8a0SDoL5raoSzBwVR1usuvb2fGykatQHWu2yEuBKhZcTOLf3CET+B2RIEe7y8FSdywq5N+VV7iEPwpc3T46TxLWLhCTBd8IfM8IuNIXyEuEpnRg6rBZqbME96OrBfguWs+T1C4XF4RP7xJutX9sgPbDL3mJ6A1cIq2fj6vWczUm0UCoSnON0M5RmWbALtX95Yr3MjQN2Iiv3vwosDrza+6MKmHtetOgorsz4JKCk7Cfafew5tzIm07M4BNG5CZGFXQrXHk2O1SycLrPhUn89xmLbr1uglsjRqInTmd264gCRUCUNb31zCGDRG4VX4JfPJh5VIIPsopP1sTY8sYSlNBRX9Dv+5Y6T+s4yrMzvItpiv4c+HqcWfRcK56GfgM/fi0TlxmpsSb8j7qBMFVYea5Q8qi4+m2NPjszXW5dRr8ognQ+bsEV6La2RtYHc44weI4RL0XTUt+56H1p/YGszZbYBo175ojrXCY597+B3MyLZUON/JDZc7vBtc9lvn3LjoZ3yiOK3TrV+/TaxL6wHyDm4XSg23mTWxRnhlT73QTCQwdoczg1g3x7vi7NX5Suh/pzywH0R/DoHBDvVR+NAYhTei82qJO6seZGilPV1VS9xZ9SBDa6jH4QEuP1LaCHr9dAB5XbdyjByHfD2eP0t8KjrVPbJW/1Z9VUuoVz3IUO/xvqolNKoeZKgL76laQvtUtAZl4f1US+hVPchQFN5XtYRm1YMMJY/QcbXEFc2qBxlqPd5xtURO2/P7TPNv1Fe1hH7Vg4xGwvuplri/6kHGuOc2NEbhQ2MUPjRG4UPjD4ped/B8yB0SAAAAAElFTkSuQmCC" alt=" " />
是整数。
后面是M个数，每一个数代表可能的正整数T（按从小到大的顺序排列）。

就是把式子化简一下，设比值为k，然后会发现T=(2K-2)/(2K-1)*N,我们要让t为正整数，同时（2K-2）/(2K-1)又不能为一，2k-1就一定为n的因子。k为正整数，所以2k-1为奇数，然后就可以n½枚举找n的奇因子了。注意一下2k-2不能为零，也就是2k-1不能为1。

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<cmath>

using namespace std;

long long n,a[1000010],cnt;

int main()

{

    freopen("math.in","r",stdin);

    freopen("math.out","w",stdout);

    scanf("%lld",&n);

    long long sn=sqrt(n);

    for(long long i=1;i<=sn;i++)

        if(!(n%i))

        {

            if((i%2)&&i!=1)

                a[++cnt]=(n/i)*(i-1);

            if((n/i)%2&&(n/i)!=1)

                a[++cnt]=i*(n/i-1);

        }

    cout<<cnt<<" ";

    sort(a+1,a+1+cnt);

    for(int i=1;i<=cnt;i++)

        printf("%lld ",a[i]);

    fclose(stdin);

    fclose(stdout);

    return 0;

}

然后我终于改了a组t2：

Description

有一副n*m的地图，有n*m块地，每块是下列四种中的一种：
墙：用#表示，墙有4个面，分别是前面，后面，左面，右面。
起点：用C表示，为主角的起点，是一片空地。
终点：用F表示，为主角的目的地，是一片空地。
空地：用 . 表示。
其中除了墙不能穿过，其他地方都能走。

主角有以下3种操作：
1.移动到相邻的前后左右的地方，花费一个单位时间。
2.向前后左右其中一个方向发射子弹，子弹沿直线穿过，打在最近的一堵墙的一面，然后墙的这面就会形成一个开口通往秘密通道。同一时间最多只能有两个开口，若出现有3个开口，出现时间最早的开口会立即消失。该操作不用时间。
3.可以从一个与开口相邻的空地跳进去，进入秘密通道，从另外一个开口正对的空地跳出来。这个过程花费一个单位时间。

地图四周都是墙，问主角最少用多少时间从C走到F。C和F
只会出现一次。

Input

第一行输入两个正整数n，m。
接下来n行，每行m个字符描述地图。

Output

输出1个整数，表示最短时间完成路途。如果无解输出nemoguce

我当时打了bfs，然后愉悦爆零。哇的一下哭出声。

这个题是把每一个非墙点与相邻点连边权为1的边，然后找到四个方向的墙，连一条边权为min（dis（此点，墙边传送门的点））+1的边（其实和最近的传送门点应该连边权为dis的边，但可以从边权为1的边走过去，不影响答案）。然后跑spfa，求起点到终点最短路即可。

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<queue>

using namespace std;

int n,m,s,t,cnt,head[250010],dis[250010],vis[250010];

char mp[510][510];

struct Edge

{

    int v,nxt,val;

}e[10000010];

void add(int u,int v,int val)

{

    e[++cnt].v=v;

    e[cnt].nxt=head[u];

    e[cnt].val=val;

    head[u]=cnt;

}

int num(int x,int y)

{

    return (x-1)*m+y;

}

void spfa()

{

    memset(dis,0x3f3f3f3f,sizeof(dis));

    queue<int>q;

    dis[s]=0;

    vis[s]=1;

    q.push(s);

    while(!q.empty())

    {

        int u=q.front();

        q.pop();

        vis[u]=0;

        for(int i=head[u];i;i=e[i].nxt)

        {

            int v=e[i].v;

            if(dis[v]>dis[u]+e[i].val)

            {

                dis[v]=dis[u]+e[i].val;

                if(!vis[v])

                {

                    q.push(v);

                    vis[v]=1;

                }

            }

        }

    }

}

int main()

{

    freopen("portal.in","r",stdin);

    freopen("portal.out","w",stdout);

    cin>>n>>m;

    for(int i=1;i<=n;i++)

        for(int j=1;j<=m;j++)

        {

            cin>>mp[i][j];

            if(mp[i][j]=='C')

                s=num(i,j);

            else if(mp[i][j]=='F')

                t=num(i,j);

        }

    for(int i=2;i<n;i++)

        for(int j=2;j<m;j++)

            if(mp[i][j]!='#')

            {

                int x1=i,x2=i,y1=j,y2=j;

                int mindis=0x3f3f3f3f;

                while(mp[x1+1][j]!='#')

                    x1++;

                while(mp[x2-1][j]!='#')

                    x2--;

                while(mp[i][y1+1]!='#')

                    y1++;

                while(mp[i][y2-1]!='#')

                    y2--;

                mindis=min(x1-i,min(i-x2,min(y1-j,j-y2)));

                mindis++;

                add(num(i,j),num(x1,j),mindis);

                add(num(i,j),num(x2,j),mindis);

                add(num(i,j),num(i,y1),mindis);

                add(num(i,j),num(i,y2),mindis);

                if(mp[i+1][j]!='#')

                    add(num(i,j),num(i+1,j),1);

                if(mp[i-1][j]!='#')

                    add(num(i,j),num(i-1,j),1);

                if(mp[i][j+1]!='#')

                    add(num(i,j),num(i,j+1),1);

                if(mp[i][j-1]!='#')

                    add(num(i,j),num(i,j-1),1);

            }

    spfa();

    if(dis[t]==0x3f3f3f3f)

        cout<<"nemoguce";

    else

        cout<<dis[t];

    fclose(stdin);

    fclose(stdout);

    return 0;

}

其实题目不是不可做，就是考试死活想不到正路上......dalao炒鸡多啊我的天，orz。

然而我还是改完了t3：

Description

有n个城市，标号为1到n，修建道路花费m天，第i天时，若gcd(a,b)=m-i+1，则标号为a的城市和标号为b的城市会建好一条直接相连的道路，有多次询问，每次询问某两座城市最早什么时候能连通。

Input

第一行输入三个正整数n,m,q，其中q表示询问个数。
接下来q行，每行两个正整数x，y，表示询问城市x和城市y最早什么时候连通。

Output

输出q行，每行一个正整数，表示最早连通的天数

并查集，按秩合并，把时间赋为点权（我之前纠结了半天怎么连边权2333），在并查集树上找max点权就行。

re了好几次，最后发现freopen文件名打错了......

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

using namespace std;

int n,m,q,fa[100010],dep[100010],siz[100010],tim[100010],x,y;

int find_fa(int x)

{

    if(fa[x]==x)

        return x;

    return find_fa(fa[x]);

}

void update(int x)

{

    if(fa[x]==x)

        return;

    update(fa[x]);

    dep[x]=dep[fa[x]]+1;

}

int query(int a,int b)

{

    int res=0;

    if(dep[a]<dep[b])

        swap(a,b);

    while(dep[a]>dep[b])

    {

        res=max(res,tim[a]);

        a=fa[a];

    }

    while(a!=b)

    {

        res=max(res,max(tim[a],tim[b]));

        a=fa[a];

        b=fa[b];

    }

    return res;

}

int main()

{

    freopen("pictionary.in","r",stdin);

    freopen("pictionary.out","w",stdout);

    scanf("%d%d%d",&n,&m,&q);

    for(int i=1;i<=n;i++)

        fa[i]=i,siz[i]=1;

    for(int i=1;i<=m;i++)

    {

        int r=n/(m-i+1);

        int p=m-i+1;

        for(int j=2;j<=r;j++)

        {

            int f1=find_fa(p*(j-1));

            int f2=find_fa(p*j);

            if(f1!=f2)

            {

                if(siz[f1]<siz[f2])

                    swap(f1,f2);

                fa[f2]=f1;

                tim[f2]=i;

                siz[f1]=max(siz[f1],siz[f2]+1);

            }

        }

    }

    for(int i=1;i<=n;i++)

        update(i);

    for(int i=1;i<=q;i++)

    {

        scanf("%d%d",&x,&y);

        printf("%d\n",query(x,y));

    }

    fclose(stdin);

    fclose(stdout);

    return 0;

}