Codeforces #480 Tutorial

Problem A,B,C:

简单的模拟，注意A中p mod q时对q=0特殊处理（注意范围）

Problem D:

Brief Intro:

给定长度为N的数组A，将A中所有连续子序列分成最少的组，使得每组任意一对数的积均为完全平方数

求最终分成组数为K的子序列个数，K属于[1,N]

Algorithm:

能推出的性质：若P,Q两数积为完全平方数，则任意一个质因子的次幂的奇偶性必然相同

那么想判断P，Q是否满足条件，只要保留每个质因子的次幂为0或1，再判断P,Q是否相同即可

下面只要考虑如何O(N^2)地判断

为了能O（1）判断新加入的数是否已经出现过，需要预处理出每一个数的上一个“自己”出现的位置

由于数的范围过广，使用map记录一个数在检索到k时最后的位置

Code:

#include <bits/stdc++.h>

using namespace std;

const int MAXN=+;
int n,dat[MAXN],res[MAXN],pre[MAXN];

vector<int> prime;
map<int,int> mp;

bool isprime(int x)
{
    int up_limit=sqrt(x);
    for(int i=;i<=up_limit;i++)
        if(x%i==) return false;
    return true;
}

void init()
{
    int up_limit=sqrt(1e8+);
    for(int i=;i<=up_limit;i++)
        if(isprime(i)) prime.push_back(i);
}

int trans(int x)
{
    int up_limit=sqrt(abs(x)),ret=x;
    for(int i=;i<prime.size();i++)
    {
        if(prime[i]>up_limit || ret== || ret==-) break;

        int t=ret,cnt=;
        while(t%prime[i]==) t/=prime[i],cnt++;
        if(cnt%==) ret=t*prime[i];
        else ret=t;
    }
    return ret;
}

int main()
{
    cin >> n;
    for(int i=;i<=n;i++)
        cin >> dat[i];

    init();

    for(int i=;i<=n;i++) //质因数分解
        dat[i]=trans(dat[i]);

    for(int i=;i<=n;i++)  //预处理pre
        if(mp.count(dat[i]))
        {
            pre[i]=mp[dat[i]];
            mp[dat[i]]=i;
        }
        else
        {
            pre[i]=-;
            mp[dat[i]]=i;
        }

    for(int i=;i<=n;i++)
    {
        bool f=true;int cnt=;
        for(int j=i;j<=n;j++)
        {
            if(dat[j]) f=false;
            if(pre[j]<i && dat[j]) cnt++,mp[dat[j]]=true; //O（1）判断
            if(!f) res[cnt]++;
            else res[]++;
        }
    }

    for(int i=;i<=n;i++) cout << res[i] << " ";
    return ;
}

Review:

1、特解：0

在看到数据范围后，总要考虑特解。

除非一段全部为0，否则忽略当前遇到的0

2、积为完全平方数的性质：

我当时只想到了传导性，反而忽略了每个质因子次幂奇偶性相同这一性质

从只考虑奇偶性   到    转化后判断相等的方法值得借鉴

3、求解一串数中不同数的个数的预处理：

求出每一个数前一次出现的位置     常用的预处理方式

Problem E:

一棵树中有N个点，每个点的权值为2^N

要舍去K个点，使得这K个点的权值和最小，且剩下的点连通

Algorithm:

显而易见的贪心策略：

反向求解，寻找n-k个要选的点

由于第n个点的权值 > 1~n-1的权值和，所以从第n个点开始贪心选取即可

为了将复杂度控制在 O(NlogN) ，使用树上倍增查找路径终点

Code:

#include <bits/stdc++.h>

using namespace std;

const int MAXN=1e6+;
vector<int> G[MAXN],res;
int n,k,f[MAXN][],vis[MAXN],dep[MAXN];

inline int read()
{
    char ch;int f=,num;
    while(!isdigit(ch=getchar())) f|=(ch=='-');
    num=ch-'';
    while(isdigit(ch=getchar())) num=num*+ch-'';
    return f?-num:num;
}

void dfs(int cur,int anc) //初始化
{
    dep[cur]=dep[anc]+;f[cur][]=anc;
    for(int i=;i<=;i++) f[cur][i]=f[f[cur][i-]][i-];

    for(int i=;i<G[cur].size();i++)
    {
        int v=G[cur][i];
        if(v==anc) continue;
        dfs(v,cur);
    }
}

int main()
{
    n=read();k=read();
    for(int i=;i<n;i++)
    {
        int x=read(),y=read();
        G[x].push_back(y);G[y].push_back(x);
    }

    dfs(n,);

    memset(vis,,sizeof(vis));
    vis[n]=vis[]=;k=n-k-;
    for(int i=n-;i>=;i--)
    {
        if(vis[i]) continue;

        int t=i;
        for(int j=;j>=;j--) //倍增找路径
            if(!vis[f[t][j]]) t=f[t][j];

        if(dep[i]-dep[t]+<=k)
        {
            k-=(dep[i]-dep[t]+);
            t=i;
            while(!vis[t]) vis[t]=,t=f[t][];
        }
        else res.push_back(i);
    }

    sort(res.begin(),res.end());
    for(int i=;i<res.size();i++) cout << res[i] << " ";
    return ;
}

Review:

1、当正向贪心难以实现时，可以尝试反向贪心

2、当要在树上O(logN)搜寻路径时，使用树上倍增法