这个day1稍微有点毒瘤吧？？

DAY1

排列

以前总是把day1t1想太复杂了翻车，差不多往正解的方向想了一下感觉不可能这么复杂这可是noipday1t1啊一定有非常简单的方法然后翻车了？？

题目转换为求二分图完全匹配数，这个怎么都是十分不好算的，容易想到容斥。

用g[i]表示起码选了i条二分图的补图中的边的匹配数。

那么答案就是

$ans=\sum_{i=0}^{n}g[i]*(n-i)!*(-1)^i$

发现这个二分图的补图长得十分有特点啊。

这是若干条不想交的链构成的图，链与链之间互不影响可以单独考虑。用f[i][j][0/1]表示选到某条链上第i个点，在这条链上选了j条边，第i个点有没有被选的方案数，可以分别dp出每条链。

再跑个背包合并所有链，容斥统计答案即可。这样就有95分啦。

 //Achen

 #include<bits/stdc++.h>

 #define For(i,a,b) for(int i=(a);i<=(b);i++)

 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)

 #define Formylove return 0

 const int N=,p=;

 using namespace std;

 typedef long long LL;

 typedef double db;

 int n,K,g[N],f[N][N][],fac[N],up;

 template<typename T>void read(T &x) {

     T f=; x=; char ch=getchar();

     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();

     if(ch=='-') f=-,ch=getchar();

     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;

 }

 int mo(LL x) { return x%p; }

 #define ANS

 int main() {

 #ifdef ANS

     freopen("permutation.in","r",stdin);

     freopen("permutation.out","w",stdout);

 #endif

     read(n); read(K);

     fac[]=;

     For(i,,n) fac[i]=(LL)fac[i-]*i%p;

     g[]=;

     For(cs,,) {

         For(i,,K) {

             f[i][][]=;

             int j=i,tp=;

             for(;;) {

                 if(j+K>n) break;

                 j+=K; tp++;

                 For(l,,tp) {

                     f[j][l+][]=f[j-K][l][];

                     f[j][l][]=(f[j-K][l][]+f[j-K][l][])%p;

                 }

             }

             Rep(l,up,) For(m,,tp)

                 (g[l+m]+=(LL)g[l]*(f[j][m][]+f[j][m][])%p)%=p;

             up+=tp;

         }

     }

     LL ans=;

     For(i,,n) {

         //printf("%d\n",g[i]);

         if(!(i&)) ans=(ans+(LL)g[i]*fac[n-i]%p)%p;

         else ans=(ans-(LL)g[i]*fac[n-i]%p+p)%p;

     }

     printf("%lld\n",ans);

     //cerr<<clock()<<endl;

     Formylove;

 }

95pt

最后5分是什么玩意？？？

发现这些链其实都一样，只有两种不同长度。又发现刚才的dp其实是一个组合数。就可以直接得到最后要的f,然后多项式exp就可以的答案了。noip之前我不会写这种玩意的。。

苹果树

最简单的一道题，很多人都切了，而且各种方法花样AC。大致有标解的lca,辉神李巨的树剖，ljq的伪tarjan。但是我的随机化应该是最好写的，写得快点大概5min就够了吧。

 //Achen

 #include<bits/stdc++.h>

 #define For(i,a,b) for(int i=(a);i<=(b);i++)

 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)

 #define Formylove return 0

 const int N=;

 using namespace std;

 typedef unsigned long long LL;

 typedef double db;

 int n,m;

 template<typename T>void read(T &x) {

     T f=; x=; char ch=getchar();

     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<<];

 void add(int u,int v) {

     nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v;

     nxt[++ecnt]=fir[v]; fir[v]=ecnt; to[ecnt]=u;

 }

 map<LL,int>mp,mp2;

 LL val[N],val2[N],ans;

 void dfs(int x,int fa) {

     for(int i=fir[x];i;i=nxt[i]) if(to[i]!=fa) {

         dfs(to[i],x);

         val[x]=(val[x]^val[to[i]]);

     }

     if(fa) {

         if(!val[x])

             ans+=m;

         else if(mp[val[x]])

             ans++;

     }

 }

 #define ANS

 int main() {

 #ifdef ANS

     freopen("tree.in","r",stdin);

     freopen("tree.out","w",stdout);

 #endif

     srand();

     read(n); read(m);

     For(i,,n) {

         int u,v;

         read(u); read(v);

         add(u,v);

     }

     For(i,,m) {

         int u,v;

         read(u); read(v);

         LL t;

         t=(((LL)rand()<<)|rand());

         while(!t) {

             t=(((LL)rand()<<)|rand());

         }

         mp[t]=;

         val[u]^=t;

         val[v]^=t;

     }

     dfs(,);

     printf("%llu\n",ans);

     //cerr<<clock()<<endl;

     Formylove;

 }

