2020ICPC 江西省大学生程序设计竞赛

2020ICPC 江西省大学生程序设计竞赛

A Simple Math Problem

分析：

• 莫比乌斯函数反演

  $$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=1}^i[gcd(i,j)==1]f(j)\\ =& \sum _{j=1}^n\sum _{i=j}^{n}f(j)\sum _{d|(i,j)}u(d)\\ =& \sum _{d=1}^n\sum _{j=1}^{[\frac{n}{d}]}\sum _{i=j}^{[\frac{n}{d}]}f(j\ast d)\ast u[d]\\ =& \sum _{d=1}^{n}u(d)\sum _{j=1}^{[\frac{n}{d}]}f(j\ast d)\sum _{i=j}^{[\frac{n}{d}]}1\\ =& \sum _{d=1}^{n}u(d)\sum _{j=1}^{[\frac{n}{d}]}f(j\ast d)\ast ([\frac{n}{d}]-j+1)\end{array}$$

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

const int N=1e5+5;
int mu[N], b[N], pri[N];
void mus(int n)
{
    int tot=0; 
    mu[1] = 1;
    for(int i=2;i<n;i++)
    {
        if(!b[i]) mu[i] = -1, pri[++tot] = i;
        for(int j=1;j<=tot && i*pri[j]<n;j++)
        {
            b[i*pri[j]] = 1;
            if(i%pri[j]==0) break ;
            mu[i*pri[j]] = -mu[i];
        }
    }
}
int f[N];
signed main()
{
    int n;
    cin>>n;
    mus(n+5);
    for(int i=1;i<=n;i++) f[i]=i%10+f[i/10];
    long long ans=0;
    for(int d=1;d<=n;d++)
    {
        for(int j=1;j*d<=n;j++)
        {
            ans+=mu[d]*f[j*d]*(n/d-j+1);
        }
    }
    cout<<ans<<endl;
    return 0;
}

Apple

分析：

• 先按 $1,2,3,...,n$ 分给 $n$ 个人，多余的分给最后一个人即可

  因此，只需判断 $m$ 是否足够分即可

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

const int N=1e5+5;

signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int T;
    cin>>T;
    while(T--)
    {
        int n,m;
        cin>>n>>m;
        if(n>=m*(m+1)/2) cout<<"possible"<<endl;
        else cout<<"impossible"<<endl;
    }
    return 0;
}

Color Sequence

分析：

• 求满足条件的连续子序列的数量

  • 条件：子序列每种颜色出现的数量都是偶数

• 再看颜色种类数 $<=21$

  会想到要去用到 $DP$，考虑状态的转移，下面只考虑一种颜色来分析：

  $f[i][j]$ 表示前j个颜色，第 $i$ 种颜色的贡献状态

  要求是偶数，$f[i][j]$ 的贡献就等于前 $j$ 个颜色 $c[i]$ 出现的次数 $/2$（奇数只跟奇数，偶数只跟偶数）

• 多种颜色同时考虑

  与单个颜色的区别：状态由单个颜色的情况，转换成多种颜色的共同状态

  每种颜色要么偶数个，要么奇数个，与具体数值没有关系

  又因为 $c_i<=20$ 故可以想到用亦或的形式来维护状态（前缀亦或求取状态）

• 最后跑一遍状态前缀和的累加即可

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

typedef long long ll;
const int N=1e6+5;
int a[N],f[N],s[N<<3];
signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int n;
    cin>>n;
    for(int i=1;i<=n;i++) cin>>a[i];
    ll pre=0;
    for(int i=1;i<=n;i++)
    {
        pre^=(1<<a[i]); // 前缀亦或求状态
        f[i]=pre; // 状态记录
    }
    int ans=0;
    s[0]=1;
    for(int i=1;i<=n;i++)
    {
        ans+=s[f[i]]; // 累加
        s[f[i]]++; // 状态前缀和
    }
    cout<<ans<<endl;
    return 0;
}

Mathematical Practice

分析：

• $m$ 个盒子，$n$ 个不同球，每个球可放进任何一个盒子里或者不放
  故 $ans=(m+1{)}^n$

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

typedef long long ll;
const int mo=998244353;
int ksm(int a,int b,int p=mo)
{
    int ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%p;
        a=a*a%p; b>>=1;
    }
    return ans;
}
signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int n,m;
    cin>>n>>m;
    cout<<ksm(m+1,n)<<endl;
    return 0;
}

Simple Math Problem

分析：

• 分两种情况
  • 在上半部分：直接计算即可，三角形的面积+行数
  • 在下半部分：容斥一下

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

typedef long long ll;
const int N=1e5+5;

signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    ll n,x,y;
    cin>>n>>x>>y;
    ll ans;
    if(x+y+1<=n)
    {
        ans=(x+y+1)*(x+y)/2+x;
        cout<<ans<<endl;
    }
    else
    {
        x=n-x-1;
        y=n-y-1;
        ans=(x+y+1)*(x+y)/2+x;
        ans=n*n-ans-1;
        cout<<ans<<endl;
    }
    return 0;
}

Travel Expense

$Floyd+$ 二分

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

typedef long long ll;
const int N=105;
int a[N],f[N][N];
bool check(int k,int d,int w)
{
    int res=0,s=1;
    for(int i=1;i<=d;i++)
    {
        s*=k;
        res+=s;
        if(res>w) return 0;
    }
    return 1;
}
signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int n,m;
    cin>>n>>m;
    memset(f,0x3f,sizeof(f));
    for(int i=1;i<=m;i++)
    {
        int u,v;
        cin>>u>>v;
        f[u][v]=f[v][u]=1;
    }
    for(int i=1;i<=n;i++) f[i][i]=0;
    for(int k=1;k<=n;k++)
    {
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
            }
        }
    }
    int q;
    cin>>q;
    while(q--)
    {
        int u,v,w;
        cin>>u>>v>>w;
        if(f[u][v]>w) { cout<<"0"<<endl; continue; }
        if(f[u][v]>30) { cout<<"1"<<endl; continue; }
        int ans=0,d=f[u][v];
        if(d==1) {cout<<w<<endl; continue; }
        //cout<<"d="<<d<<endl;
        int l=0,r=w;
        while(l<=r)
        {
            int mid=l+r>>1;
            if(check(mid,d,w)) l=mid+1, ans=mid;
            else r=mid-1;
        }
        cout<<ans<<endl;
    }
    return 0;
}

Zoos’s Animal Codes

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

signed main()
{
    string sr,st;
    cin>>sr>>st;
    cout<<sr+st<<endl;
    return 0;
}