洛谷10月月赛R1·浴谷八连测R1·提高组

SAC E#1 -T1 一道中档题 Factorial

题目大意：求n!在k进制下末尾0的个数。
做法：我们对k进行质因数分解，然后用k的质因数去除 n!，然后取num2/num1的最小值就可以了。
看数据范围，对k的质因数分解是要用rho大整数分解来做的，然而数据好像比较弱，筛到1e7就可以过的QAQ。
对于rho的应用见我前边一篇博客 欧拉函数专项

#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#define LL long long
using namespace std;
LL n,k,ans,cnt,m,MOD;
bool np[10000009];
LL prime[10000009],num[10000009],q[10000009],r;
void su()
{
    np[1]=1,np[0]=1;
    for(int i=2;i<=10000001;i++)
    {
        if(!np[i]) prime[++cnt]=i;
        for(int j=1;j<=cnt&&1ll*(prime[j]*i)<=10000001;j++)
        {
            np[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
void pre()
{
    LL p=k;
    for(int i=1;i<=cnt;i++)
    {
        LL x=prime[i],f=0;
        while(p%x==0&&p)
        {
            p/=x;
            num[x]++;
            f=1;
        }
        if(f) q[++r]=prime[i];
        if(!p) break;
    }
}
LL count(LL x)
{
    LL s=0,p=x;
    while(p<=n)
    {
        s+=n/p;
        p*=x;
    }
    return s;
}
int main()
{
    scanf("%lld%lld",&n,&k);

    su();
    pre();
    ans=1e18;
    for(int i=1;i<=r;i++)
    {
        ans=min(ans,count(q[i])/num[q[i]]);
    }

    if(ans==1e18) ans=0;
    printf("%lld",ans);
    return 0;
}

SAC E#1 - T2一道简单题 Sequence2

清晰版题目描述
小强拿到一个3×n的数组，要在每一列选一个数（或者不选），满足以下条件：
1.如果在第一行选，那它必须大于等于上一个数
2.如果在第二行选，那么必须小于等于上一个数
3.如果在第三行选，对于连续的一段在第三行选的数，必须满足方向相同（都小于等于上一个数或者都大于等于上一个数）

20分：暴力枚举每一个位置 不选，选1，选2，选3
60分：用dp[i][j]表示第 i 列选第 j 行达到的最长长度。
设 dp i,j 表示：第 i 个位置选择第 j 个序列的数，前 i 个位置
所能构成的最长子序列的长度。特别的，dp i,2 表示第 i 个位
置选择第三个序列的数，且大于等于前一个数；dp i,3 表示
第 i 个位置选择第三个序列的数，且小于等于前一个数。
dp i,1 = max(dp j,k ),(1≤ j < i,1≤ k ≤ 4,a j,k ≤ a i,1 ) + 1
dp i,2 = max(dp j,k ),(1 ≤ j < i,1 ≤ k ≤ 4,a j,k ≥ a i,2 ) + 1
dp i,3 = max(dp j,k ),(1≤ j < i,1 ≤ k ≤ 4,k≠4,a j,k ≤ a i,3 ) + 1
dp i,4 = max(dp j,k ),(1≤ j < i,1≤ k ≤ 4,k≠3,a j,k ≥ a i,4) + 1
边界条件：dp 0,j = 0。
状态数 O(n)，转移复杂度 O(n)，总复杂度 O(n^2 )。
100分：在60分的基础上，在选取最大值的时候，用树状数组来优化 （需要离散化）。
在这里要用8个树状数组，为了节省代码量，我们只用树状数组来计算前面的最大值 ( 所以需要加上一些特殊处理 ) 例如当需要a j,k>=a i,2时我们可以按照m-a i,2来查询（a j,k>=a i,2 即 m-a j,k <= m-a i,2）。

60分裸的dp

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define M 100009
#define LL long long
using namespace std;
int n,a[M][5],dp[M][5],ans;
int main()
{
    scanf("%d",&n);
    for(int j=1;j<=3;j++)for(int i=1;i<=n;i++) 
    scanf("%d",&a[i][j]);
    for(int i=1;i<=n;i++) a[i][4]=a[i][3];
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<i;j++)
        {
            for(int k=1;k<=4;k++)
            {
                if(a[j][k]<=a[i][1]) dp[i][1]=max(dp[i][1],dp[j][k]+1);
                if(a[j][k]>=a[i][2]) dp[i][2]=max(dp[i][2],dp[j][k]+1);
                if(a[j][k]>=a[i][3]&&k!=4) dp[i][3]=max(dp[i][3],dp[j][k]+1);
                if(a[j][k]<=a[i][4]&&k!=3) dp[i][4]=max(dp[i][4],dp[j][k]+1);
            }
        }
    }
    for(int i=1;i<=4;i++) ans=max(ans,dp[n][i]);
    printf("%d",ans);
    return 0;
}

100分

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define M 100009
#define LL long long
using namespace std;
int n,a[M][5],dp[M][5],ans;
int te[10][M*3],t[M*3],tot,m;
inline int read()
{
    register int x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
int query(int id,int x)
{
    int s=0;
    for(int i=x;i;i-=(-i)&i)
     s=max(te[id][i],s);
    return s;
}
void update(int id,int x,int p)
{
    for(int i=x;i<=m;i+=i&(-i))
     te[id][i]=max(te[id][i],p);
}
int main()
{
    scanf("%d",&n);
    for(int j=1;j<=3;j++)
     for(int i=1;i<=n;i++) 
      a[i][j]=read(),t[++tot]=a[i][j];

    sort(t+1,t+tot+1);
    m=unique(t+1,t+tot+1)-t-1;
    for(int j=1;j<=3;j++)
     for(int i=1;i<=n;i++)
      a[i][j]=lower_bound(t+1,t+m+1,a[i][j])-t;//离散化 

    for(int i=1;i<=n;i++)
    {
        dp[i][1]=max(dp[i][1],query(1,a[i][1])+1);
        dp[i][1]=max(dp[i][1],query(3,a[i][1])+1);
        dp[i][1]=max(dp[i][1],query(5,a[i][1])+1);
        dp[i][1]=max(dp[i][1],query(7,a[i][1])+1);

        dp[i][2]=max(dp[i][2],query(2,m-a[i][2]+1)+1);
        dp[i][2]=max(dp[i][2],query(4,m-a[i][2]+1)+1);
        dp[i][2]=max(dp[i][2],query(6,m-a[i][2]+1)+1);
        dp[i][2]=max(dp[i][2],query(8,m-a[i][2]+1)+1);

        dp[i][3]=max(dp[i][3],query(1,a[i][3])+1);
        dp[i][3]=max(dp[i][3],query(3,a[i][3])+1);
        dp[i][3]=max(dp[i][3],query(5,a[i][3])+1);

        dp[i][4]=max(dp[i][4],query(2,m-a[i][3]+1)+1);
        dp[i][4]=max(dp[i][4],query(4,m-a[i][3]+1)+1);
        dp[i][4]=max(dp[i][4],query(8,m-a[i][3]+1)+1);

        update(1,a[i][1],dp[i][1]);
        update(3,a[i][2],dp[i][2]);
        update(5,a[i][3],dp[i][3]);
        update(7,a[i][3],dp[i][4]);

        update(2,m-a[i][1]+1,dp[i][1]);
        update(4,m-a[i][2]+1,dp[i][2]);
        update(6,m-a[i][3]+1,dp[i][3]);
        update(8,m-a[i][3]+1,dp[i][4]);
    } 
    for(int i=1;i<=4;i++) ans=max(ans,dp[n][i]);
    printf("%d",ans);
    return 0;
}