2016中国大学生程序设计竞赛 - 网络选拔赛 题解(待续)

最新推荐文章于 2022-02-11 11:23:33 发布
最新推荐文章于 2022-02-11 11:23:33 发布
阅读量564 收藏
点赞数
分类专栏： 比赛题解 高斯消元
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/nike0good/article/details/52207308
版权

HDU 5832 A water problem

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<sstream>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
char n[10000000+10];

int main()
{
//  freopen("A.in","r",stdin);
//  freopen(".out","w",stdout);
    int kcase=1;
    while(~scanf("%s",n)) {
        printf("Case #%d: ",kcase++);
        int t=0;
        int len=strlen(n);
        Rep(i,len) {
            t=t*10+n[i]-'0';
            t%=10001;
        }
        if (!t) puts("YES") ; else puts("NO");
    }


    return 0;
}

HDU 5833 Zhu and 772002

There are n numbers a1,a2,…,an. The value of the prime factors of each number does not exceed 2000, you can choose at least one number and multiply them, then you can get a number b.

How many different ways of choices can make b is a perfect square number. The answer maybe too large, so you should output the answer modulo by 1000000007.
(1≤n≤300),(1≤ai≤10^18).

想到xor方程就行

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<sstream>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase) printf("Case #%d:\n",kcase);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 

#define MAXN (2000+10)
int p[MAXN],tot;
bool b[MAXN]={0};
void make_prime(int n)
{
    tot=0;
    Fork(i,2,n)
    {
        if (!b[i]) p[++tot]=i;
        For(j,tot)
        {
            if (i*p[j]>n) break;
            b[i*p[j]]=1;
            if (i%p[j]==0) break;  
        }
    }
}
ll A[400];
int n;
ull a[400][10];
bool cmp(int i,int j) {
    Rep(k,10) if (a[i][k]^a[j][k]) return a[i][k]<a[j][k];
    return 0;
}
void swap2(int i,int j) {
    Rep(k,10) if (a[i][k]^a[j][k]) swap(a[i][k],a[j][k]);
}
bool is_0(int i) {
    Rep(k,10) if (a[i][k]) return 0;
    return 1;
}
void upd_min(int K,int i) {  // a[k]=min(a[k],a[k]^a[i]);   
    int p=305;
    Rep(k,10) a[p][k]=a[K][k]^a[i][k];
    if (cmp(p,K)) swap2(p,K);
}
int gauss(int n) {
    For(i,n) {
        Fork(j,i+1,n) if (cmp(i,j)) swap2(i,j);
        if (is_0(i)) return i-1;
        For(k,n) if (i^k) {
            upd_min(k,i);
//          a[k]=min(a[k],a[k]^a[i]);           
        }
    }               
    return n;
}
ll pow2(ll a,ll b) {
    if (b==1) return a%F;
    if (!b) return 1%F;
    ll p=pow2(a,b/2);
    p=mul(p,p);
    if (b&1) p=mul(p,a);
    return p;
}
int main()
{
//  freopen("B.in","r",stdin);
//  freopen(".out","w",stdout);
    int T=read();
    make_prime(2000);
    For(kcase,T) {
        n=read();
        For(i,n) scanf("%I64d",&A[i]);
        MEM(a)
        For(i,n) {
            ull t=1;
            For(j,tot) {
                int p1=(j-1)/50,p2=(j-1)%50;
                if (!p2) t=1;
                bool b=0;
                while(A[i]%p[j]==0) {
                    A[i]/=p[j];
                    b=!b;
                }
                if (b) a[i][p1]|=t;
                t*=2;
            }
        }
        int p0=n-gauss(n);
        ll ans=pow2(2,p0);
        ans=sub(ans,1);

        Pr(kcase);
//      printf("%I64d\n",ans);
        cout<<ans<<endl;
    }


    return 0;
}

HDU 5838 Mountain