多项式

披着斗地主外皮的多项式？？瑟瑟发抖。

模拟只给50，我大概写模拟都要写3h+把。。。写了个30分的暴力还写挂了，没把不同花色的同种对子区分开。。。

还是noip之前不改系列。

 //Achen

 #include<bits/stdc++.h>

 #define For(i,a,b) for(int i=(a);i<=(b);i++)

 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)

 #define Formylove return 0

 const int N=,p=;

 using namespace std;

 typedef long long LL;

 typedef double db;

 int T,n,cnt[N],cc[N],C[][],pr[];

 LL ans;

 char s[N];

 template<typename T>void read(T &x) {

     T f=; x=; char ch=getchar();

     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();

     if(ch=='-') f=-,ch=getchar();

     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;

 }

 void solve1() {

     ans=;

     For(i,,) cc[i]=;

     For(i,,) cc[cnt[i]]++;

     if(n==) ans=;

     else if(n==) {

         if(cc[]==) ans=;

         else ans=;

     }

     else if(n==) {

         if(cc[]==) ans=;

         else if(cc[]==) ans=;

         else ans=;

     }

     else if(n==) {

         ans++; //4 single

         if(cc[]==) ans+=; //4:1 3+1:C(4,3) 2+2:3 2+1+1:C(4,2)

         if(cc[]==) ans+=; //3+1 3+1 2+1+1:C(3,2)

         if(cc[]==) ans++; //2+1

         else if(cc[]==) ans+=; //2+2 2+1+1 2+1+1

     }

 }

 void dfs2(int pos,int num,int lim) {

     if(pos==) {

         ans=(ans+num)%p;

         return;

     }

     //if(cnt[pos]==2) dfs2(pos+1,num*2%p,0); // 1+1 or 2

     //else dfs2(pos+1,num,0); //

     if(cnt[pos]==) dfs2(pos+,num,);

     for(int j=pos+;cnt[j]&&j<=;j++) { //shun zi

         if(j>=pos+&&j>lim) {

             int N=num;

             For(k,pos,j) {

                 N=(LL)N*C[cnt[k]][]%p;

                 cnt[k]--;

             }

             dfs2(pos,N,j);

             For(k,pos,j) cnt[k]++;

         }

     }

     for(int j=pos;cnt[j]==&&j<=;j++) { //lian dui

         if(j>=pos+) {

             dfs2(j+,num,);

         }

     }

 }

 #define ANS

 int main() {

 #ifdef ANS

     freopen("polynomial.in","r",stdin);

     //freopen("polynomial.out","w",stdout);

 #endif

     read(T);

     For(i,,) C[i][]=;

     pr[]=;

     For(i,,) pr[i]=pr[i-]*%p;

     For(i,,) For(j,,i) C[i][j]=(C[i-][j-]+C[i-][j])%p;

     while(T--) {

         read(n);

         memset(cnt,,sizeof(cnt));

         For(i,,n) {

             scanf("%s",s);

             int len=strlen(s);

             if(len==) {

                 if(s[]=='J') cnt[]++;

                 else if(s[]=='Q') cnt[]++;

                 else if(s[]=='K') cnt[]++;

                 else if(s[]=='A') cnt[]++;

                 else {

                     if(s[]-''==) cnt[]++;

                     else cnt[s[]-'']++;

                 }

             }

             else {

                 if(s[]=='') cnt[]++;

                 else cnt[]++;

             }

         }

         For(i,,n) scanf("%s",s);

         ans=;

         if(n<=) solve1();

         else dfs2(,,);

         printf("%lld\n",ans);

     }

     //cerr<<clock()<<endl;

     Formylove;

 }

30pt

DAY2

day2其实并不难，应该是容易AK的，但是DCOI实在太菜了。。

我写t2的时候把线段树的时空复杂度算错了，就想写分块水水，码力太弱（本质是没想清楚），花了太多时间最终却没调出来，只交了暴力，还导致本来很简单的t3根本没时间思考，写了暴力就交。

序列

