[hzxjhs模拟2(T1,T3)][欧拉函数][dp]

3 篇文章 0 订阅
1 篇文章 0 订阅

T1

【解题思路】
一个人的步长为ai,那么他所走过的格子一定是gcd(a[i],n)的倍数(如果不信请画一画),然后许多人的步长肯定会有重复的,所以为了防止重复,我们可以这样想。
枚举n的约数,如果有一个人的步长(这里我们把他当作步长与n的gcd),是这个约数d的倍数,那么这个人所走过的格子数为eluer_phi(n/d),也就是gcd(i,n)==d的个数,,为什么?
假设n=12,
ai:2,3,6
枚举n的约数,2满足,且gcd(ai,n)|2这句话的意思是说gcd(ai,n)能被2整除,也就是说2是gcd(ai,n)的倍数,所以,,,(我小学没有学好,,没有理解到整除。。。),那么所有gcd(i,n)==2点都能到达,那么euler_phi(n/d)
本来gcd(i,n)=d,除去d,那么i/d与n/d是互质的。。。,i/d又小于n/d,也就是euler_phi(n/d)的值。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
using namespace std;
int euler_phi(int n){
    int q=sqrt(n+0.5);
    int ans=n;
    for(int i = 2; i <=q; i++)if(n%i==0){
        ans=ans/i*(i-1);
        while(n%i==0)n/=i;
    }
    if(n>1)ans=ans/n*(n-1);//为质数
    return ans;
}
int gcd(int a,int b){
    if(b==0)return a;
    return gcd(b,a%b);
}
int gc[55],n,m;
int main(){
    freopen("run.in","r",stdin);
    freopen("run.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i = 1; i <=m;i++){
        scanf("%d",&gc[i]);
        gc[i]=gcd(gc[i],n);
    }
    int q=sqrt(n+0.5);
    int ans=0;
    for(int i = 1; i <= q;i++)
        if(n%i==0){
            for(int j = 1; j <= m; j++)
                if((i%gc[j])==0){
                    ans+=euler_phi(n/i);
                    break;//这个约数已经算过了
                }
            if(i*i!=n)
                for(int j = 1; j <=m; j++)
                    if(((n/i)%gc[j])==0){
                        ans+=euler_phi(i);
                        break;
                    }
        }
    printf("%d\n",n-ans);
    return 0;
}

T3

这里写图片描述
【题解】
这里写图片描述
代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
using namespace std;
const int N=3050;
struct ii{
    int to,ne;
    ii(int to = 0,int ne = 0):to(to),ne(ne){ }
}ed[N*2];
int v[N];
int head[N],n,k,f[N][N]; 
void dfs(int u,int fa){
    for(int i =head[u];i;i=ed[i].ne)
        if(ed[i].to!=fa){
            for(int j = 0; j<k;j++)//这里的j从0-k-1,那么下面的j-1就是这里的
                f[ed[i].to][j]=f[u][j]+v[ed[i].to];
            dfs(ed[i].to,u);
            for(int j=1;j<=k;j++)
                f[u][j]=max(f[u][j],f[ed[i].to][j-1]);//f[ed[i].to][j-1]就是他的儿子算出来的值
        }
}
int main(){
    freopen("tree.in","r",stdin);
    freopen("tree.out","w",stdout);
    scanf("%d%d",&n,&k);
    for(int i = 1; i <= n; i++)
        scanf("%d",&v[i]);
    for(int i = 1; i <=  n-1;i++){
        int x,y;
        scanf("%d%d",&x,&y);
        ed[i*2-1]=ii(y,head[x]);
        head[x]=i*2-1;
        ed[i*2]=ii(x,head[y]);
        head[y]=i*2;
    }
    dfs(1,0);
    printf("%d",f[1][k-1]+v[1]);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值