同BZOJ 2669([cqoi2012]局部极小值-状态压缩+dp)
http://blog.csdn.net/nike0good/article/details/52207269

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<sstream>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (772002)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase) printf("Case #%d: ",kcase);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define MAXN (100)
int n,m;
int p2[30];
ll f[MAXN][1<<9];
int g[MAXN],g2[MAXN],dir[8][2]={{1,0},{-1,0},{0,1},{0,-1},{1,1},{1,-1},{-1,1},{-1,-1}};
bool inside(int x,int y) {
    return 1<=x&&1<=y&&x<=n&&y<=m;
}
int id[MAXN][MAXN];
int idx(int i,int j) {
    return (i-1)*m+j-1;
}
int x[MAXN],y[MAXN],siz=0;
int b[30];
int calc() {
    memset(b,-1,sizeof(b));
    For(i,n) For(j,m) id[i][j]=idx(i,j);
    MEM(g) MEM(g2)
    Rep(i,siz) {
        int X=x[i],Y=y[i];
        g2[i]=1<<idx(X,Y);
        b[id[X][Y]]=i;
        Rep(di,8) {
            int a=X+dir[di][0],b=Y+dir[di][1];
            if (inside(a,b)) {
                int p=1<<idx(a,b);
                g[i]|=p;
            }
        }
    }

    MEM(f)
    f[0][0]=1;

    int S=1<<siz;
    int all=(1<<(n*m))-1;
    Rep(i,n*m) {
        Rep(st,S) if (f[i][st]) {
            int now=all;
            Rep(j,siz)
                if (st&p2[j]) {
                    now-=now&g2[j];
                } else {
                    now-=now&g[j];
                }
            int tot=0;
            int tot1=0,tot2=0;
            Rep(p,n*m) if (now&p2[p]) {
                tot++;
                if (-1 == b[p]) tot2++;
                else tot1++;
            }
            int lef=n*m-i;
            int choseX=siz-tot1;
            int useBlank=i-choseX;
            tot2-=useBlank;
            tot2=min(tot2,lef-tot1);
            if (tot2>0) 
                upd(f[i+1][st],mul(tot2,f[i][st]));


            Rep(p,n*m) {
                if (now&p2[p]) {
                    if (b[p]==-1) {
                    } else {
                        upd(f[i+1][st|p2[b[p]]],f[i][st]);
                    }
                }
            }

        }
    }
//  Rep(p,S) cout<<f[n*m][p]<<endl;
    return f[n*m][S-1];
} 
bool ma[6][6];
ll ans=0;
void dfs(int dep,int l) {
    if (dep==1) upd(ans,calc());
    else ans=sub(ans,calc());

    For(i,n) For(j,m) if (idx(i,j)>=l && !ma[i][j]){
        bool fl=0;
        Rep(di,8) {
            int a=i+dir[di][0],b=j+dir[di][1];
            if (!inside(a,b)) continue;
            if (ma[a][b]) { fl=1; break;
            }
        }
        if (!fl) {
            ma[i][j]=1;
            x[siz]=i,y[siz]=j;
            ++siz;
            dfs(-dep,idx(i,j)+1);
            ma[i][j]=0;
            --siz;
        }
    }
}
int main()
{
//  freopen("F.in","r",stdin);
//  freopen("F.out","w",stdout);
    p2[0]=1;
    For(i,25) p2[i]=p2[i-1]*2;
    int kcase=1;
    while(cin>>n>>m) {
        siz=0;
        For(i,n) {
            char s[10];
            scanf("%s",s);
            Rep(j,m) if (s[j]=='X') {
                x[siz]=i; y[siz]=j+1;
                ++siz; ma[i][j+1]=1;
            } else ma[i][j+1]=0;
        }

        Pr(kcase++);
        ans=0;
        if (siz<=9) dfs(1,0);       
        cout<<ans<<endl;
    }
    return 0;
}
nike0good
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
国际大学生程序设计竞赛中山大学内部选拔真题解（一）光盘
06-12
本压缩包文件“国际大学生程序设计竞赛中山大学内部选拔真题解（一）光盘”包含了中山大学在历次ICPC选拔赛中所使用的部分真题及解答，对于参赛者或者对算法和编程有兴趣的学生来说，这是一份极其宝贵的资源。...
国际大学生程序设计竞赛例题解 7 中山大学ICPC集训队内部选拔赛试题 2005-2006年
09-16
国际大学生程序设计竞赛例题解 7 中山大学ICPC集训队内部选拔赛试题 2005-2006年 郭嵩山 电子工业出版社
2016 ICPC 大连网络赛 部分题解
ddmxynl793173的博客
09-11 139
　　先讲1007，有m个人，n种石头，将n种石头分给m个人，每两个人之间要么是朋友关系，要么是敌人关系，朋友的话他们必须有一种相同颜色的石头，敌人的话他们必须所有石头的颜色都不相同。另外，一个人可以不拥有任何一种石头。求m个人的所有关系是不是都能用n种石头表示出来。比赛当时找的关系是n种石头可以表示n+1个人的关系。但是一直WA，因为考虑不周。 　　我们考虑这样的一种情况，我们把人分为...
2016 年青岛网络赛题解
LzyRapX
09-17 2525
2016 年青岛onlines题解1001 题意：给出一个整数nn, 找出一个大于等于nn的最小整数mm, 使得mm可以表示为2a3b5c7d2^a3^b5^c7^d.题解：快速幂预处理出所有形为2a3b5c7d2^a3^b5^c7^d即可, 大概只有几k个. 1002 题意：给 nn，求 ∑nk=11k2\sum_{k=1}^{n}\frac{1}{k^2}。题解：nn 很大的时候(大概大于100
2016沈阳站网络赛题解
acmer的随记
09-19 517
A 操作：在一个矩形内的每个格子中，第i个物品都加上ai个 求：某个矩形内每种物品的奇偶性 思路：奇偶性可以用异或来表示：第i个物品被加上奇数次可以表示为与(1<<i)(1<<i)的异或。所以，题意转化如下。 操作：矩形内每个格子异或x 求：矩形内异或和 （参照CF198 DIV1 D思路） 操作：分解成四个角落的点的异或 求：分解成四个[1,1,xi,yi]的四个矩形的【异或前缀】
2020年湖南中医药大学“华为杯”大学生程序设计竞赛——正式赛(赛后补题)
lunaol的博客
12-03 2847
目录A-幸福小组B-菜鸡驿站C-TC的火柴D-n盏灯E-电子手表F- TC的steam账号总结总结 A-幸福小组 题目描述 俗话说：男女搭配，干活不累。还据说，如果在做一件事情时男生和女生人数相等，学习和工作效率可以达到最佳【瞎说的】。 为了更好地学习计算机专业知识，某个班分为n个学习小组，如果某一个组中男生和女生人数相等，这样的小组称为“幸福小组”。 现在分别给出这n个小组中男生和女生人数，请你编写一个程序统计一共有多少个“幸福小组”。 输入 单组输入。第1行输入一个正整数n，表示学习小组数量。（n&lt
2020年第十一届蓝桥杯省赛试题及详解（Java本科A组）
future277809183的专栏
02-11 5070
结果填空 (满分5分) 结果填空 (满分5分) 结果填空 (满分10分) 结果填空 (满分10分) 结果填空 (满分15分) 程序设计（满分15分） 程序设计（满分20分） 程序设计（满分20分） 第一题: 门牌制作 【问题描述】 小蓝要为一条街的住户制作门牌号。 这条街一共有 2020 位住户，门牌号从 1 到 2020 编号。 小蓝制作门牌的方法是先制作 0 到 9 这几个数字字符，最后根据需要将字 符粘贴到门牌上，例如门牌 1017 需要依次粘贴字符 1、 0、 1、 7，即需要 1
ACM的那些故事
Geecky的博客
06-14 8364
我的ACM参赛故事 作者：PHPPLAY 来自豆瓣
ACM 回忆
Stormjing的博客
05-16 2591
关于acm 前言： 即将进入研二，ACM的事情也渐渐远去，记忆终将模糊，但那段奋斗永远让人热血沸腾。开个贴讲讲ACM与中南的故事，当是宣泄，也当是一种宣传吧。多年不写作文，也不知如今的文笔能否打动更多的有志青年为中南ACM梦想而奋斗。 用时间序好了，不定时更新。一层楼对应一个时期或一年，有补充就在那层楼的回复中。也许真写起来故事不长，就这样吧。 1947年 ACM(Association for ...
关于ACM，关于CSU
weixin_34122604的博客
10-12 817
原文地址：http://tieba.baidu.com/p/2432943599 前言： 即将进入研二，ACM的事情也渐渐远去，记忆终将模糊，但那段奋斗永远让人热血沸腾。开个贴讲讲ACM与中南的故事，当是宣泄，也当是一种宣传吧。多年不写作文，也不知如今的文笔能否打动更多的有志青年为中南ACM梦想而奋斗。 用时间序好了，不定时更新。一层楼对应一个时期或一年，有补充就在那层楼的中。也许...
对边bfs的最短路问题
weixin_34354945的博客
07-23 83
在我们的日常做题中，我们经常会碰到一些情况，他的信息都是依照边给出的，并要求我们在这些边的基础上进行bfs以此来找出某种最短路。由于对边进行操作，我们这里就需要用到邻接表。 先行分享一个邻接表的模板 #include<stdio.h> #include<algorithm> using namespace ...
ACM大学生程序设计竞赛在线题库最新精选题解_14651101.pdf
04-08
ACM赛题
竞赛资料源码-国际大学生程序设计竞赛中山大学内部选拔真题解（二）.zip
02-05
竞赛准备：适用于各种教育部认可的竞赛，如全国电子设计大赛、全国大学生智能汽车竞赛等，他们可以借助这些资料了解竞赛的规则、要求和技巧。 学习与项目开发：可以用作毕设、课设、作业和竞赛项目的开发基础，可以...
国际大学生程序设计竞赛例题解（一）.zip
11-08
程序设计竞赛相关代码、设计文档、使用说明，供学习参考 程序设计竞赛相关代码、设计文档、使用说明，供学习参考 程序设计竞赛相关代码、设计文档、使用说明，供学习参考 程序设计竞赛相关代码、设计文档、使用说明...
人偶师(Gauss消元xor版-多解)
nike0good |Oier&ACMer | 熟能生巧
03-09 1553
Problem 2 人偶师(alice.cpp/c/pas) 【题目描述】 n点m双向边的图，每个点有2个状态：开和关。每次操作改变一个点的状态，以及与其有边直接相连的点的状态。问开启所有点至少需要多少次操作。 【输入格式】 第一行2个整数n,m。 第二行n个整数，第i个数表示第i点的状态，0为关，1为开。 第3..m+2行，每行2个整数a,b，表示a和b直接相连，同一条边不会出现多次
BZOJ 1013([JSOI2008]球形空间产生器sphere-gauss消元练习)
nike0good |Oier&ACMer | 熟能生巧
04-21 1481
1013: [JSOI2008]球形空间产生器sphere Time Limit: 1 Sec  Memory Limit: 162 MB Submit: 1181  Solved: 654 [Submit][Status][Discuss] Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在，你被困在了这个n维球体中，你只知道球面上n+1个点的坐标
CMO 2012回馈(gauss整数解)
nike0good |Oier&ACMer | 熟能生巧
06-22 1470
本题是CMO(数学 Olympics) 2012 第二题 所以还是很坑的……（出题人是shuxuedi） 反正这题总算是写了一个远远长于正解的打表（Gauss消元-判断有无整数解/无解，已确定答案是否可行） #include #include #include #include #include #include #include #include #include using name
POJ 1222(Gauss消元xor版)
nike0good |Oier&ACMer | 熟能生巧
02-06 1455
EXTENDED LIGHTS OUT Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 5182   Accepted: 3403 Description Lights Out就是下图的游戏,给你一个5*6的矩阵.  你的目标是把灯全关上.
2024四川大学生程序设计竞赛题解
最新发布
07-18
2024年的四川大学生程序设计竞赛的具体题目和官方题解通常会在比赛结束后由主办方公布，或者参赛者可以在官方网站、论坛、博客等地方找到分享的内容。 如果你想了解往届的竞赛题解，通常会涉及算法、数据结构、数学...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值