这种玩意套路的做法就是枚举每个起点i，然后分别统计后面a[j]-(j-i+1)大于0的和小于0的和，前面类似。

也就是每个位置的权值是a[i]-i，每次查询大于或小于某个权值的数的权值和。

树状数组肯定是可以做的，常数优秀一点说不定能过呢。

发现起点一位一位的移时，每次查询的那个数是一点点增大或减小的，而ai又很小，拿cnt数组记一下，每次把刚刚越过线的部分算一算就可以On做了。

 //Achen

 #include<bits/stdc++.h>

 #define For(i,a,b) for(int i=(a);i<=(b);i++)

 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)

 #define Formylove return 0

 const int N=4e6+;

 typedef long long LL;

 typedef double db;

 using namespace std;

 LL n,a[N],cnt[N],pos;

 LL ans[N],now,sum,sumcnt,nowcnt;

 template<typename T>void read(T &x) {

     char ch=getchar(); T f=; x=;

     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("a.in","r",stdin);

     freopen("a.out","w",stdout);

 #endif

     read(n);

     For(i,,n) read(a[i]);

     Rep(i,n,) {

         pos=-i;

         sum+=a[i]-i; sumcnt++;

         if(a[i]-i>=pos) now+=a[i]-i,nowcnt++;

         ans[i]=(now-nowcnt*pos)+((sumcnt-nowcnt)*pos-(sum-now));

         cnt[a[i]-i+n]++;

         now-=pos*cnt[pos+n];

         nowcnt-=cnt[pos+n];

     }

     sum=now=sumcnt=nowcnt=;

     memset(cnt,,sizeof(cnt));

     For(i,,n) {

         pos=n-i;

         sum+=a[i]-i; sumcnt++;

         if(a[i]-i<=pos) now+=a[i]-i,nowcnt++;

         ans[i+]+=(nowcnt*pos-now)+((sum-now)-(sumcnt-nowcnt)*pos);

         cnt[a[i]-i+n]++;

         now-=pos*cnt[pos+n];

         nowcnt-=cnt[pos+n];

     }

     LL rs=ans[];

     For(i,,n) rs=min(rs,ans[i]);

     printf("%lld\n",rs);

     Formylove;

 }

迷宫

线段树时空复杂度都怎么算怎么ok啊，而且超级好写的啊。。

就是有点卡常，我要加register才能卡过最后两个点。

 //Achen

 #include<bits/stdc++.h>

 #define For(i,a,b) for(register int i=(a);i<=(b);i++)

 #define Rep(i,a,b) for(register int i=(a);i>=(b);i--)

 #define Formylove return 0

 #define inf 1e8

 const int N=;

 typedef long long LL;

 typedef double db;

 using namespace std;

 int n,m,q,a[][N];

 template<typename T>void read(T &x) {

     char ch=getchar(); T f=; x=;

     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();

     if(ch=='-') f=-,ch=getchar();

     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;

 }

 void GM(int &x,int y) { if(x>y) x=y; }

 #define lc x<<1

 #define rc ((x<<1)|1)

 #define mid ((l+r)>>1)

 int sg[N<<][][];

 void upd(int x,int l,int r) {

     memset(sg[x],/,sizeof(sg[x]));

     if(l==r) {

         For(i,,n) if(a[i][l]) sg[x][i][i]=;

         For(i,,n) Rep(j,i-,) if(a[i][l]&&a[i-][l]) {

             GM(sg[x][i][j],sg[x][i-][j]+);

             GM(sg[x][j][i],sg[x][i-][j]+);

         }

     }

     else {

         For(i,,n) For(j,,n) {

             For(k,,n) if(a[k][mid]&&a[k][mid+])

                 GM(sg[x][i][j],sg[lc][i][k]+sg[rc][k][j]+);

         }

     }

 }

 void build(int x,int l,int r) {

     if(l==r) { upd(x,l,r); return; }

     build(lc,l,mid); build(rc,mid+,r);

     upd(x,l,r);

 }

 void update(int x,int l,int r,int p1,int p2) {

     if(l==r) {

         a[p2][p1]^=;

         upd(x,l,r); return;

     }

     if(p1<=mid) update(lc,l,mid,p1,p2);

     else update(rc,mid+,r,p1,p2);

     upd(x,l,r);

 }

 int tp[][],rs[][],ok;

 void merge(int x,int l,int r) {

     if(!ok) {

         ok=;

         For(i,,n) For(j,,n)

             tp[i][j]=sg[x][i][j];

     }

     else {

         memset(rs,/,sizeof(rs));

         For(i,,n) For(j,,n) {

             For(k,,n) if(a[k][l-]&&a[k][l])

                 GM(rs[i][j],tp[i][k]+sg[x][k][j]+);

         }

         memcpy(tp,rs,sizeof(rs));

     }

 }

 void qry(int x,int l,int r,int ql,int qr) {

     if(l>=ql&&r<=qr) { merge(x,l,r); return; }

     if(ql<=mid) qry(lc,l,mid,ql,qr);

     if(qr>mid) qry(rc,mid+,r,ql,qr);

 }

 #define ANS

 int main() {

 #ifdef ANS

     freopen("maze.in","r",stdin);

     freopen("maze.out","w",stdout);

 #endif

     read(n); read(m); read(q);

     For(i,,n) For(j,,m) read(a[i][j]);

     build(,,m);

     For(cs,,q) {

         int o,x,y,z,w;

         read(o); read(x); read(y);

         if(o==)

             update(,,m,y,x);

         else {

             read(z); read(w);

             ok=;

             qry(,,m,y,w);

             if(tp[x][z]>=inf) puts("-1");

             else printf("%d\n",tp[x][z]);

         }

     }

     Formylove;

 }

