Hello 2019

Hello 2019

Hello 2019

C. Yuhao and a Parenthesis

处理两组数据

1. 一组数据是放左边的括号，且该括号放左边要合法。

2. 一组数据是放右边的括号，且该括号放右边要合法。

合法的意义

令左括号为1，右括号-1，为放左边的括号从左往右的任意的前缀和不低于0

令左括号为-1，右括号1，为放右边的括号从右往左的任意的前缀和不低于0

将两种数据分别放入multiset，进行匹配，计算答案。

#include <iostream>
#include <string>
#include <vector>
#include <set>
using namespace std;
const int N=1e5+7;
int n,sum;
vector <int> ho1,ho2;
multiset<int> sa1,sa2;
string a[N];
int check(string s,int sz){
    int cnt1=0,cnt2=0;
    int ok1=1,ok2=1;
    for(int i=0;i<sz;i++){
        if(s[i]=='(') cnt1++;
        else cnt1--;
        if(cnt1<0) ok1=0;
    }
    for(int i=sz-1;i>=0;i--){
        if(s[i]==')') cnt2++;
        else cnt2--;
        if(cnt2<0) ok2=0;
    }

    if(cnt1>0&&ok1){
        sum=cnt1;
        return 1;
    }else if(cnt1<=0&&ok2){
        sum=cnt1;
        return 1;
    }else{
        return 0;
    }

}
int main(){
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        if(check(a[i],a[i].size())){
            if(sum>0) sa1.insert(sum);
            else sa2.insert(-sum);
        }
    }
    int ans=0;
    for(auto i:sa1){
        auto it=sa2.find(i);
        if(it==sa2.end()) continue;
        sa2.erase(it);
        ans++;
    }
    int cnt=0;
    for(auto i:sa2){
        if(i==0) cnt++;
    }
    cout<<ans+cnt/2<<endl;
}

---

D. Makoto and a Blackboard

结论

该期望满足积性函数，即f(a,b)=f(a)f(b),a,b互质。

那么可以通过唯一分解得到,$x={\prod }_k{p}_k^b$，$p$为质数。

$E(x)={\prod }_{k}E(p_k^b)$

那么只需算质因子$p_k^b$的期望，就变成了普通的用记忆化来算贡献了。

状态划分为，$f(p,s)$为变化$p$次后为$a^s$的期望。

期望公式为：$f(p,s)={\sum }_{i=0}^{i=s}f(p-1,i)\times \frac{1}{i+1}$

#include <iostream>
#include <cstring>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
ll n,k;
ll dp[10007][100],inv[100];
ll ksm(ll x,ll p){
    ll res=1;
    while(p){
        if(p&1) res=res*x%mod;
        p/=2;
        x=x*x%mod;
    }
    return res;
}
//变化p次后为num^s的期望。
ll dfs(ll p,ll s,ll num){
    if(p==0){
        return ksm(num,s);
    }
    //为什么不存num，因为num不参与转移。
    if(dp[p][s]!=-1) return dp[p][s];
    ll ans=0;
    for(int i=0;i<=s;i++){
        ans=(ans+dfs(p-1,i,num)*inv[s+1]%mod)%mod;
    }
    return dp[p][s]=ans;
}
ll DP(ll p,ll s,ll num){
    //初始化
    memset(dp,-1,sizeof(dp));
    
    return dfs(p,s,num);
}
int main(){
    cin>>n>>k;
    memset(dp,-1,sizeof dp);
    //线性处理逆元
    inv[1]=1;
    for (ll i=2; i<=60; i++)
        inv[i]=(mod-mod/i)*inv[mod%i]%mod;
    ll ans=1;
    //简单的唯一分解
    for(ll i=2;i*i<=n;i++){
        ll cnt=0;
        while(n%i==0){
            n/=i;
            cnt++;
        }
        if(cnt==0) continue;
        
        ans=ans*DP(k,cnt,i)%mod;
    }
    
    if(n!=1) ans=ans*DP(k,1,n)%mod;
    cout<<ans<<"\n";
}