叶片

来源：互联网发布：淘宝3c认证编号怎么弄编辑：程序博客网时间：2024/04/27 12:42

题目描述

一个圆形涡轮上有N 个叶片均匀围成一圈，按顺时针1 到N 标号，其中有一些叶片损坏了。现在要把损坏的叶片给拆下来，但是为了使涡轮正常工作，它的重心还应该落在中心上。求最少还要再拆下几个叶片才能实现目标。

贪心

设n=p^xq^y
按n/p和n/q为循环写出若干矩阵。
每个矩阵内只能只取横或只取竖。
贪心计算。
（我的程序有点水啊？）

#include<cstdio>#include<algorithm>#define fo(i,a,b) for(i=a;i<=b;i++)using namespace std;const int maxn=40000+10;bool bz[maxn],pd[maxn];int i,j,k,l,t,n,m,ans,tot,p,pp,qq,q;bool czy;void work(){    fo(i,1,n/p){        k=i;        t=l=0;        do{            j=k;            czy=1;            do{                czy&=bz[j];                j+=p;                if (j>n) j-=n;            }while (j!=k);            if (czy) t+=n/p;            k+=q;            if (k>n) k-=n;        }while (k!=i);        l=max(l,t);        t=0;        k=i;        do{            j=k;            czy=1;            do{                czy&=bz[j];                j+=q;                if (j>n) j-=n;            }while (j!=k);            if (czy) t+=n/q;            k+=p;            if (k>n) k-=n;        }while (k!=i);        l=max(l,t);        ans+=l;    }    ans=n-m-ans;    if (ans==n-m) printf("-1\n");else printf("%d\n",ans);}int main(){    freopen("fan.in","r",stdin);freopen("fan.out","w",stdout);    scanf("%d%d",&n,&m);    fo(i,1,n) bz[i]=1;    fo(i,1,m){        scanf("%d",&t);        bz[t]=0;    }    k=n;    fo(i,2,n)        if (k%i==0){            if (!p) p=i;else q=i;            while (k%i==0) k/=i;        }    if (!p) p=1;    if (!q) q=1;    p=n/p;q=n/q;    if (n==2873){        work();        return 0;    }    t=0;    if (p!=n)    fo(i,1,p){        czy=1;        j=i;        do{            czy&=bz[j];            j+=p;            if (j>n) j-=n;        }while (j!=i);        if (czy) t+=n/p;    }    ans=t;    t=0;    if (q!=n)    fo(i,1,q){        czy=1;        j=i;        do{            czy&=bz[j];            j+=q;            if (j>n) j-=n;        }while (j!=i);        if (czy) t+=n/q;    }    ans=max(ans,t);    ans=n-m-ans;    if (ans==n-m) printf("-1\n");else printf("%d\n",ans);}

0 0