动态数点

从n^2暴力优化到标解蛮容易的，二分长度，st表处理gcd和min值。只是似乎又要卡常了。

看李巨一个排序的复杂度爆踩std。

按a从小到大排序，然后拓展。拓展过的点就不用再拓展了，发现每个点最多被拓展两边，左边一遍右边一遍（3 15 5）。

辉神就更强了，从左到右依次拓展但是每次拓展的中心是上一次拓展的右边界的下一个，这样最坏的情况下也只能卡成nlogv的，而实际远不到上限，跑得飞快。

 //Achen

 #include<bits/stdc++.h>

 #define For(i,a,b) for(int i=(a);i<=(b);i++)

 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)

 #define Formylove return 0

 const int N=;

 typedef long long LL;

 typedef double db;

 using namespace std;

 int n,b[N],vis[N],a[N];

 int tot,ans;

 template<typename T>void read(T &x) {

     char ch=getchar(); T f=; x=;

     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();

     if(ch=='-') f=-,ch=getchar();

     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;

 }

 struct node {

     int pos,a;

     friend bool operator <(const node&A,const node&B) {

         return A.a<B.a;

     }

 }p[N];

 #define ANS

 int main() {

 #ifdef ANS

     freopen("point.in","r",stdin);

     freopen("point.out","w",stdout);

 #endif

     read(n);

     For(i,,n) {

         read(a[i]);

         p[i].a=a[i];

         p[i].pos=i;

     }

     sort(p+,p+n+);

     For(i,,n) {

         int l=p[i].pos,r=p[i].pos;

         if(vis[l]) continue;

         while(l>&&a[l-]>=p[i].a&&a[l-]%p[i].a==)

             vis[--l]=;

         while(r<n&&a[r+]>=p[i].a&&a[r+]%p[i].a==)

             vis[++r]=;

         b[l]=r-l;

         ans=max(ans,r-l);

     }

     For(i,,n) tot+=(b[i]==ans);

     printf("%d %d\n",tot,ans);

     For(i,,n) if(b[i]==ans) printf("%d ",i);

     puts("");

     Formylove;

 }

李巨做法的代码

