文章标题

最新推荐文章于 2022-08-04 15:18:28 发布
最新推荐文章于 2022-08-04 15:18:28 发布
阅读量447 收藏 1
点赞数
分类专栏： dp 文章标签： dp 快速幂
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sue2333/article/details/53190536
版权

T1

LGTB 与序列
LGTB 得到了一个序列，他想在这个序列中选择一个最长的连续子序列，使得这个子序列的最大公约数等于
1。请告诉他他能得到的最大长度，如果没有这样的序列，输出��1
输入
输入第一行包含一个整数n 代表序列大小
接下来一行，包含n 个整数a1, a2, …, an，代表序列
对于50% 的数据，1 n 1000
对于100% 的数据，1 n 105 1 ai 109
输出
输出包含一个整数l，代表最长的连续子序列，如果无解请输出��1
样例
样例输入
2
7 2
样例输出
2
样例输入
3
2 2 4
样例输出
-1

一道水题…不过我不太懂n=1的时候答案居然是1？excuse me？
代码：

#include<iostream>
#include<cstdio>
#define LL long long
using namespace std;
const int N=1e5+5;
int n,s[N],ans;
int gcd(int a,int b)
{
    if(!b)return a;
    else return gcd(b,a%b);
}
void worrk()
{
    int hs=s[1];
    for(int i=2;i<=n;i++)
    {
        hs=gcd(hs,s[i]);
        if(hs==1){ans=n;return;}
    }
}
int main()
{
    freopen("seq.in","r",stdin);
    freopen("seq.out","w",stdout);
    scanf("%d",&n);ans=-1;
    for(int i=1;i<=n;i++)
        scanf("%d",&s[i]);
    worrk();
    printf("%d",ans);
    return 0;
}

T2

LGTB 与桌子
LGTB 新买了一张n m 的矩(桌) 阵(子)，他想给某些1 1 的小矩形染色，使得染色之后，原矩阵的每
个n n 的子矩阵中都包含恰好k 个被染色了的小矩形
他想知道有多少种染色方案能让他满足上述要求
因为答案可能很大，请输出方案数模1000000007 (109 + 7) 后的值
输入
输入第一行包含三个整数n, m, k 意义如题面所示
对于15% 的数据，1 n m 20, n m
对于40% 的数据，1 n 10, n m 1000
对于100% 的数据，1 n 100, n m 1018, 0 k n2
输出
输出包含1 个整数，代表LGTB 的方案数
样例
样例输入
5 6 1
样例输出
45

这道题…考试时想出了正解….好开心….然后…一个小地方忘了mod…结果….爆了…TAT我要哭死哭死哭死,…

下面来讲讲做法：
1.因为使n*n的涂色数一样，所以第1列与n+1列等效，2与n+2等效，同理一直等效下去。等效意义为涂色数相同。

2.由此易得方程f[i][j]=∑(q=1~min(n,j))f[i][j-q]*C(n,q)^h[i]。
下面来解释一下这个方程：
f[i][j]表示前i列共涂了j个色块，而q则表示当前列要涂q个，所以有C(n,q)种涂法，而h[i]表示与当前列等效的个数，即m/n或m/n+1(i<=m%n)。表示的意思为：要涂q个，而与它等效的h[i]个同样每种都有C(n,q)种，故相乘。

3.C(n,q)^h[i]可以预处理，因为h[i]只有两种。
C(n,q)=(n-q+1)/k*C(n,q-1)，这里在mod的时候要用到逆元，一个循环处理一下就ok，接下来就快速幂了。当然其他推方法也可以，我个人认为这种推的不容易错。

4.复杂度大概O(n^4)，实际没有这么大，中间很多地方for不到√所以可以放心写，最慢的一组大概一秒多，还是蛮快啦。

代码：

