2019/10/27 TZOJ

1001 Gaussian Prime

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=3798

.......................................我是真的一言难尽

1002 Sum of Factorials

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=2696

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define N 1000005
int a[15];
void init()
{
    a[0]=a[1]=1;
    for(int i=2;i<10;i++)
        a[i]=a[i-1]*i;
}
int main()
{
    init();
    int n;
    while(~scanf("%d",&n),n>=0)
    {
        if(n==0)
        {
            printf("NO
");
            continue;
        }
        int all=0;
        for(int i=9;i>=0;i--)
        {
            all+=a[i];
            if(all>n) all-=a[i];
        }
        if(all==n) printf("YES
");
        else printf("NO
");
    }
}

1003 Billboard

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6016

黑板报h*w，第i个广告是1*wi。优先上、左。线段树存剩余最大容量。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define M 200005
int tre[M<<2];
void build(int ind,int left,int right,int w)
{
    tre[ind]=w;
    if(left==right) return;
    int mid=(left+right)>>1;
    build(ind*2,left,mid,w);
    build(ind*2+1,mid+1,right,w);
    tre[ind]=max(tre[ind*2],tre[ind*2+1]);
}
int finda(int ind,int left,int right,int a)
{
    if(left==right)
    {
        tre[ind]-=a;
        return left;
    }
    int mid=(left+right)>>1,ans;
    if(tre[ind*2]>=a) ans=finda(ind*2,left,mid,a);
    else ans=finda(ind*2+1,mid+1,right,a);
    tre[ind]=max(tre[ind*2],tre[ind*2+1]);
    return ans;
}
int main()
{
    int h,w,n;
    while(~scanf("%d%d%d",&h,&w,&n))
    {
        int minn=min(h,n);
        build(1,1,minn,w);
        while(n--)
        {
            int a;scanf("%d",&a);
            if(tre[1]<a) printf("-1
");
            else printf("%d
",finda(1,1,minn,a));
        }
    }
}

1004 Monkey Party

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6054

咕

1005 Happy Necklace

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6064

一点都不happy的项链，f(2)=3, f(3)=4, f(4)=6, f(n)=f(n-1)+f(n-2)，矩阵快速幂。

　　1    1    0

A = 0    0    1

　　1    0    0

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
struct p{
    ll rol,col;
    ll matris[3][3];
};
p Matris(p a,p b)
{
    p tem;
    memset(tem.matris,0,sizeof(tem.matris));
    for(ll i=0;i<3;i++)
        for(ll j=0;j<3;j++)
            for(ll k=0;k<3;k++)
            {
                tem.matris[i][j]+=a.matris[i][k]*(b.matris[k][j]%mod);
                tem.matris[i][j]%=mod;
            }
    return tem;
}
ll quick_mi(ll n)
{
    p res,tep;
    res.matris[0][0]=6,res.matris[0][1]=4,res.matris[0][2]=3;
    for(int i=1;i<=2;i++)
        for(int j=0;j<=2;j++)
            res.matris[i][j]=0;
    tep.matris[0][0]=1,tep.matris[0][1]=1,tep.matris[0][2]=0;
    tep.matris[1][0]=0,tep.matris[1][1]=0,tep.matris[1][2]=1;
    tep.matris[2][0]=1,tep.matris[2][1]=0,tep.matris[2][2]=0;
    while(n)
    {
        if(n&1) res=Matris(res,tep);
        tep=Matris(tep,tep);
        n/=2;
    }
    return res.matris[0][0]%mod;
}
int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
        ll n;scanf("%lld",&n);
        if(n==2) printf("3
");
        else if(n==3) printf("4
");
        else if(n==4) printf("6
");
        else printf("%lld
",quick_mi(n-4));
    }
}

1006 CA Loves GCD

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6068

我不爱

1007 Squarefree number

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6071

欧拉筛把10^6内的素数标下，如果除了两个或以上就No。素数除完了可能还很大，就要看是不是完全平方数，否则n就是一个很大的素数或者两个大于10^6的素数乘积。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define N 1000005
int prim[N+5];
int vis[N+5],cnt=0;
void init()
{
    vis[0]=vis[1]=1;
    for(int i=0;i<N;i++)
        if(!vis[i])
        {
            prim[cnt++]=i;
            for(ll j=(ll)i*i;j<N;j+=i) vis[j]=1;
        }
}
int main()
{
    init();
    int t;scanf("%d",&t);
    for(int kk=1;kk<=t;kk++)
    {
        int flag=1;
        ll n;scanf("%lld",&n);
        for(int i=0;i<cnt;i++)
        {
            if(n%prim[i]==0)
            {
                int num=0;
                while(n%prim[i]==0) n/=prim[i],num++;
                if(num>=2)
                {
                    flag=0;
                    printf("Case %d: No
",kk);
                    break;
                }
            }
        }
        if(n>1000000&&flag)
        {
            int q=(int)sqrt(n);
            if((ll)q*q==n) printf("Case %d: No
",kk);
            else printf("Case %d: Yes
",kk);
            continue;
        }
        if(flag) printf("Case %d: Yes
",kk);
    }
}

1008 Gym Class

http://www.tzcoder.cn/acmhome/problemdetail.do?&method=showdetail&id=6081

选自己前面包括自己的最小id作为评分，记录一下要加上的评分minn。没有要求的入度为0进队列。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAXN=1e9;
const int N=100005;
vector<int> vec[N];
int inde[N];
priority_queue<int> q;
void init()
{
    memset(inde,0,sizeof(inde));
    for(int i=1;i<=N;i++) vec[i].clear();
    while(!q.empty()) q.pop();
}
int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
        ll ans=0,minn=MAXN;
        init();
        int n,m;scanf("%d%d",&n,&m);
        while(m--)
        {
            int a,b;scanf("%d%d",&a,&b);
            inde[b]++;
            vec[a].push_back(b);
        }
        for(int i=1;i<=n;i++)
            if(!inde[i]) q.push(i);
        while(!q.empty())
        {
            ll p=q.top();q.pop();
            minn=min(p,minn);
            ans+=minn;
            for(int i=0;i<vec[p].size();i++)
            {
                inde[vec[p][i]]--;
                if(inde[vec[p][i]]==0) q.push(vec[p][i]);
            }
        }
        printf("%lld
",ans);
    }
}