NOIp2018集训test-10-4/test-10-5 (联考四day1/day2)的更多相关文章

  1. NOIp2018集训test-10-6/test-10-7 (联考五day1/day2)

    昨天考完月考,明天初赛,dcoi2017级今天终于开始停课准备noip了,大概没有比本弱校停课更晚的学校了吧.本来就够菜了,怕是要凉透哦. DAY1 T1石头剪刀布 据说爆搜随便做,但是我觉得我的O( ...

  2. NOIp2018集训test-9-22(am/pm) (联考三day1/day2)

    szzq学长出的题,先orz一下. day1 倾斜的线 做过差不多的题,写在我自己的博客里,我却忘得一干二净,反而李巨记得清清楚楚我写了的. 题目就是要最小化这个东西 $|\frac{y_i-y_j} ...

  3. NOIp2018集训test-10-21 (联考六day1)

    今天被高一狂踩,两个手抖,t1一个1写成2,t3一个+=写成=,所谓失之毫厘谬以千里,直接丢了50分. 完全背包 看到背包体积如此之大物品体积如此之小容易很想到贪心,肯定要先加很多很多的性价比最高的最 ...

  4. NOIp2018集训test-9-15(联考二day1)

    T1.矩阵游戏 水题.每一行最后乘的数为x[i],每一列为y[i],暴力算第一行的列的贡献,每一行的列的贡献是公差为所有列的贡献之和的等差数列,然后每一行再乘上行的贡献求和即为答案. //Achen ...

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

    又被辉神吊打了.今天不仅被辉神李巨吊打,还给基本上给全班垫底了. 看到T3就知道是十进制快速幂,全机房考试的当时应该就我会,结果我tm没找到递推. Orz lyc BM直接水过,Orz wys六个fo ...

  6. 给定整数a1、a2、a3、...、an，判断是否可以从中选出若干个数，使得它们的和等于k（k任意给定，且满足-10^8 <= k <= 10^8）。

    给定整数a1.a2.a3.....an,判断是否可以从中选出若干个数,使得它们的和等于k(k任意给定,且满足-10^8 <= k <= 10^8). 分析:此题相对于本节"寻找满 ...

  7. 2018.10.26 浪在ACM 集训队第四次测试赛

    2018.10.26 浪在ACM 集训队第四次测试赛 题目一览表 来源 考察知识点 完成时间 A 生活大爆炸版 石头剪刀布  NOIP 提高组 2014   模拟???  2018.11.9 B 联合 ...

  8. gitlab 迁移、升级打怪之路：8.8.5--> 8.10.8 --> 8.17.8 --> 9.5.9 --> 10.1.4 --> 10.2.5

    gitlab 迁移.升级打怪之路:8.8.5--> 8.10.8 --> 8.17.8 --> 9.5.9 --> 10.1.4 --> 10.2.5 gitlab 数据 ...

  9. var a =10 与 a = 10的区别

    学习文章------汤姆大叔-变量对象 总结笔记 变量特点: ①变量声明可以存储在变量对象中.②变量不能直接用delete删除. var a =10 与 a = 10的区别: ①a = 10只是为全局 ...

随机推荐

  1. fiddler增加ip以及响应时间列

    最近打算看一下移动端app的响应等请求,这里打算用fillder来查看appium的模拟出发请求的操作来查看结果, 所以我们需要在左侧的面板增加我们所需要的ip,响应时间等数据以方便我们查看 fidd ...

  2. Python 十大装 B 语法解析

    Python 是一种代表简单思想的语言,其语法相对简单,很容易上手.不过,如果就此小视 Python 语法的精妙和深邃,那就大错特错了.本文精心筛选了最能展现 Python 语法之精妙的十个知识点,并 ...

  3. Oracle系列教程之一软件安装与卸载

    1.安装软件 软件下载地址:https://www.oracle.com/technetwork/cn/database/enterprise-edition/downloads/index.html ...

  4. psecurity配置

    <?xml version="1.0" encoding="UTF-8"?><beans xmlns="http://www.spr ...

  5. 【LeetCode 10】正则表达式匹配

    题目链接 [题解] 看到这个题解 写的代码. 就是加个备忘录法.优化一下暴力的做法. 匹配的时候如果遇到号的话,就两种可能.不再考虑它前面一个字符了. 跳过这个或者.或者继续用前面那个字符匹配. 即d ...

  6. 安装nodejs nvm

    安装nodejs sudo apt-get install nodejs sudo apt-get install npm 安装nvm https://www.runoob.com/w3cnote/n ...

  7. win7 cmd 常用命令

    进入不同的分区  d: 查看之栏目树:tree 查看当前目录下的子目录:dir 切换不同的目录:cd

  8. windows 驱动开发 DDK与WDK WDM的区别

    1.首先,先从基础的东西说起,开发WINDOWS下的驱动程序,需要一个专门的开发包,如:开发JAVA程序,我们可能需要一个JDK,开发WINDOWS应用程序,我们需要WINDOWS的SDK,现在开发W ...

  9. JAVA学习之进制转换练习

    public static void main(String[] args) { toBin(60); toBa(60); toHex(60); } /** 十进制-->二进制 */ publi ...

  10. jQuery 封装的ajax

    jquery封装的ajax 具体操作: $.get(url  [,data]  [,fn回调函数]   [, dataType]); data:给服务器传递的数据,请求字符串 .json对象 都可以设 ...