NOIp2018集训test-9-7(pm) (联考一day1)

又被辉神吊打了。今天不仅被辉神李巨吊打，还给基本上给全班垫底了。

看到T3就知道是十进制快速幂，全机房考试的当时应该就我会，结果我tm没找到递推。

Orz lyc BM直接水过，Orz wys六个for循环出递推，Orz辉神手推n^2递推。

不敢说话，我去背BM板子了。

这样下去大概NOIp继续被辉神李巨吊打，可以退役回去学常规了。

1 餐馆 (restaurant)

完全背包裸题。

//Achen
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define Formylove return 0
#define For(i,a,b) for(int i=(a);i<=(b);i++)
#define Rep(i,a,b) for(int i=(a);i>=(b);i--)
const int N=5007; 
typedef long long LL;
typedef double db;
using namespace std;
int T,n,m,S,t[N],w[N];
LL f[N],ans;

template<typename T>void read(T &x)  {
    char ch=getchar(); x=0; T f=1;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=-1,ch=getchar();
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; x*=f;
}

#define ANS
int main() {
#ifdef ANS
    freopen("restaurant.in","r",stdin);
    freopen("restaurant.out","w",stdout);
#endif
    read(T);
    while(T--) {
        memset(f,0,sizeof(f));
        read(n); read(m); read(S);
        For(i,1,n) {
            read(t[i]); read(w[i]);
        }
        For(i,1,m) { read(t[n+i]); read(w[n+i]); } 
        For(i,1,m) For(j,1,n) {
            int x; read(x);
            if(x>0) t[n+i]+=x*t[j];
        }
        For(i,1,n+m) For(j,t[i],S) 
            f[j]=max(f[j],f[j-t[i]]+w[i]);
        ans=0;
        For(j,0,S) ans=max(ans,f[j]);
        printf("%lld
",ans);
    }
    Formylove;
}

2 烯烃 (olefin)

树形dp裸题，最最基础的换根操作。

//Achen
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define Formylove return 0
#define For(i,a,b) for(int i=(a);i<=(b);i++)
#define Rep(i,a,b) for(int i=(a);i>=(b);i--)
const int N=1e5+7;
typedef long long LL;
typedef double db;
using namespace std;
int id,T,n,m,tag[N],ans[N];

template<typename T>void read(T &x)  {
    char ch=getchar(); x=0; T f=1;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=-1,ch=getchar();
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; x*=f;
}

int ecnt,fir[N],nxt[N<<1],to[N<<1],fa[N];
void add(int u,int v) {
    nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v;
}

int f[N];
void dfs(int x) {
    f[x]=0;
    for(int i=fir[x];i;i=nxt[i]) {
        dfs(to[i]);
        f[x]=max(f[x],f[to[i]]+1);
    }
}

void get_max(int &a,int b) { if(b>a) a=b; }

void dfs2(int x,int ff) {
    int se=0,fi=x==1?0:ff+1;
    for(int i=fir[x];i;i=nxt[i]) {
        se=max(se,f[to[i]]+1);
        if(se>fi) swap(fi,se); 
    }
    for(int i=fir[x];i;i=nxt[i]) if(tag[to[i]]) {
        if(fi==f[to[i]]+1) get_max(ans[tag[to[i]]],se+f[to[i]]+1);
        else get_max(ans[tag[to[i]]],fi+f[to[i]]+1);
    }
    for(int i=fir[x];i;i=nxt[i]) {
        if(f[to[i]]+1==fi) dfs2(to[i],se);
        else dfs2(to[i],fi);
    }
}

void init() {
    ecnt=0;
    memset(ans,0,sizeof(ans));
    memset(fir,0,sizeof(fir));
    memset(tag,0,sizeof(tag));
}

#define ANS
int main() {
#ifdef ANS
    freopen("olefin.in","r",stdin);
    freopen("olefin.out","w",stdout);
#endif
    read(id);
    read(T);
    while(T--) {
        init();
        read(n); read(m);
        For(i,2,n) {
            read(fa[i]);
            add(fa[i],i);
        }
        For(i,1,m) {
            int x;
            read(x);
            tag[x]=i;
        }
        dfs(1);
        dfs2(1,0);
        For(i,1,m-1) printf("%d ",ans[i]);
        if(m) printf("%d
",ans[m]);
    }
    Formylove;
}

3 三米诺 (tromino)

十进制快速幂，我不会找递推，只能用背BM板子然后水过去。

//Achen
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define Formylove return 0
#define For(i,a,b) for(int i=(a);i<=(b);i++)
#define Rep(i,a,b) for(int i=(a);i>=(b);i--)
const int N=40007,p=998244353;
typedef long long LL;
typedef double db;
using namespace std;
char s[N];
LL n,f[15]={0,1,3,10,23,62,170,441,1173,3127,8266};

template<typename T>void read(T &x)  {
    char ch=getchar(); x=0; T f=1;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=-1,ch=getchar();
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; x*=f;
}

struct jz {
    LL a[6][6];
    friend jz operator *(const jz&A,const jz&B) {
        jz rs;
        For(i,0,5) For(j,0,5) {
            rs.a[i][j]=0;
            For(k,0,5) 
                (rs.a[i][j]+=A.a[i][k]*B.a[k][j]%p)%=p;
        }
        return rs;
    }
}bs,rs;

void ksm() {
    int len=strlen(s);
    int ok=0;
    Rep(i,len-1,0) {
        int c;
        c=s[i]-'0';
        if(!ok) {
            c--; 
            if(c<0) c+=10;
            else ok=1;
        }
        For(i,1,c) rs=rs*bs;
        jz pr=bs;
        For(i,1,9) bs=bs*pr;
    }  
}

#define ANS
int main() {
#ifdef ANS
    freopen("tromino.in","r",stdin);
    freopen("tromino.out","w",stdout);
#endif
    scanf("%s",s);
    For(i,0,5) For(j,0,5) rs.a[i][j]=bs.a[i][j]=0;
    For(i,0,5) {
        rs.a[i][i]=1;
        if(i<5) bs.a[i+1][i]=1;
    }
    bs.a[5][5]=1;
    bs.a[4][5]=2;
    bs.a[3][5]=6;
    bs.a[2][5]=1;
    bs.a[0][5]=p-1;
    ksm();
    LL ans=0;
    For(i,0,5) 
        (ans+=rs.a[i][0]*f[i+1]%p)%=p;
    printf("%lld
",ans);
    Formylove;
}