#include<iostream>
#include<cstdio>
#define LL long long
using namespace std;
const int N=105,M=1e4+5,mod=1e9+7;
int n,k,w,nie[N];
LL f[2][M],c[M],r[2][M],m;
int ys;
void ksm(int i,int cc)
{
    LL d=c[i],ans=1;
    while(cc)
    {
        if(cc&1)ans=(ans*d)%mod;
        cc>>=1;d=(d*d)%mod;
    }
    r[0][i]=ans%mod; //C(n,i)^(m/n)
    r[1][i]=((LL)ans*(LL)c[i])%mod;//C(n,i)^(m/n+1)
}
void pre()
{
    ys=(LL)m%n;int c1=(LL)((LL)m/n)%(mod-1);
    c[0]=1;r[0][0]=r[1][0]=1;
    nie[1]=1;
    for(int i=2;i<=n;i++)
    nie[i]=(LL)((LL)(mod-mod/i)*(LL)nie[mod%i])%mod;//逆元
    for(int i=1;i<=n;i++)
    {
        c[i]=(c[i-1]*(LL)(n-i+1)%mod*(LL)nie[i])%mod;
        ksm(i,c1);
    }
}
void worrk()
{
    int pre=0,now=1;f[0][0]=f[1][0]=1;
    for(int i=1;i<=n;i++)
    {
        int yy=min(k,i*n),t=0;
        if(i<=ys)t=1;
        for(int j=0;j<=yy;j++)
        {
            f[now][j]=0;
            int y=min(j,n);
            for(int q=0;q<=y;q++)
            f[now][j]=(f[now][j]+((LL)f[pre][j-q]*(LL)r[t][q])%mod)%mod;
        }
        swap(now,pre);
    }
    cout<<f[pre][k]<<endl;
}
int main()
{
    freopen("table.in","r",stdin);
    freopen("table.out","w",stdout);
    scanf("%d",&n);cin>>m;scanf("%d",&k);
    pre();
    worrk();
}

直接开数组是不会超的，习惯开滚动了。

T3

LGTB 与正方形
LGTB 最近迷上了正方形，现在他有n 个在二维平面上的点，请你告诉他在这些点中选4 个点能组成四条边
都平行于坐标轴的正方形的方案数有多少
输入
输入第一行包含一个整数n，代表点的数量
接下来n 行每行包含两个整数xi, yi，代表点的坐标
对于10% 的数据，1 n 50
对于30% 的数据，1 n 1000
对于100% 的数据，1 n 105，0 xi, yi 105
数据保证没有两点在同一位置
输出
输出包含一个整数代表方案数
样例
样例输入
5
0 0
0 2
2 0
2 2
1 1
样例输出
1
样例输入
9
0 0
1 1
2 2
0 1
1 0
0 2
2 0
1 2
2 1
样例输出
5

这道题先贴暴力代码吧，30分好像是。
排序查找。

#include<iostream>
#include<cstdio>
#include<algorithm>
#define LL long long
using namespace std;
const int N=1e5+5;
int n,xi,yi,tx[N],ty[N];
LL ans;
struct node{
    int x,y,num;
};
node xx[N],yy[N];
bool comp1(const node&a,const node&b)
{
    if(a.x==b.x)return a.y<b.y;
    else return a.x<b.x;
}
bool comp2(const node&a,const node&b)
{
    if(a.y==b.y)return a.x<b.x;
    else return a.y<b.y;
}
LL pd(int x1,int x2)
{
    int cj=xx[x2].y-xx[x1].y;
    int xw1=xx[x1].x,xw2=xx[x2].x;
    int yw1=xx[x1].y,yw2=xx[x2].y;
    int y1=ty[xx[x1].num]+1,y2=ty[xx[x2].num]+1;
    LL t1=0,t2=0;
    while(yy[y1].y==yw1)
    {
        if(yy[y1].x-xw1==cj)t1++;
        y1++;
    }
    while(yy[y2].y==yw2)
    {
        if(yy[y2].x-xw2==cj)t2++;
        y2++;
    }
    return t1*t2;
}

