Hello 2019
C. Yuhao and a Parenthesis
处理两组数据
-
一组数据是放左边的括号,且该括号放左边要合法。
-
一组数据是放右边的括号,且该括号放右边要合法。
合法的意义
令左括号为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 = ∏ k p k b x=\prod_kp_k^b x=∏kpkb, p p p为质数。
E ( x ) = ∏ k E ( p k b ) E(x)=\prod_k E(p_k^b) E(x)=∏kE(pkb)
那么只需算质因子 p k b p_k^b pkb的期望,就变成了普通的用记忆化来算贡献了。
状态划分为, f ( p , s ) f(p,s) f(p,s)为变化 p p p次后为 a s a^s as的期望。
期望公式为: f ( p , s ) = ∑ i = 0 i = s f ( p − 1 , i ) × 1 i + 1 f(p,s)=\sum_{i=0}^{i=s} f(p-1,i)\times\frac 1{i+1} f(p,s)=∑i=0i=sf(p−1,i)×i+11
#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";
}