首页 > 代码库 > HDU 4407
HDU 4407
这题不难吧,如果正在做组合的题。。。
使用容斥原理求解出(1~x)的与p互素的和,这是很容易的,很明显,首先要把p分解质因数。
而对于第二个操作,记录下他的转换的顺序,当要执行第一个操作时,遍历一次记录下的操作转换就可以了。
呃,这题虽然想到,但是,我的WA。看了网上的,思路和我的一样,我自己COPY别人的代码下来对拍数据,无误。。但就是过不了。能想到的边界数据也试过了,还是过不了。。求各路大神指错。
我的代码:
#include <iostream>#include <algorithm>#include <cstring>#include <cstdio>#define Np 700 #define LL __int64using namespace std;bool isprime[Np];LL prime[Np],np;LL dsorp[Np],dp;struct Update{ LL pos,val;};Update up[1500];LL upp;void initial(){ memset(isprime,true,sizeof(isprime)); np=0; isprime[1]=false; for(LL i=2;i<Np;i++){ if(isprime[i]){ prime[np++]=i; for(LL j=i*i;j<Np;j+=i){ isprime[j]=false; } } }}void dfs(LL i,LL num,LL tmp,LL &s,LL n){ if(i>=dp){ if(num==0){ s=(1+n)*n/2; } else{ LL sc=n/tmp; if(num&1){ s=s-((sc+1)*sc/2*tmp); } else{ s=s+((sc+1)*sc/2*tmp); } } return; } dfs(i+1,num,tmp,s,n); dfs(i+1,num+1,tmp*dsorp[i],s,n);}LL work(LL x,LL y,LL p){ dp=0; for(LL i=0;i<np&&prime[i]*prime[i]<=p;i++){ if(p%prime[i]==0){ dsorp[dp++]=prime[i]; while(p%prime[i]==0) p=p/prime[i]; } } if(p>1) dsorp[dp++]=p; LL ansy,ansx; ansy=ansx=0; dfs(0,0,1,ansy,y); dfs(0,0,1,ansx,x-1); return ansy-ansx;}LL gcd(LL a,LL b){ if(b==0) return a; return gcd(b,a%b);}bool coprime(LL a,LL p){ if(gcd(a,p)==1LL) return true; else return false;}int main(){ LL T,n,m,op,x,y,p,c; initial(); scanf("%I64d",&T); while(T--){ scanf("%I64d%I64d",&n,&m); upp=0; for(LL i=1;i<=m;i++){ scanf("%I64d",&op); if(op==1LL){ scanf("%I64d%I64d%I64d",&x,&y,&p); LL ans=work(x,y,p); for(LL k=0;k<upp;k++){ if(up[k].pos>=x&&up[k].pos<=y){ bool flag1=coprime(up[k].pos,p); bool flag2=coprime(up[k].val,p); if(flag1){ if(flag2){ ans=ans-up[k].pos+up[k].val; } else ans=ans-up[k].pos; } else{ if(flag2) ans=ans+up[k].val; } } } printf("%I64d\n",ans); } else{ scanf("%I64d%I64d",&x,&c); LL i; for(i=0;i<upp;i++){ if(up[i].pos==x){ up[i].val=c; } } if(i==upp){ up[upp].pos=x; up[upp].val=c; upp++; } } } } return 0;}
别人的代码:http://blog.csdn.net/dgq8211/article/details/8031257
#include <map>#include <stdio.h>typedef __int64 LL;using namespace std;#define N 650bool np[N];int prime[120],fac[9],F_top,p;LL ans;void get_prime(){ int top = -1; for(int i=2;i<N;i++) if(!np[i]) { prime[++top] = i; for(int j=i+i;j<N;j+=i) np[j] = true; }}void Div(int n){ F_top = 0; for(int i=0;prime[i]*prime[i]<=n;i++) { if(n % prime[i]) continue; while(n % prime[i] == 0) n /= prime[i]; fac[F_top++] = prime[i]; } if(n != 1) fac[F_top++] = n;}LL PreSum(int n){ return (LL)n*(n+1)/2;}void dfs(int i,int cnt,int m,int n,int num,int x) // C(n,m).{ if(cnt == m) { LL sum = num * PreSum(x/num); m&1 ? ans -= sum : ans += sum; return ; } if(num*fac[i] > x || n-i < m-cnt) return ; dfs(i+1,cnt+1,m,n,num*fac[i],x); dfs(i+1,cnt,m,n,num,x);}LL sovle(int n){ ans = PreSum(n); for(int m=1;m<=F_top;m++) dfs(0,0,m,F_top,1,n); return ans;}int gcd(int a,int b){ return b ? gcd(b,a%b) : a;}int main(){ int z,n,Q,cmd,a,b; get_prime(); scanf("%d",&z); map<int,int> M; while(z--) { M.clear(); scanf("%d%d",&n,&Q); while(Q--) { scanf("%d",&cmd); if(cmd == 1) { scanf("%d%d%d",&a,&b,&p); Div(p); ans = sovle(b) - sovle(a-1); for(map<int,int>::iterator it=M.begin();it!=M.end();it++) if((*it).first >= a && (*it).first <= b) { if(gcd((*it).first,p) == 1) ans -= (*it).first; if(gcd((*it).second,p) == 1) ans += (*it).second; } printf("%I64d\n",ans); } else { scanf("%d%d",&a,&b); M[a] = b; } } } return 0;}
HDU 4407
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。