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[];
 void init()
 {
     a[]=a[]=;
     for(int i=;i<;i++)
         a[i]=a[i-]*i;
 }
 int main()
 {
     init();
     int n;
     while(~scanf("%d",&n),n>=)
     {
         if(n==)
         {
             printf("NO\n");
             continue;
         }
         int all=;
         for(int i=;i>=;i--)
         {
             all+=a[i];
             if(all>n) all-=a[i];
         }
         if(all==n) printf("YES\n");
         else printf("NO\n");
     }
 }

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<<];
 void build(int ind,int left,int right,int w)
 {
     tre[ind]=w;
     if(left==right) return;
     int mid=(left+right)>>;
     build(ind*,left,mid,w);
     build(ind*+,mid+,right,w);
     tre[ind]=max(tre[ind*],tre[ind*+]);
 }
 int finda(int ind,int left,int right,int a)
 {
     if(left==right)
     {
         tre[ind]-=a;
         return left;
     }
     int mid=(left+right)>>,ans;
     if(tre[ind*]>=a) ans=finda(ind*,left,mid,a);
     else ans=finda(ind*+,mid+,right,a);
     tre[ind]=max(tre[ind*],tre[ind*+]);
     return ans;
 }
 int main()
 {
     int h,w,n;
     while(~scanf("%d%d%d",&h,&w,&n))
     {
         int minn=min(h,n);
         build(,,minn,w);
         while(n--)
         {
             int a;scanf("%d",&a);
             if(tre[]<a) printf("-1\n");
             else printf("%d\n",finda(,,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+;
 struct p{
     ll rol,col;
     ll matris[][];
 };
 p Matris(p a,p b)
 {
     p tem;
     memset(tem.matris,,sizeof(tem.matris));
     for(ll i=;i<;i++)
         for(ll j=;j<;j++)
             for(ll k=;k<;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[][]=,res.matris[][]=,res.matris[][]=;
     for(int i=;i<=;i++)
         for(int j=;j<=;j++)
             res.matris[i][j]=;
     tep.matris[][]=,tep.matris[][]=,tep.matris[][]=;
     tep.matris[][]=,tep.matris[][]=,tep.matris[][]=;
     tep.matris[][]=,tep.matris[][]=,tep.matris[][]=;
     while(n)
     {
         if(n&) res=Matris(res,tep);
         tep=Matris(tep,tep);
         n/=;
     }
     return res.matris[][]%mod;
 }
 int main()
 {
     int t;scanf("%d",&t);
     while(t--)
     {
         ll n;scanf("%lld",&n);
         if(n==) printf("3\n");
         else if(n==) printf("4\n");
         else if(n==) printf("6\n");
         else printf("%lld\n",quick_mi(n-));
     }
 }

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+];
 int vis[N+],cnt=;
 void init()
 {
     vis[]=vis[]=;
     for(int i=;i<N;i++)
         if(!vis[i])
         {
             prim[cnt++]=i;
             for(ll j=(ll)i*i;j<N;j+=i) vis[j]=;
         }
 }
 int main()
 {
     init();
     int t;scanf("%d",&t);
     for(int kk=;kk<=t;kk++)
     {
         int flag=;
         ll n;scanf("%lld",&n);
         for(int i=;i<cnt;i++)
         {
             if(n%prim[i]==)
             {
                 int num=;
                 while(n%prim[i]==) n/=prim[i],num++;
                 if(num>=)
                 {
                     flag=;
                     printf("Case %d: No\n",kk);
                     break;
                 }
             }
         }
         if(n>&&flag)
         {
             int q=(int)sqrt(n);
             if((ll)q*q==n) printf("Case %d: No\n",kk);
             else printf("Case %d: Yes\n",kk);
             continue;
         }
         if(flag) printf("Case %d: Yes\n",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=;
 vector<int> vec[N];
 int inde[N];
 priority_queue<int> q;
 void init()
 {
     memset(inde,,sizeof(inde));
     for(int i=;i<=N;i++) vec[i].clear();
     while(!q.empty()) q.pop();
 }
 int main()
 {
     int t;scanf("%d",&t);
     while(t--)
     {
         ll ans=,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=;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=;i<vec[p].size();i++)
             {
                 inde[vec[p][i]]--;
                 if(inde[vec[p][i]]==) q.push(vec[p][i]);
             }
         }
         printf("%lld\n",ans);
     }
 }