void worrk()
{
    for(int i=1;i<n;i++)
    {
        int w=i+1;
        while(xx[w].x==xx[i].x)
        {
            ans+=pd(i,w);
            w++;
        }
    }
    cout<<ans<<endl;
}
int main()
{
    freopen("square.in","r",stdin);
    freopen("square.out","w",stdout);
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&xi,&yi);
        xx[i].x=yy[i].x=xi;
        xx[i].y=yy[i].y=yi;
        xx[i].num=yy[i].num=i;
    }
    sort(xx+1,xx+n+1,comp1);
    sort(yy+1,yy+n+1,comp2);
    for(int i=1;i<=n;i++)
    {
        tx[xx[i].num]=i;
        ty[yy[i].num]=i;
    }
    worrk();
}

正解懒得写了…
贴上来自何神的解题报告：

迷之方法（很有道理）
将原坐标通过x分类，设每个分类为P[x]
如果P[x]中数的个数小于sqrt(n)，直接暴力枚举里面元素（x相同，判断y获取正方形）
另外的P[x]暴力枚举这些集合，判断其中那些是y相等的
对于第一类小于sqrt(n)的点的集合，每个集合里面最多sqrt(n)个元素，所以该时间复杂度为O(n sqrt(n))
对于第二类，由于每一种size(P[x])>=sqrt(n)，这种P[x]有最多sqrt(n)种，枚举这sqrt(n)种P[x]，O(sqrt(n)^2) = O(n)，枚举某一个集合中的所有元素，判断能否构成正方形就可以了该类O(n sqrt(n))
总时间复杂度O(n sqrt(n)) GGGG

迷之标程：

