牛客小白月赛13

牛客小白月赛13

---

C.小A买彩票

题意：

数据范围：n<=30

解法：

code：

#include<bits/stdc++.h>
using namespace std;
#define int long long
int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}
int ppow(int a,int b){
    int ans=1;
    while(b){
        if(b&1)ans=ans*a;
        a=a*a;
        b>>=1;
    }
    return ans;
}
const int maxm=35;
int d[maxm][maxm*4];
signed main(){
    int n;cin>>n;
    d[0][0]=1;
    for(int i=1;i<=n;i++){
        for(int j=i;j<=i*4;j++){
            for(int k=1;k<=4;k++){
                d[i][j]+=d[i-1][j-k];
            }
        }
    }
    int tot=ppow(4,n);
    int sum=0;
    for(int i=3*n;i<=4*n;i++){
        sum+=d[n][i];
    }
    int g=gcd(sum,tot);
    sum/=g,tot/=g;
    cout<<sum<<'/'<<tot<<endl;
    return 0;
}

---

I.小A取石子

题意：

数据范围：n,k<=1e5

解法：

这题只是在nim博弈的板题上加了一个可以去掉k个石子的条件
复习一下nim博弈必败的条件：石子数量异或和为0
第一遍遍历计算异或和
第二遍遍历枚举去掉石子的是那堆
只要存在异或和不为0就是必胜

code：

#include<bits/stdc++.h>
using namespace std;
const int maxm=1e5+5;
int a[maxm];
signed main(){
    int n,k;cin>>n>>k;
    int sum=0;
    for(int i=1;i<=n;i++)cin>>a[i],sum^=a[i];
    for(int i=1;i<=n&&!sum;i++){
        if(a[i]>=k){//判断去掉k个的情况
            int temp=sum^a[i]^(a[i]-k);
            if(temp)sum=temp;
        }
    }
    if(sum)puts("YES");
    else puts("NO");
    return 0;
}

---

J.小A的数学题

题意：

答案对1e9+7取模

数据范围：n,m<=1e6

解法：

答案为sigma[d²f(d)]

f(d)的计算：当d=1的时候是O(n)的,当d不为1的时候是O(log(n))的
从1到n枚举d即可，只有当d=1的时候比较慢是O(n)，其他都是O(logn)，总复杂度应该是O(nlogn)

ps：
这题容斥也可以计算出f(d)

code：

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxm=1e6+5;
const int mod=1e9+7;
int notprime[maxm];
int prime[maxm],cnt;
int mob[maxm];
void Minit(){
    mob[1]=1;
    for(int i=2;i<maxm;i++){
        if(!notprime[i]){
            prime[cnt++]=i;
            mob[i]=-1;
        }
        for(int j=0;j<cnt;j++){
            if(prime[j]*i>=maxm)break;
            notprime[prime[j]*i]=1;
            mob[prime[j]*i]=i%prime[j]?-mob[i]:0;
            if(i%prime[j]==0)break;
        }
    }
}
signed main(){
    Minit();
    int n,m;cin>>n>>m;
    int k=min(n,m);
    int ans=0;
    for(int i=1;i<=k;i++){
        for(int j=i;j<=k;j+=i){
            ans+=i*i%mod*mob[j/i]%mod*(n/j)%mod*(m/j)%mod;
            ans%=mod;
        }
    }
    cout<<ans<<endl;
    return 0;
}

---