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=;
 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=; T f=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; 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,,sizeof(f));
         read(n); read(m); read(S);
         For(i,,n) {
             read(t[i]); read(w[i]);
         }
         For(i,,m) { read(t[n+i]); read(w[n+i]); }
         For(i,,m) For(j,,n) {
             int x; read(x);
             if(x>) t[n+i]+=x*t[j];
         }
         For(i,,n+m) For(j,t[i],S)
             f[j]=max(f[j],f[j-t[i]]+w[i]);
         ans=;
         For(j,,S) ans=max(ans,f[j]);
         printf("%lld\n",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+;
 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=; T f=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
 }

 int ecnt,fir[N],nxt[N<<],to[N<<],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]=;
     for(int i=fir[x];i;i=nxt[i]) {
         dfs(to[i]);
         f[x]=max(f[x],f[to[i]]+);
     }
 }

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

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

 void init() {
     ecnt=;
     memset(ans,,sizeof(ans));
     memset(fir,,sizeof(fir));
     memset(tag,,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,,n) {
             read(fa[i]);
             add(fa[i],i);
         }
         For(i,,m) {
             int x;
             read(x);
             tag[x]=i;
         }
         dfs();
         dfs2(,);
         For(i,,m-) printf("%d ",ans[i]);
         if(m) printf("%d\n",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=,p=;
 typedef long long LL;
 typedef double db;
 using namespace std;
 char s[N];
 LL n,f[]={,,,,,,,,,,};

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

 struct jz {
     LL a[][];
     friend jz operator *(const jz&A,const jz&B) {
         jz rs;
         For(i,,) For(j,,) {
             rs.a[i][j]=;
             For(k,,)
                 (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=;
     Rep(i,len-,) {
         int c;
         c=s[i]-'';
         if(!ok) {
             c--;
             if(c<) c+=;
             else ok=;
         }
         For(i,,c) rs=rs*bs;
         jz pr=bs;
         For(i,,) 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,,) For(j,,) rs.a[i][j]=bs.a[i][j]=;
     For(i,,) {
         rs.a[i][i]=;
         if(i<) bs.a[i+][i]=;
     }
     bs.a[][]=;
     bs.a[][]=;
     bs.a[][]=;
     bs.a[][]=;
     bs.a[][]=p-;
     ksm();
     LL ans=;
     For(i,,)
         (ans+=rs.a[i][]*f[i+]%p)%=p;
     printf("%lld\n",ans);
     Formylove;
 }