#include <set>
#include <cmath>
#include <stack>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cstdlib>
#include <numeric>
#include <vector>
#include <ctime>
#include <queue>
#include <list>
#include <map>
#define pi acos(-1.0)
#define INF 0x3f3f3f3f
#define clr(x)  memset(x,0,sizeof(x));
#define clrto(x,siz,y)  for(int xx=0;xx<=siz;xx++)  x[xx]=y;
#define clrset(x,siz)  for(int xx=0;xx<=siz;xx++)  x[xx]=xx;
#define clr_1(x) memset(x,-1,sizeof(x));
#define clrmax(x) memset(x,0x3f,sizeof(x));
#define clrvec(x,siz) for(int xx=0;xx<=siz;xx++)  x[xx].clear();
#define fop2   freopen(".in","r",stdin); //freopen(".out","w",stdout);
#define fop   freopen("in.txt","r",stdin);//freopen("out.txt","w",stdout);
#define myprogram By_135678942570
#define clrcpy(x,siz,y)  for(int xx=0;xx<siz;xx++)  x[xx]=y[xx];
#define pb push_back
using namespace std;
struct point
{
    int x,y;
}P[100011];
int SIZ=400;
struct hash_map
{
    const static int mod=100007;
    int head[mod];
    struct hash_tables
    {
        long long key;
        int val;
        int next;
    }ele[100007];
    int N;
    int getHash(long long x)
    {
        return x%mod;
    }
    void init()
    {
        memset(head,255,sizeof(head));
        N=0;
    }
    int fint(long long x)
    {
        for(int i=head[getHash(x)];i!=-1;i=ele[i].next)
            if(ele[i].key==x)
                return i;
        return -1;
    }
    void insert(long long x)
    {
        int tmp=getHash(x);
        ele[N].key=x;
        ele[N].val=0;
        ele[N].next=head[tmp];
        head[tmp]=N++;
    }
    int& operator [](long long x)
    {
        int tmp=fint(x);
        if(tmp==-1)
        {
            insert(x);
            return ele[N-1].val;
        }
        else return ele[tmp].val;
    }
}HT;
vector<int>row[100111];
vector<int>GREAT;
vector<int>LESS;
int main()
{
    // fop;
    freopen("square.in", "r", stdin);
    freopen("square.out", "w", stdout);
    HT.init();
    int n;
    scanf("%d",&n);
    for(int i=0;i<n;i++)
    {
        scanf("%d%d",&P[i].x,&P[i].y);
        row[P[i].x].pb(P[i].y);
        HT.insert(P[i].x*1000003ll+P[i].y);
    }
    for(int i=0;i<=100000;i++)
        if(row[i].size()>=SIZ)
            GREAT.pb(i);
        else if(row[i].size()>0)
            LESS.pb(i);
    int res=0;
    for(int i=0;i<LESS.size();i++)
    {
        int rowid=LESS[i];
        for(int j=0;j<row[rowid].size();j++)
            for(int k=j+1;k<row[rowid].size();k++)
            {
                int x1=rowid;
                int y1=row[rowid][j],y2=row[rowid][k];
                int x2=x1+abs(y1-y2);
                int x22=x1-abs(y1-y2);
                int f1=HT.fint(x2*1000003ll+y1);
                int f2=HT.fint(x2*1000003ll+y2);
                if(f1!=-1&&f2!=-1)
                {
                    res++;
                }
                if(x22>=0&&row[x22].size()>=SIZ)
                {
                    f1=HT.fint(x22*1000003ll+y1);
                    f2=HT.fint(x22*1000003ll+y2);
                    if(f1!=-1&&f2!=-1)
                        res++;
                }
            }
    }
    for(int i=0;i<GREAT.size();i++)
        sort(row[GREAT[i]].begin(),row[GREAT[i]].end());
    for(int i=0;i<GREAT.size();i++)
        for(int j=i+1;j<GREAT.size();j++)
        {
            int x1=GREAT[i];
            int x2=GREAT[j];
            int len=abs(x1-x2);
            int ii=0,jj=0;
            while(ii<row[x1].size()&&jj<row[x2].size())
            {
                if(row[x1][ii]==row[x2][jj])
                {
                    int f1=HT.fint(x1*1000003ll+row[x1][ii]+len);
                    int f2=HT.fint(x2*1000003ll+row[x2][jj]+len);
                    if(f1!=-1&&f2!=-1)
                        res++;
                    ii++,jj++;
                }
                else if(row[x1][ii]>row[x2][jj])
                    jj++;
                else ii++;
            }
        }
    printf("%d\n",res);
}

题外话：离noip还有两天！fighting！！
Do what I love.
Fight for what I want.

