叶片

题目描述

一个圆形涡轮上有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);
}