梳子很怕冷x
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
使用POI生成文章标题
Cloud的专栏
09-24 7122
生成文章标题有两种方式，一种是自定义标题，还有一种是使用现有样式来生成，现在网上大部分是第一种方式，但这种方式有很多弊端，我们先来看第一种 第一步，引入依赖 <dependency> <groupId>org.apache.poi</groupId> <artifactId>poi</artifactId> <version>3.9</version> </dependency>
Markdown 跳转到本文章标题
小白学CS
10-30 4719
Markdown 跳转到本文章标题
python爬虫爬取微信公众号的阅读数、喜爱数、文章标题和链接等信息
weixin_46569877的博客
08-04 9554
python爬取微信公众号
火车头插件：根据文章标题生成图片；标题生成相关配图
qq_26057339的博客
05-12 1742
本插件可以根据文章标题自动配图；根据文章标题生成一张好看的插图； 插件售价200；不限制ip,不限域名。一次购买终身使用。 没有第三方接口的调用不用担心插件失效。 插件效果： 插件优势： 用户可以自定义文字的颜色 用户可以自定义文字的大小 用户可以自定义图片背景 用户可以自定义图片格式 用户可以自定义图片储存的目录 用户可以自定义图片在文章哪个位置插入 插件分为 火车头插件和服务器接口： 服务器接口负责处理图片和储存图片； 火车头插件处理文章内容 服务器插件有运行环境要求 php 版本低于
RiPro主题文章标题美化教程 WordPress美化教程
03-02
/*文章标题优化*/ .sjblog-avatar {  width: 50px;  float: left; } .sjblog-avatar img {  float: left;  width: 50px;  border-radius: 50%;  margin-right: 20px; } .sjblog-tgroup{ ...
Wordpress自定义函数控制文章标题长度
09-29
主要为大家介绍了Wordpress自定义函数控制文章标题长度,需要的朋友可以参考下
Python3 文章标题关键字提取的例子
01-20
1.读取所有文章标题； 2.用“结巴分词”的工具包进行文章标题的词语分割； 3.用“sklearn”的工具包计算Tf-idf（词频-逆文档率）; 4.得到满足关键词权重阈值的词 结巴分词详见：结巴分词Github sklearn详见：文本...
PHPCMS内容页文章标题添加页码的方法教程
09-29
如何让PHPCMS的内容页标题显示页码呢？让每篇文章的标题都不一样，这样不但提高了网页体验率，同时也可以达到SEO优化的效果，让百度、谷歌等多收录一些。
ChatGPT4.0知识问答、DALL-E生成AI图片、Code Copilot辅助编程，打开新世界的大门.txt
04-30
ChatGPT4.0知识问答、DALL-E生成AI图片、Code Copilot辅助编程，打开新世界的大门
基于matlab实现DOA 估计和自适应波束形成.rar
最新发布
04-30
基于matlab实现DOA 估计和自适应波束形成.rar
基于C++的线程安全容器。.zip
04-30
基于C++的线程安全容器。.zip
华为数字化转型实践28个精华问答glkm.pptx
04-30
华为数字化转型实践28个精华问答glkm.pptx
本周-综合案例.zip
04-30
本周-综合案例.zip
基于Swift简单易上手的iOS开发框架.zip
04-30
基于Swift简单易上手的iOS开发框架.zip
liba52-0-32bit-0.7.5+svn613-1.19.x86_64.rpm
04-30
liba52-0-32bit-0.7.5+svn613-1.19.x86_64
271ssm-mysql-jsp 动漫之家系统.zip（可运行源码+数据库文件+文档）
04-30
本次的设计主要是通过对动漫系统开发的背景、现状进行了分析，总结出了本次动漫之家系统开发的意义。根据此次开发的目的和意义，本次的系统开发选择了SSM框架、HTML5以及idea平台来进行动漫之家系统的开发，通过MySQL来进行数据库的开发。通过对整个动漫之家系统进行功能需求的调查研究，通过对此次的系统开发进行可行性的分析。通过实体图模型以及功能结构模型来对本次的系统开发进行了整体的开发。在整个系统开发完毕之后，通过截图说明的方式来进行系统功能的介绍，最后通过系统测试来对本次系统的完整性进行测试，最终通过本次的开发，整个动漫之家系统可以实现很好的运行，起到了为动漫爱好者提供动漫资讯的功能运行。 在前端的系统开发上，主要是为了给动漫爱好者们提供一个在线交流、在线观看动漫、在线购买动漫周边的综合性服务平台。通过这个平台，可以通过注册成为会员后，在动漫内容下进行留言互动来实现更好的动漫交流与观后感的分享，可以提高站内用户对于该网站的使用兴趣。而后台则主要为该动漫之家系统的管理员提供管理服务，后台的管理中，管理员能 关键词：动漫之家；论坛网站；SSM框架；MySQL数据库
基于matlab数字图像处理的黄豆数量识别(GUI界面).zip代码48
04-30
基于matlab数字图像处理的黄豆数量识别(GUI界面)，基于matlab数字图像处理的黄豆数量识别(GUI界面)，基于matlab数字图像处理的黄豆数量识别(GUI界面)
wordpress 改写文章标题
03-20
您可以使用以下代码来改写 WordPress 文章标题： function custom_title( $title ) { $new_title = '新标题'; return $new_title; } add_filter( 'the_title', 'custom_title' ); 将“新标题”替换为您想要的...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值